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ABSTRACT 


A general modelling framework was suggested by Min 
and Ray [l7 ] for emulsion polymerization reactions, A simplified 
model based on the framework has been proposed to simulate 
the emulsion polyiaerization reaction in a continuous stirred 
tank reactor. The model uses the population balance equa,tion 
approach. The method of moments has been used to solve the 
population balance equations. The model is capable of 
calculating conversion, rate of polyLiorizat ion, number of 
particles, particle size distribution and the molecular 
weight distribution of resultant polsaner latex. The variables 
studied were residence time and initiator concentration. 

The results obtained using this model were compared with 
the experimental data and the theoretical results of DeG-raff 
and Poehlein [21]. The proposed model is found to predict 
the aforesaid quantities well at low residence times. The 
predictions of the model are in moderate agreement with 
data at high residence times. The theoretical approach by 
DeGraff and Poehlein behaves in the opposite way. 



CHAPTER 1 


IHTRODUCTIOH 

The 50 year old technique of emulsion polymerization 
is still being widely used in polymer industry* It often 
has significant advantages over homogeneous polymerization. 
Its final product is a latex and hence p- ssesses the obvious' 
•advantage of direct application in cases like paints and 
coatings* Apart from this, there are distinct advantages 
in operating the reactor. The viscosity of a polymer latex 
in an emulsion polymerization reactor is low and hence the 
temperature control is relatively easy. In additioni a 
considerably high molecular weight product can be obtained 
at much higher reaction rates. For these reasons, emulsion 
polymerization has become commercially important in polymer 
industry. At present it is a major processing step in the 

manufacture of polymer products such as ^nthetic rubbers, 

•<* 

paints, inks, coatings, adhesives, high impact strength 
co-polymers etc. 

The commercial technique of emulsion polymerization 
can be divided into two classes: (1) batch and (2) continuous 

Batch systems are relatively simple and were developed 
first. These systems, althou^ less efficient, are used 
for low production rates and for speciality applications. 

The advantages of this technique are the relatively simple 
equipment to be handled and the fact that a small amount 
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of material can be easily processed. This type of poly- 
merization has been extensively studied in literature, 

ITumerous articles are available in literature regarding the 
experimental and theoretical aspects of emulsion polymerization 
in batch. The more important of these works are by anith 
and Ewart [IJ, Stockmayer [2], O’Toole [3], Garden [4-‘10], 
Ugelstad [11], Medvedev [12], Min and Ray [15] • These and 
many other works have been comprehensively reviewed by 
iil'efxander and Nappel [14], Vander Hoff [15, 16], IViih and Ray [17]. 

As the annual output of the polymeric materials 
increases, more and more attention is being given to the 
use of continuous reactors. Although continuous stirred 
tank reactor (CSTR) is widely used commercially for emulsion 
polymerization, very little published information is available 
about those reactors. Experimental data is limited since 
the equipments are more costly and less flexible, experimental 
run times are longer and larger amounts of raw materials 
are consumed. Similarly the modelling of these reactors 
involves some conceptual problems which are not present in 
batch systems. 

The earliest of the models developed for emulsion 
polymerization in a CSTR is by Greshberg and LongfielE. [18], 
They have employed a residence time distribution model coupled 
with the Smith-Bwart kinetics. Omi and co-workers [19] 
and Momura and co-workers [20] have employed Smith— Ewart 
kinetics for different reaction systems* They have also used 
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different models for initiation of particles. DeG-raff 
and Poehlein [21] have employed Stochnayer^o [2] theory 
coupled with residence time distribution for modelling a 
single stirred tank reactor. Recently population balance 
approach has been employed for modelling of continuous 
reactors, 

Stevens and Funderburk [22] have employed a simple 
population balance approach to model a CSTR, Further work 
using this approach was done by Thompson and Stevens [23] > 
Thompsonand co-workers [24] and Cauley and Thompson [ 25] , 

These authors have usod population balance equations to 
generate radical number distribution. They have used size, 
independent particle growth rate, Cauley and Thompson [25] 
have further assumed instant aneous termination of radicals 
inside the particles. They have not developed material 
balance equations for various reactive species. The conversion 
and molec\iLar weights of the polymer could not be calculated, 

Min and Ray [17] have developed a comprehensive 
model for emulsion polymerization in batch and CSTR, They 
have taken into consideration a number of mechanisms like 
desorption of radicals from particles, chain transfer, 
initiation in aqueous phase, end -group stabilization of 
polymer particles, initiation by precipitation etc. They 
have, however, not presented any numerical result from 
this comprehensive model for OSTR. 
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The latest work on modelling of CSTR is by Eirilov 
and Ray [26], They have used the general model of Min and 
Ray [17] for a specific case by making a number of simplifi- 
cations. They have employed simple population balance 
equation for particle size distribution and have used 
Stocki-layer * s [2] relation for calculating the average number 
of radicals per particle. They have calculated average 
particle size and conversion obtained in the reactor. They 
did not calculate either the particle size distribution or 
the molecular weights of the polymer product. 

In the present work, an attempt has been made to develop 
and solve modelling equations for emulsion polymerization 
in a CSTR, The population balance approach has been used. 

The popiilation balance equations employed are more oompreheiji- 
sive. As a result, the molecular weight distribution of 
the resultant product can also be generated in addition to 
conversion and the particle size distribution. 



CHAPTER 2 


?M-§I£4£L.lPP.^-.AIg_.. .P0PIILATI0N BALANCE EQUATIONS 

Before starting the actual formulation of modelling 
equations, it is necessary to tinder stand the physical picture 
and details of various mechanisms operating in an emtilsion 
polymerization reactor. The currently accepted [14-17] 
picture is as follows; 

A typical emulsion polymerization recipe consists of 
100 parts of water, 50 parts of monomer, 2 parts of emulsifier, 
0,1 part of initiator and small amounts of chain transfer 
agents and other additives. The soap is present in micelles, 
on the surface of monomer droplets and in the aqueous phase. 

The monomer is present mainly as droplets and a small amouni? 

(1 per cent) is solubilized in the micellar soap, A small 
portion of the monomer may also be dissolved in the aqueous 
phase. Generally, the initiator is present in the aqueous 
phase* In early stages of the reaction, the system consists 
of 3 types of pai’ticles dispersed in continuous aqueous 
phase* (l) emulsified monomer droplets of size 10000 A and 
at a concentration of 10 /ml, (2) soap micelles with size 50- 

Q -- TO 

100 A and at a concentration of 10 /ml and (3) emulsifier 
stabilized monomer swollen polymer particles of size 500- 
1000 A, and at a concentration of 10^^-10^^/ml in later 


stages 
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As the reaction proceeds the polymer particles grow 

'I ' 

in size, consimiing the monomer present in the whole system, 
thus requiring more and more amount of soap to stabilize 
them, The soap is supplied by micelles which decrease in 
number and finally disappear. Similarly, the monomer 
droplets shrink in size as they supply monomer to the polymer 
particles and finally disappear, leaving the reaction 
system consisting of polymer particles dispersed in the 
aqueous phase , 

Now some of the mechanisms operating in the reaction 
mixture will be discussed, 

Nucleation of polymer particles is one of the least 
understood phenomena. It can occut in two different ways. 

The initiator present in the aqueous phase thermally 
decomposes into constituent radicals- (reaction temperature 
approximately 50--70°C) and thes0 radicals grow by adding 
the monomer units from the monomer dissolved in the aqueous 
phase. During this process the radicals and the active 
polymer chains exhibit Brownian motion. During their motion 
they can collide with micelles, get absorbed into them and 
start the polymerization reaction thereby convei’ting the 
micelles into polymer particles. They can also collide with 
and get absorbed in existing polymer particles. The radicals 
which are not captured by micelles or particles continue to 
add on monomer units arid in the process decrease their 
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aCLueous phase soluhility. When a critical limit of chain 
size is reached they will precipitate and initiate new particles* 
Both these mechanisms of particle nucleation operate simul- 
taneously* The former is more important for monomer scarcely 
soluble in water, while the latter is predominant for 
monomers which are relatively more soluble in water. 

The particles thus initiated consist of monomer, 
polymer chains and the live polymer radicals. Polymerization 
continues in the particles, the monomer being supplied by 
monomer droplets. Radicals absorb into and desorb from the 
particles. Propagation, termination as well as chain transfer 
can take place inside particles. Soap required for stabili- 
zation is supx'ilied by the micelles until they disappear. 

As more and more amount of monomer is consumed the monomer 
droplets decrease in size and number. Ultimately they 
disappear leaving only the polymer particles containing live 
radicals and dissolved monomer. As the reaction proceeds to 
higher range of conversion the monomer concentration in a 
particle decreases, viscosity of reaction medium increases and 
the reaction can become diffusion controlled. Gel effect 
may play an important role at higher conversion. Apart from 
these mechanisms, particles can also coagulate with each 
other to a limited extent. The extent of coagulation depends 
upon the particle size, amount of soap adsorbed on each 
particle, reaction tempera'ture etc. 
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The polsnner particles consist of monomer, polymer 
chains and live radicals as stated ear3.ier. However, the 
exact structure of these particles is not yet known. The 
monomer and polymer may solubilize each other to form a 
homogeneous mixture. If that is the case, reaction will 
occur at a constant rate throughout the homogeneous particle. 

On the other hand, ¥il3.iams [27] has postulated a non-uniform 
morx)hology for the particles, comprising of a polymer rich 
core and a monomer rich shell. The reaction is assumed to 
take place on the growing monomer polsmier interface, Faper[28] 
has suggested exactly the opposite structure consisting of 
a monomer with core and a polsrmer rich shell. 

This, in short, is a physical picture of the emulsion 
polymerization process. The aforesaid mechanisms operate 
simultaneously in a batch or continuo^is emiilsion polymerization 
reactor. The objeQtive of the present work is to use the 
knowledge of these basic mechanisms to model the character- 
istics of a xjolymer latex produced in a continuous stirred 
tank reactor. 

The problem can be approached in two different ways: 

In the residence time distribution approach the 
segregated flow model of a OSTR is employed. In this model 
the reactor is assumed to consist of a number of small 
independent reaction units, Each reaction unit is assumed 
to reside for different times inside the reactor. The reaction 



units ai-e ass-umed to operate as isolated batch reactors, 
ilhe residence time distribution of different units combined 
with the knowledge of batch behaviour is believed to represent 
the CSTR characteristics. This method was most recently 
used by De&raff and Poehlein [21] for the modelling of 
emulsion polymerization CSTR, Here a particle was identified 
as a reaction unit. They assumed the usual exponential 
residence time distribution. Another approach for modelling 
is the use of population balance equations which will be 
employed in this work. A short discussion about the theory 
of population balance equations is available in ref. 29. 

A brief introduction to this method is given below, 

T'or any particulate system, a particle phase j^ace 
consisting of a minimum number, say p, of independent co- 
ordinates, attached to a particle distribution, which 
completely describes the properties of the distribution, is 
postulated. The particle phase space may be divided into 
two subregions consisting of m internal and three external 
co-ordinates. Hx-ternal co-ordinates refer to the three 
spatial co-ordinates of the particle. Internal co-ordinates 
quantitatively represent the state of the particle independent 
of its spatial position. 

let us consider a (m+3) dimensional particle distri- 
bution function n(Z,t) defined over region Z consisting of 
three external co-ordinates and m independent internal 
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co-ordinates. In a dynamic particulate ^stem, individual 
particles change their position in the particle phase 
space. If the change of particle co-ordinates is gradual 
and Continuous in any direction the rate of such 'motion' 
is called as convective particle velocity along that 
co--ordinate axis. Thus the particle phase-space velocity 
vector is defined as 


'V'ssV 4* 


^int + ^ext 


• • V 6 + v^ 6^ + v^6_ + 

mm 1 1 2 2 3 3 


(2 , 1 ) 


where v's are the components of the pai-ticle velocity 

and 5's are the unit vectors along the various co-ordinate 

axes. In general., the external velocity "v^can be calculated 

ext 

from the fluid velocity and momentum exchange relationship 
between fluid and particle. The internal particle velocity 
"^nt assumed to be unique function of state of the particle, 
There is also a possibility of particles suddenly 
appearing or disappearing in particle phase space along any 
co-ordinate. This possibility is represented by the birth 
and death functions of particle distribution. In general, 
these birth and death fLinctions can be related to the state 
of the total system and the position of a particle in the 
particle phase ^ace. 

The population balance equation for particles in some 
fixed subregion of particle phase space will be based on 
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the usual conservation law. 

Accuoaulation = Input - Output + net generation 

( 2 . 2 ) 

Consider the subregion to move convectively with 
particle phase space velocity v. The population balance 
equation can be written as 


n dZ = 


(B-D) dZ 


( 2 . 3 ) 


where B and D represent birth and death functions. 
Expanding L.H. S, 


ndZ = 


dZ + n 


a?. I 


^ H 'Y' 

[“^ + • (-g:^ n)]dZ 


where x are the co-ordinates of system. 

Replacing by (v^^^ + and combining with 

2.3 gives the following expressions 

) ^ ^ ' ^"^int*^^ s= 0 


The equation is valid for any arbitrary region 2^. 
Hence the integrand vanishes for all Z. 

Therefore + g, ("^ext’^^ ^ ^ 

( 2 . 4 ) 


This is a general expression which can be used to 
describe any particulate system, provided the number of 
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particles in tlie system is siiCficiently large so that 
n(Z,t) can retain its statistical significance. 

This eq,uation will now be applied to a continuous 
well stirred particulate system of volume l/‘(t) , having a 
population density function n(Z,t)> p iniDut and q. output 
streams. Since in a well stirred system, tracing of external 
coordinates of particles is not of interest, this equation 
can be integrated over the range of external co-ordinates V or 
VO lime occupied by the system 

( (:^ H- D-B = ° 

V 

Therefore + V '^nt*^ > ^'^xt'^ dV = 0 

^ (2,5) 

The integral can be transformed into a surface 
integral of population flux over the surfaces of the system. 
The surfaces are: the cross-sections of input-output streamg, 
free interface of liquid and to tad. fluid -particle interface. 
The final expression is 


[9 V 


ext 


•n dV =s n 


Z Q.m + T 


n 


(2.6) 


P <1 

where Q*s are the flow rates of input-output streams. 


Combining 2.5 and 2.6 and simplifying 


^ + vw: 


int 


•n + n 


= b-d + 1 - z: 


On 


(2.7) 
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This expression represents the macroscopic population 
balance equ.ation for a particiilate , well mixed system. With 
the knowledge of input-output streams, expressions for 
velocity v^^ , birth and death functions this equation can 
be used to model emulsion polymerization in a CSTR mathe- 
matically, The values of population density n(t) can be 
obtained and used to calculate the characteristics of the 
resultant polymer latex. 

In next chapter, this equation will be used to develop 
the design equations for an emulsion iDolymerization CSTR, 



CHAPTER 3 


PO MUIATIOCT OP.MOPELL IHG EQHATIOFS 


Prom the previous chapters it can be seen that an 
emulsion polymerization latex is a particulate system 
comprising of a continuous aq.u.eous phase and a particulate 
phase consisting of polymer particles, micelles and monomer 
droplet a, A set of equations, consisting of ma,terial balance 
equations for continuous phase and population balance 
equations for particuHate phase, completes the mathematical 
representation of physical pol 3 merization sj/stem. Energy 
balance equations for the system are generally unnecessary 
because of low latex viscosity and high heat transfer rates, 
which result in an isothermal operation. 

In the following pages development of the modelling 
equations will be given. This presentation closely follows 
the approach given by Min and Ray [17]. Per convenience, 
nomenclature for this chapter is provided at the end of the 
chapter. 


ginetie Schenie t The usual kinetic scheme involving initiation, 
propagation, termination as well as chain transfer to monomer 
and transfer agents has been used. The schene is given 
below: 


kj 


1, Initiation ; [l*^ 


w 


ki 




R+R — ^l^^l 


R+R = 
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24 Prflja^atioQ^ 

3. Tern illation : 


k. 


P +M 
n 


P + P 
n m 


P + P 
n m 


to 


k 


td 


-> 


4* Chain transfer: 



+ P. 


n 



+ P.. 


n 


^n+1 

^n+m 


M 4-M 
n m 


^ft 


k. 


-> ^V^l 


■fm 


-^M + P. 

^ n 1 


Particulate Phase : 

For characterizing a polymer latex, the polymer particles 
are of prime importance. Very little, if ai:iy, reaction takes 
place in the monomer droplets or micelles. Hence it is not 
necessaryto develop separate population balance equations 
for these species. It will be assxmied that these distributions 
can be represented by average micelle size (volume) v^^ and 
average droplet size (volume) v^ respectively. 

With the polymer recipe and average residence time 
of the reactor usually employed in practice, the number 
density of polymer particles is of the order of 10^‘^/cm^ 
emulsion. Hence we can attach an appropriate distribution 
function to the particles and use macroscopic population 
balance equation as given in the last chapter. 

^ + IT. “ If(logV) = B-D + 21 -T^ 

(2.7) 

Consider a reactor of volume V, with an input and 
output stream, both having a volumetric flow rate of q,. The 
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The contents of reactor are a s sinned to he completely 
mixed. Let r^^ represent the net generation term (B-D) 
resulting in discrete changes of variables due to the 
various mechanisms* Let n be the distribution function 
under consideration. The aforesaid equation, when applied 
to such eystem becomes 


^ n 

■p 


^int 


m 


+ V-'v',. „4.*n 




+ i n 


feed 


- n 

g U 


(3.1) 


where Q represents the mean residence time and m 
depends upon the number of mechanisms to be included in the 


model. 


Radical llumber Distribution : 

Consider a distribution function f(i,v,t) which 

represents number of particles per cm^ of emulsion^ having 

volume V to v+dv and i radicals, at any time t^ The knowledge 

of this distribution is necessary to calculate the rate of 

polymerization and The population balance equation 

avg 

corresponding to this function is 

m 

« I f(i,v,t} (3,2) 

Pirst term represents rate of change of f(i,v,t) with 
respect to time, r^ in the second teim represents the rate 
of change of volume of particle with respect to time,[ie . 



f(i,v,t) 


’eed 
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The size of a particle "varies due to the propagation 
reaction taking place inside the pai’ticle and due to particle 
swelling. This quantity will depend upon the model assumed 
for particle structure, For homogeneous particle morphology 
Garden [4] has calculated -this quantity hy correlating the 
propagation mechanism with 'the monomer concentration and 
polymer and monomer densities. He obtained following 
relationship for the conditions when the monomer droplets 
are present in the system and the aqueous phase is saturated 
(with monomer) . 


dv 


= r. 


k-n . . d 

P m 


V 




.i 


A 


(1-9) 


(3.3) 


The rate terms r^ in the BHS of 3.2 depend upon the 
particular mechanisms to be taken into consideration. These 
terms will be discussed below. 

First "berm, r^^ represents the change in f (i,v,t) due 
to absorption of radicals into the particles. Radicals can 
either propagate and precipitate as oligomers in the aqueous 
phase or can be absorbed into the particles. The relative 
importance of these two mechanisms depends upon the monomer 
solubility and the propagation rate constant in the aqueous 
phase, Gardon [4] has obtained an expression for rate of 
radical absorption by computing the probability of collision 
of a radical with micelles and particles. The final 
expression is 
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R- 

la 


2 f k, [1] 

_ _<i. ■‘•w 

\ + Ap 

_ 2 f l^d iilw 
~JV/D' 


when S >> 4/L 
when S /C.C.Aj'L 


where 1 represents the average length travelled, by a growing 
chain before it reaches the critical length at which it 
precipitates and is given by 


l> ^ (2D t 


prec'^ 


0.5 


Here represents the time required by radical to reach 

critical chain length for precipitation and is a function of 
the propagation rate constant and monomer concentration in 
the aqueous phase. D is the diffusion coefficient of 


radicals in the aqueous phase. 

Rate rp in eqn.3.2 is then given by the following 
expression 

^1 ~ \a ^ ^A - f(i,v,t)) (3.4) 

where the first term represents the formation of f (i^V|t) 
by absorption of a radical into f(i-l,v,t) type of particles 
and the second term represents destruction of f(i,v,t) type 
of particle by absorption of a radical, thus transforming 
it into a f(ifl,v,t) type particle. 

Second rate term (rp) represents the effect of formation 
of a radical pair inside the particle by the decomposition 
of monomer-soluble initiator. 
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^2 " ^ ‘ (f(i+2*v,t) -f(i,v,t)) (3*5) 

Radicals can also desorb from the particle* This 
effect can be accounted for by the following expression; 

k ap - 

^3 = [ (i+l) f(i+l,T,t) - i f(i,v,t)] (3.6) 

I'irst term represents the formation of f(i,v,t) type 
particle by radical desorption from a f(i+l,v,t) type particle 
and second term represents the destruction of f(i,v,t) by 
desorption of a radical. 

Rate of change of f(i,v,t) by termination of two 
radicals inside the particle, is given by the following 
expression; 

^tc 1 ” 

^4 - 2 %“v^ [(i+2)(i+l) f(i+2,v,t) -(i) (i~l)f (i,v,t)i 

-A. 

( 3 . 7 ) 

The rate of formation of .new particles by radical 
absorption in micelles is given by 

"^5 = ^ ( 3 . 8 ) 

6 functions on ri^it hand side of the expression 
indicate that new particles, thus formed, will have volume 
v^ and will contain one radical. 

Particle Size Distribution ; 

let P(v,t) be the particle size distribution function. 

This distribution is quite important since the particle size 
is one of the characterizing properties of a polymer latex. 
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Consider the balance equation for ^(v^t), wliich represents 
the number of iDarticles having volume between v and v+dv at 
any time t. It is related to radical number distribution 
by the following equation: 


I'(v,t) = TZ f(i,v,t) (3.9) 

i=0 

By combining 3.2 to 3.9 ibe population balance 
equation for F(v,t)caa be ^derived to be 

I''(y.t) 




+ (r^ I‘'(v,t)) 


Rate of change of 
l'(v,t) with respect 
to time 


Rate of change of P(v,t) due to 
+ volume growth of particles 


R 


fir \ 

i' V 

Rate of change of I‘(v,t) 
duo to generation of new 
particles of size v 


+ 


© 






feed 


-B(v,t)] 


m 


Change of I'‘(v, t) due to 
+ input and output 

(3.10) 


where 


k.^ .. d .. - 

„ P m <p 

"v = Tq • ^ 


"avg 


(3.11) 


and [i.G. average number of radicals per particle] is 


given by 


V 


( r i f(i,v,t) d? 


m i^O 


avg 


,3^ OC 


( 3 . 12 ) 


I ^21 f(i,v,t) d v 


v„ i=0 
m 
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Live Polyme r Chai n Lengbli Distribtit i.on : 

In order to calciHate the molecular weight distribution 
the chain length distribution of the polymer radicals present 
inside the particle is needed. This distribution is used to 
solve the dead polymer chain length distribu.tion which is in 
turn used to obtain the m.w.d. of the polymer product. 

Consider the particles having volume between v and v+dv and i 
radicals at any time t. These i radicals will in general have 
different chain lengths. Amongst these radicals consider the 
once having chain length n. Let L(n,i,v,t) represent the 
total number of radicals of chain length n present in all the 
particles having volume between v and v+dv and having i radicals, 
at any time t. The population balance for such radicals is given 


by the following expression. 


,i,v,t) 1 


+ I [L(n 


feed 


2 — 

i=l 

l(n,i,v,t)] 

( 3 . 13 ) 


lirst term represents change in I(n,i,v,t) with respect 
to time. The second term represents change in L(n,i,v,t) due 
to change in size of a particle, where r^ represents the 
rate of change of volume with respect to tine. 

The rate terms r^ on right hand side, as before, 
depend on specific assumptions made. Firstly consider the 
effect of absorption of radicals. It has been assumed that 
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only primary radicals can absorb into the particle and 
that the particle nucleation by precipitation in aqueous phase 
has been negloigible. Let g(n,i^v*t) represent number of 
radicals of .chain length n in a f(i,v,t) type particle. 

When a radicalis absorbed in a f(i-l,v,t) type particle it 
converts tho particle into a f(i,v,t) type, making in the 
process, a contribution of g(n,i~l,v,t) radicals to L (n, i,v, t) . 
Similarly when a primary radical is absorbed in a f(i,v,t) 
type particle g(n,i,v,t) nuiaber of radicals will be lost 
from l(n,i,v,t). When a primary radical is absorbed in a 
particle it adds only to L(l,i,v,t) . Thus tho total effect, 
per cm^ of emulsion can be represented as; 




•g(n,i-l,v,t) 



g(n,i,v,t) 


Ap 

f(i~l,v,t) 

^^1 


- - f(n,i,v,t) +f (i-l,v, t) 5nj 

■ ' (3.14) 

From the definition 

L(n,i,v,t) , N 

g(n,i,T,t) (5.15) 

The initiator dissolved in a polymer particle can 
decompose to give a pair of radicals of chain length one. 
each, A contribution of g(n,l,-2,v,t) is made to I(n,i,v,t) 
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per generation of a pear due to such an event. In addition 
to this tho radicals foraod themselves contribute to 
addition of radicals with chain length one. The same process 
occurring in a f(i,v,t) type particle wo\iLd result in a 
loss to l(n,i,v,t) 

^2 ~ ^ipf^^^P [ (s(^»i“-2,v,t) + 26n^)f (i“2,v,t) 

-g(n,i,v,t)* f(i,v,t)] 

= [ I]p‘'^*Wj^[Ii(n,i-2,v,t) - L(n,i,v,t) +26n^f (i-2,v,t) ] 

(3.16) 

The effect of radical desorption from particles is 
given by 

^2 = ICq —[ (i+l)f (i+l,v,t) •g(n,i+l,v,t) (1- 
- i f(i,v,t) g(n,i,v,t)] 

= -~[(i) L(n,i+l,v,t) - (i) L(n, i,v,t)] (3.17) 

Here the f irst term represents the production of L(n,i,v,t) 
by desorption of radicals from f(i+l,v,t) except when a 
radical of length n desorbs. The , second term represents loss 
due to desorption from radicals from f(i,v,t}type particles. 

The effect of the kinetic steps inside the polymer 
particle depends upon the structure of particle. Assuming 
homogeneous particle moxphology, the effect of propogation 
and transfer reaction can be expressed as 
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= ^p[^]p f •g(n-l,ijVjt) - kp[M]p f (i,v, t) g(3a*i|V^ t) 

** kp^[Tr]p g(n,i,v^t) + kp^[Tr]p«i f (i>v,t)6np 

-^fJ^'IJp g(n,i,v^t) + kp^[M]p i f(i,v,t) Sn^ 

*= i*^p[J^]p[l'(n~l,i,v,t) - L(ii,i,v,t)] +k^^[ Tr] p[ ~L(n, i,v, t) 

+f (i,v,t) •i.6np] +k^[M]p[f (i,v,t)ri*6np-L(n,i,v,t)] 

( 3 . 18 ) 


Iho first two terms represent the effect of propagation; 
third and fourth term represent chain transfer to transfer 
agent and the fifth and sixth term represent chain transfer 
to monomer. The chain transfer reactions involve addition 
of one chain to f(l,i,v,t). This effect has been accounted 
for in fourth and sixth term. 

The effect of termination. reactions can be written as 


- (i+1)) •L(n,i+2#v,t) 


•5 Lv 2 

iiliizil L(n,i,v*t)] 

i(n,i,v,t)] 


“ “ 1 ^ 


(3.19) 


first and third terms in the initial expression 
represent change due to termination between any two radicals. 
Third term represents the termination between a radical of 
chain length n and remaining (i+l) radicals. The possibility 



25 


of terniination between two radicals of sane cliain length, n, 
which will be a second order tei-m, has been neglected. 

The effect due to radical absorption in micelles is 
given by the following expression; 

= % ( 2 ^) 

•*>r ill 

Dead Polymer Chain Length Distribution; 

This distribution is related to the chain length of 


dead polymer molecules inside the particles and is thus 
closely connected with m.w.d. of the resultant polymer latex. 
Consider the particles having volume between v and v+dv and 
possessing i active radicals at any time t. Apart from these 
i active radicals, each particle contg,ins some dead polymer 
chains. Let G(n,i,v,t) represent total number of these dead 
chains contained in all f (i,v,t) type particles at any time t. 
The population balance equation for this function can be 
written as 

-it vv a 

+ ^ rG'(n,i,v,t) \. “G(n,i,v,t)] (3.2l) 

. i feed 

First term represents the change in Gr(n, i,v,t) with 
respect to time and the second term expresses the change due 
to change in volume of particle. 
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Rate terms on tlie ri^it hand side depend upon the 
particular mechanisms to he considered in the model. Deri- 
vation of these terms is similar to that for corresponding 
terms in live polymer chain length distribution. These 
terms will be discussed briefly. 

The first rate term r^ represents the effect of radical 
absorption on G(n|i,v,t) . Lt g*(n,i,v,t) represent number 
of dead polymer chains per pai’ticle* Hence G(n,i,v,t) = 
g‘(n,i,v,t) f (l,v,t) , When a new radical enters a polymer 
particle of type f(i-l,v,t), it converts that particle into 
f(i,v,t) type particle and a contribution of g' (nji-^l^v^t) is 
made towards G(n,i,v,t), Hence the total contribution of this 
mechanism towards change in G(n,i,v,t) is given by 

f(i,v,t) g*(n,i,v*t) 

5* G(n— 1, i»v,t) — G(n,i,v,t)] (3*22) 

Similarly the effect of radical initiation in particles 
and desorption from part ides can be shown to be given by 
following expression. It is assmed here that dead 
polymer chains can not desorb from particles. 

Tg « fj^Ttllp [G(n,i-2,v,t) - G(n,i,v,t)] ( 3 . 23 ) 

y » G(n,i-f-l<v,t) - i G(n#i,v,t)] 

^ ( 3 . 24 ) 
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'Ihe effect of the kinetic steps is more important in 
case of G(n,i,v,t) and can be given by following expression: 

^4 ~ ^ft^^^^P P ^(njfiyVjt) 

h, +k. , - 

+ 2 -^^ [ (i+2) (i+1) G(n,i+2,v,t) - (i) (i-l) G(n,i,v,t)3 

k, - n-1 

+ [f(i+2,v,t) ^2 g(ni,i+2,v,t) g(n-m,i+2,v,t ) ] 

k - 

+ 2“^ [g(n-»i+2,v,t) f(i+2,v,t)] (3.25) 

In this expression the first two terms represent 
formation of dead polymer chains due to chain transfer with 
transfer agent and with monomer. Third term expresses the 
change in G(n, i,v,t) by termination of any two radicals in 
f(i+2,v,t) type particle (i.e. formation) and in f(i,v,t) type 
particle (destruction), Pourt^i term represents the additional 
contribution when two terminating radicals form a polymer 
chain of lengfh n. Last term represent the possibility of 
disproportionation involving a live polymer radical of chain 
length n. 

Apart from these population balance equations, balances 
which define the environment of particles must also be written. 
These equations include material balances for various reactive, 
species and balances for number of droplets, micelles etc. 

Em ulsifier Balances: 

The total amount of soap present per cm^ emifLsion is 
given by following expression: 
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a(s)„ 1 -- It - 

=i[[S]T (3.26) 

This soap may only be present in aq.ueous phase as 
micelles or may be aboohbed on the particles and droplets. 
The process of transfer of emulsifier will be assumed to 
be extremely fast. Under these conditions the balance for 
number of micelles yields 
pi ^ V 

m (3.27) 

The terms on I’ight hand side represent total soap, the 
soap present in the aqueous phase, the soap adsorbed on 
particles and the soap adsorbed on the monomer droplet 
respectively. 


Total Monomer Balance; 


The balance for the total amount of monomer present in 


;he reactor is given below; 


d Vp[M]„ Vp- - 

~aT^= e-COTi 




I - 

^ i f (i,v,t> 


(3.28) 


o i=0 


The first two terms represent the effect due to input 
and output terms. The third term represents the loss of 
monomer due to propagation in particles. Here the small 
contributions due to chain transfer and propagation in 
aqueous phase are neglected. 
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Me nomer Droplets ijSalance ; 

Monomer droplets balance can be written as follows. 

As in the previous case the dynamic of transfer of monomers 
will be assumed to be extremely fast, 

- \oX “• \ol (3.29) 

The terms on fight hand side of equation represent 
total monomer content, monomer in particles, and monomer in 
micelles. 

Aqueous Phase Balances: 

w# w ir-tr . - ,;ii>i..iiiiijr- 

The balance for concentration of the initiator in 
aqueous phase is given below 

^ - V -- I - 

C3.30> 

Lastly an e:35>ression for the monomer concentration 
in polymer (i.o, [M]p or * 9 ) is required. The nature of the 
equation will, in general, depend upon the model for particle 
structure considered, for a homogeneous particle morphology 
the free energy of mixing of polymer can be equated to the 

surface free- energy [ 4 ]. The resultant expression for 9 is: 

. 2 Y V T 

( 1 - 9 ) ln( 9 ) + 'f (1«-^)^ « (5.51) 

So far the population balance equations for particulate 
phase and the material balance for aqueous phase have been 
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discussed. These balances are assumed to represent the 
given system completely. These equations can be applied 
to a practical system and the results can be used to obtain 
desired characteristics of the resultant polymer latex. 
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HOMEFCLATUBE 

surface coverage by the emulsifier, cm^/mole 
surface area of micelles, cm^ 
surface area of particle of size vy cm® 
total surface area of raicelles, cm®/cm^ emulsion 
total stirface area o'f particles, cm®/cm^ emulsion 
density of polymer,, gms/can^ 
density of monomer, gms/cm^ 

efficiency of initiator decomposition, dimensionless 

number of particles containing i radicals and possessing 
volume between v to v+dv at time t^ molos/cm^ emulsion 
t'(v,t) number of particles having volume between v and v+dv 
at time t^ moles/cm^ emulsion 

g(n, i,v,t) number of live chains of chain length n present 

in a particle containing i radicals and possessing 
volume between v to v+dv at time t, dimensionless 
g* (n,i,v,t) number of dead chains of chain length n present in 

particle containing i radicals and possessing volume 
between v to v+dv at time t, dimensionless 
©(njijVjt) number of dead chains of chain length n present 

in all the particles having i radicals and possessing 
volume V to v+dv, at time t, moles/eo^ emulsion 
i number of radicals in a pai’ticle, dimensionless 

^aTg average number of radicals per particle, dimensionless 

[1]’, initiator concentration in aqueous phase, mol es/cm ^ , water 

It ^ ■ 

[I]p initiator concentration in polymer particles, moles/cm^ 


GE 

ap 

Aj, 

dl-n 


m 
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2 ?atG constant for dissociation of initiator in 
polymer phase, 1/cms 

rate ^constant for dissociation of initiator in 
aqueous phase, 1/cms 

kp rate constant for propagation, an ^ /mole, sec. 

^to ^ate constant for termination by combination, cm^/niole. sec . 

k^d rate constant for termination by disproportionation, 

crn^/mole sec 

k^^ 1 ‘ate constant for chain transfer to transfer agent, 
cm^/mole sec. 

^'fm 3:ate constant for diain transfer to monomer, cm^/mole sec, 

k^ rate constant for radical desorption, cms/sec 

L(n,i,'v,t) number of live chains of chain length n in all 
the particles having i radicals and possessing 
volurae between v and v+dv at time, t, moles/cm^ 
emulsion 

m number of micelles, mole s/cm ^ emulsion 

[M]p concentration of monomer in the particles, moles/cm^ 

[M]j^ concentration of monomer in the micelles, moles/cm^ 

[M]rji total concentration of monomer in the reactor, moles/cm^ 
emulsion 

n chain length of polymer, dimensionless 

r radius of particle, ems 

Rp rate of initiation, moles/ecc. 

rate of absorption of radicals in particles, and 
micelles, moles/ sec. cm^ 
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[S]^ total emulsifier concentration, moles/an^ emulsion 
^wc critical miscelle concentration, moles/cm^ water 
t time, sec. 

[Tr]p concentration of chain transfer agent inside the 
particle, mole s/cm ^ 

V volume of particle, cm^ ' 
volumo of micelle, cm^ 

v^ Volume of monomer droplets, cm^ 

Vp volume of reactor, 

volume of afiueous phase per cm^ emulsion, cm^/cm^ 

emulsion 

9 monomer volume fraction in the particles, dimensionless 

y monomer polymer interaction parameter, dimensionless 

Y surface tension, dynes/cms 

© mean residence time, sec. 

‘^ig Exoneckor delta function, 1 for i=J, 0 for 

6(V"-v ) Dirac delta function: 1 for vsv , 0 for vj^v , 



CHAPTER 4 

IgTIIOI) OP SOLUTION OP THE MODELLING- EQUATIONS 

In the previous chapter, the modelling equations for 
emuLsion polymerization in a CSTR were presented. These 
equations contain various material balance equations and 
the population balance equations. The material balance 
equations are ordinary differential or nonlinear algebraic 
equations and are relatively easy to handle. However, the 
various population balance equations involve different distri- 
bution functions like P(v) , f(i,v) , L(n,i,v) and G-(n,i,v) 
and are complex coupled partial differential -difference 
equations. Special care is required to decide the line of 
approach to solve these equations. In this chapter norma- 
lisation procedures have been introduced and for convenience 
the additional nomenclature for this chapter is listed at' 
the end of the chapter. 

The aim of solving these equations is to calculate the 
values of the aforesaid distribution functions. However, 
generation of these functions numerically for different values 
of attributes may take exorbitant time. To avoid this, the 
problem can be approached by making use of the concept of 
moments of a distribution^ 

Let f(x) represent a continuous distribution for the 
. Then the kth moment of f(x) is defined as 


attribute x 
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f 


k 



(4.1) 


0 

Similarly if H(i) represents a distribution defined 
only for discrete values of i, then the kth moment of the 
discrete distribution is given by 


i-i max 

H i^ H(i} (4.2) 

i-i min 


It can be seen from these definition that, with the 
knowledge of a distribution function, all of its moments can 
be calculated. Obviously, a continuous distribution fiinction, 
say f (x) , containing infinite information cannot be exactly 
specified by a finite set of moments. However, it is possible 
to represent a distribution function, in terms of finite 
number of moments to an accuracy sufficient for modelling 
purposes. One of the methods to achieve this is to represent 
the distribution function in a series expansion of orthogonal 
polynomials like laguerre polynomial [17,30]. The Laguerre 
polynomial iias the following orthogonal property, 

i ‘ 

(x)“ exp(x) 6nm (4.3) 

0 


where 6nm 


where 


si if n=m 

3 0 if 

3=0 



(4.4) 


(4.5) 
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Using the definition and orthogonal property of Laguerre 
polynomial, as stated above, a normalized distribution 
funccion G(y) can be represented in terms of summations 
involving Laguerre polynomial [40]. 
a ^ 

G(y) - Ln“(y) (4^6) 

I a+1 n=0 

Lhe coefficients Cn°^ are related to the moments of 
the normalized distribution function G(y) by folfLowing 
expression; 

C“ = ^ (-1) J J:h_^ I (4.7) 

j=0 - {jW U. Va+J+1 ^ 


a 


Equation 4.7 follows from the orthogonal properties of 
. Following Min [30] , the paraneter a is chosen so that 
and C 2 are zero, which results in rapid convergence. Lhe 
resulting expression for a is the following: 

“2 


a 


-1 + 




G, 




(4.8) 


The normalized function G(y), used here is related to 
the absolute distribution function f(x) as follows: 
f(x)dx = G(y) dy f ^ 

as 

f , 

and y » ^ — (4.9) 

£^(cc+X) 

The moments of the two distributions, f(x) and G(y)., are 
related to each other by the following expression 
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■•k 


r . '1 - 

ffTcSoT'^ ®: 


k-1 


( 4 . 10 ) 


It can be seen from the preceding lines that most of 
the informa/tion contained in a continuous distribution 
function can be obtained from a finite set of moments of that 
distribution. The number of moments required to specify 
a differential distribution depends on the complexity of its 
shape. It has been found that only 3 or 4 moments are often 
sufficient to represent a unimodal distribution satisfactorily. 

Even in practical cases, it is often impossible to 
measure complete molecular weight distribution or particle 
size distribution. Generally it is convenient to measure 
and use the leading moments of these distributions for quality 
control. This method involves the use of or poly- 

dispersity index to represent m,w,d, and average size to 
specify the particle size distribution. 


Eor the reason stated above, the modelling equations 
given in the previous chapter involving differential distri- 
bution functions will be transformed to a set of equations 
involving finite number of moments of those distributions. 

It should be noted that the method to represent a 
distribution in terms of its moments is useful only for 
continuous distributions. Obviously a discrete distribution 
containing a large number of closely spaced points can be 
approximated by a continuous distribution. Consequently, in 
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the present study the live peljoner chain distribution as 
well as the dead polymer chain length distributions which 
are actually defined for integer values of n only^ are 
approximated by a Continuous distribution in n and the 
aforesaid method is then made use of. 

Howeveri a distribution of polymer particles in i 
(the number of radicals contained in the particle) is defined 
only at a few points. It has been found from experiments 
•fehat the number of live radicals contained in a particle 
rarely exceeds 10. As a result, the use of moments to 
represent the discrete distribution in i is not justifiable. 

Hence this distribution will be generated numerically at the 
discrete points at which it actually exists (i.e. at i = 0,1,2.*). 
Since the number of points at which it is defined is small 
this technique involves the solution of a relatively small 
number of equations. 

To start with^the following assumptions will be 
made for mathematicai simplicity, 

1, A distribution in the volume of particles is 
generated* The fractional moments in volume (i.e, two third 
and one third moment) are used to calculate the total surface 
area and the average diameter of the particles. However, 
when 't* is connected with the distribution fmctions in a 
more complex way (like the term arising due to termination 
of radicals in particle) v^^e average value of the volume 
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of particles is used for mathematical simplicity. This 
approximation is likelj/' to be more exact when the particle 
size distribution is relatively narrow, 

2, The terms representing the chain transfer’ to 
transfer agent and initiation in particles will be maintained 
while developing the sets of equations involving moments 

of live and dead polymer chain length distribution. However, 
the values of [Tr]p and [I]p can be evaluated by addition 
of two material balance equations and thus these mechanisms 
can be easily accounted for. 

3. The value of monomer concentration in particles ( 9 ) 
can be calculated from equation 3.30. This equation shows 
that ( 9 ) is in general a function of particle size. However, 
from experimental studies it has been found that the concen- 
tration of monomer in a particle is almost constant and 
independent of its volume. The values of the concentrations 
of monomer in a polymer particle are reported in literature 
for different monomers. Hence,- these values can be directly 
used to solve remaining modelling equations. Most of the 
results from this study have been obtained by this method. 
However^ the results obtained without making this assumption 
are also given for comparison. 

The modelling equations given in the earlier chapters 
are not diraensionless. Hence the numerical values of various 
variables depend upon the system of units to be used. These 
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equations should firstly be transformed into sets of 
nondimensional modelling equations* This approach facilitates 
the use of these equations, since the variables are not 
confined to any particular system of units* Similarly it 
facilitates the numerical solution of these equations, since 
the range of numerical values of various dimensionless 
variables is relatively narrow. The dimensionless population 
balance equations should then be transformed into the sets 
of equations involving moments of the dimensionless distri- 
butions. 


This procedure has been explained in the following 
lines, for the case of live polymer chain length distribution, 
Consider the population balance equation ( 3(il0) 
involving live polymer chain length distribution; 




■] 


+ - 3 - (r .Lfn i V t) )~ R T £— 

[L(n,i~l,v,t) - L(n,i,v,t) + f (i-l,v,t) ‘6^^] 


+%[I]p [lCn,i-2,v,t) - l(n,i,v,t) + 26^p*f (i-2,v, t) ] 


0 Jr 


[ (i) L(n,i+l,v,t) - (i) L(n,i,v,t)] + kp[M]p 


[L(n-l,i,v,t) - L(n,i,v,t)] + kp^[Tr]p [i f(i,v,t) 6^^ 


.L(n,i,v,t)] + kpjMjp [i f(i,v,t) 6^^ - L(n,i,v, t) ] 

iiiiiXil I,(n,i+2,v,t) - L ( n , i , V, t) ] 


~Ttv 


I [I{n,i,v,t) t -L(n,i,v,t)] 
^ feed 


A. 

+ R . T V —a” \ 6 t 6 . n d (v— V ) 

3a(Ap+J^) nl il '■ m' 

( 3 . 10 ) 
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At steady state the derivatives with respect to 

time can be equated to zeroi 

Similarly all the particles are assumed to have volume 

^avg* Hence the derivatives with respect to v also reduce 

to zero, and the average Volume of the particles V can 

avg 

then be substituted in place of volume, v, in these equations. 
Using the definitions of various dimensionless variables and 


dimensionless constants given at the end of this chapter 


and after 
0 sa 


further simplifications this equation reduces to 

ii(n,i-l) -ECn.i) + f(i-l) 

^ 2 / 3 “ 

+C! 2 [I]p*v[i(n,i- 2 ) - E(n,i) + 26^^ f(i-2)] 

(v7^''^^[i L(n,i+l) - i l(n,i)] 

+ C 3 _Q^( 9 )[L(n~l,i) - l{n,i)] 


+ 0^3 Cfr]p T(i) - l(n,i)] 

Q. - 

^102^^^ I(n,i)] + (i+1) 

L(n,i+2) ~ (i)(i-l) L(n,i)] + ^ [L(n,i) t -L(n,i)] 

^ feed 

( 4 . 11 ) 


Ihis dimensionless population balance equation is then 
converted into a set of equations containing moments of 
the distribution in n. To do this the equation is multiplied 
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by n and. is smimed over all possible values of n* After 
furthei' simplif ieatioiis the following • equation is obtained; 

0 .. (i) 4-f(i-l)] 

(F2/3 + la) ^ 

[bjj;(i-2) - Ljj.(i) + 2 f(i~2)] 

+ C^(v)^/^ [i L^(i+1) - (i) Lj^(i)] 




m=l m k-m 


+Cg3[lr]p[i ?(i) 


^Ijj.(i)] + ^102^^^ t (i+ 1 ) (i)EJi+ 2 ) 


.(i){i-l) Lj^(i)] + i [Lj^(i) - L^(i)] 


feed 


C7CI] 


i' 2/3+ ” 


(4*12) 


other modelling equations can be rendered dimensionless 
and the population balance equations can be transformed to a 
set of equations involving moments of dimensionless distri- 
butions, by using the aforesaid approach. The resultant set 
of modelling equations is given on the following page. 
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Pa r t icle Size Distribution ; 

0 = I .n I, , . 0, iS^ 


0 ^ feed ^ 

R adical Fumber Pistribution t 


l!:=0 to n 


(4.13) 


0 = Or- 


t f(i-l) ~ f(i)] + C^(vr^/^[(i+l)f(i+l) 


^^2/3”^ 21) 
G, 


« i f(i)] + -^ [ (i+2) (i+l) f(i+2) -(i) (i-l) f(i)] 

( v) 


^ 1 [|(i) 


m' 


(P2/3+ S) 

+ 02{v)[f]j [f(i-2)-f(i)] 


■Q feed 

i = 0 to m 


-f(i)] 


(4.14) 


Live P oly m e r . Chain L e ngth. Distribution : 


0 5= C . 


[i]^w 


273 


{I'2/3 + 5) 


— [i^i-D - Lj.(i) + f(i-l)] 


+02CLp(^) + 2 f(i-2)] 

+ Sjj(i+1) - CD \(D] + Cj^ol 


ms»l 


m k-m 


C 


4- "5(4) - Sjj.(i)3 + j^^[Ci+l) (i) Lj^(i+2) 


(v) 




0 


feed 


+ ^ [4] - s 6^-). iCiaO to n; i=0 to m (4.12) 
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Dead P oly mer Chain 


07[i],.,(v)2/3._ 


'-''tL -^J V V ; ’ - 


Gij(i+l)-(i)G)^(l)] + C2[I]p(Y) 

, ' 1= c - _ 

[Gr^^(i-2) - Gjj.(i)] + “^ [(i+2)(i+l) ^(i+2) 

, 0 I ^ 

-•(i)(i-l) &jj(i)] + [f(i+2) XI — — g (i+2)g, (i+2)] 

(^) p=ofe:^^ ^ 

^ ^102(^)CSk(i)] [(i+1) Ek(i+2)] +Cg3[Tr]p Lj^(i)] 


^ I J - ^ 


feed 


k=0 to n 
i=0 to m 


(4.15) 


Emulsifier Balance: 


0-kis^] 

9 ^ feed ^ 


(4.16) 


Initiator Balance ; 


ii]„ + ^ [[in \ - [I] j 


feed 


(4.17) 


Monomer Balance? 


[[% n°53'?> ^o- lavs 


(4.18) 


Micello Balance; 


0 =* -m + C^Q [S]^) "" %i *“ ^ 2/3 ” ^52 ^d 


(4.19) 
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Monomer Droplets Balance 


° = - 055-^i'-(?) 

E xpression fo r (m) ; 

0 = (1-9) + in(^) + »# (1-9)2 + 0 


OjjgCi) 


( 4 . 20 ) 


'65 


F 


( 4 . 21 ) 


1/3 


It shou3.d be noted that the expressions 4.12, 4.14, 4.15 
are functions of k and i, Simila,rly the expression 4.13 is 
a function of k* Thus, these e^cpressions indeed represent 
sets of equations depending on different values of k and i. 

The particle size distribution is continuous and the 
distribution in n can be assumed to be similar to a continuous 
distribution* Hence the different moments are used in the 
aforesaid sets of equations to represent them. However, it 
can be seen from these expressions that different fractional 
moments with respect to volume (i.e, and ^ 2 ./'^ need to 

be evaluated separately. This can be done by expressing the 
differential distribution in volume in terms of first 5 integer 
moments as described previously in equations 4*4 to 4.10. Any 
fractional moment can then be calculated from the differential 
distribution. The resulting equation expressing the fractional 
moment in terms of first 5 integer moments is given by 


^ la 4- t -f 1 


m 


n 


^ n=3 


X<f ^(-1) 


k 


n 


k=0 


JL?L. t' g Hh 1 I g-ft-fk-fl 


^n~k+l \a+k+l ' \ k+l 


( 4 . 22 ) 
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TliG dis tr ib\itiorL in i, however, cnnnot be approximated 
to behave like a continuous distribution, as explained 
previously. Hence, the equations inolving i (i.e, 4.12, 4,14 
and 4.15) are solved for different discrete values of i. 

These equations, however, involve advanced terms calculated 
at (i+l) and (i+2). Thus the set of n equations will involve 
n +2 variables. However , as stated earlier, the number of 
radicals in a iDarticlo is generally limited to a small value. 

In other words, after a sufficiently high value of i the 
toms eviiluated at (i+l) and (i+2) tend to become zero. Thus 
the advanced terms requiring evaluation of (i+l) and (i+ 2 ) can be 
approxiraat ed to zero and the set of n equations can be solved 
for remaining n variables. 

Me thod of Solution : 

The modelling equations given earlier can be divided 
into 3 groups. The first group of equations contains various 
material balance o<jua.tions, the equation for ( 9 ), the equations 
involving moments of particle size distribution and the 
equations involving radical nirnber distribution. Equations 
4 . 13 , 4*14, 4.18, 4 . 19 , 4.20 and 4.21 are interdependent 
and should be solved simultaneously. Equations 4.16 and 4.17 
can bo solved separately. The solution of this first group 
of equations gives the values of the various dimensionless 
concontration teims, radical number distribution and moments 
of particle size distribution. Using these values the second 
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group oi equations involving moments of live polymer chain 
length distribution can be solved to obtain for 

different values of i and k* These answers can further be 
utilized to solve the third group of equations involving 
dead polymer chain length distribution to obtain (i) for 
different values of i and k. 


The basic him of solving these modelling equations 
was to evaluate the properties of the output stream from the 
reactor. These properties include conversion, rate of 
polymerisiation in the reactor, number of particles per cm^ 
of emulsion, particle size distribution and the molecular 
weight of the resul-tant polymer. The conversion, rate and number 
of particles per cm^ emulsion can be calculated by knowing 
the values of and 


The expressions^ to obtain these are given below; 
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( 4 . 23 ) 

( 4 ; 24 ) 


N 
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( 4 . 25 ) 


The particle size distribution of the latex can be 
generated from its moments Pjj. by using equation 4^22. 
average volume and average diameter of these latex particles 
can also b© evaluated from these moments. 
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(4i26) 

(4.27) 

(4.28) 


111© number average as well as weight average molecular 
weight and the polydlspersity index can be calculated from 
moments of dead polymer chain length distribution* These 
expressions are given below: 


where 
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(4.29) 

(4.30) 

(4.31) 

(4.32) 


The entire molecular weight distribution of the 
resultant polymer can also be generated from these moments 
by using the procedure outlined earlier in equations 4.4 to 
4 . 10 and 4 • 22 « 

JDescrintion of Computer Program t 

The equations in group 1 are subdivided in two sub- 
groups, first containing equations 4.16 to 4.21 and 4.13 
while the second comprises of equation 4.14. For the solution 
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of the firsc subgroup the feed conditions and initial guesses 

for vario-oles as well as initial guess for i is needed* 

avg 

Ihis group ox ec^uations can then be solved to give va-lues 
of diiferenc dimensionless concentrations and the moments 
of par bide size distribution. These values are then fed 
bo the second subgroup which now becomes a set of linear 
equations. The starting value of m=5 is used. The values 
of J(m+1) and f(m+2) are assumed to bo zero and the equations 
are solved to get values of f (i) ; i=l,m. An increment of 
two is then made in m and (m-i-2) equations are again solved, 
assuraing f(m-i-3) and f(m+4) to be equal to zero. The values 
of f(i), i=*l, ni+2 aro then obtained. The values of f(m+l) 
and f(ia-f2) are checked to be less than a small constant z 
(lO” was used) , If these values are not less than z an increment 
of two is again made in m and the set of equations is again 
solved. The process is repeated until the last two variable 
values of a set are less than z, thus checking the assumption 
made in the solution of previous set. The value of i„„_ is 
thon calculated , This improved value of i„,,„ is then 

dVg 

compared with the previous value of and is used together 

with the values of different variables as the improved guesses 
to solve the equations in the first subgroup. The process 
is repeated till the difference between two values of i^^^ 
is less than a particular convergence criterion. 

The equations in the first subgroup are highly nonlinear 
in nature* The method suggested by Brown has been employed 
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to solve these equatiojis. The details of this method were 
taken from reference [31]. It has "been found that it is 
not possible to give good initial guesses for these equations, 
withouG which che method does not tend to converge to the 
solution. Because of this reason an objective function 
Consisting of stmis of absolute values of RHS of the aforesaid 
equations is formed. This objective function is then minimized 
with respect to the aforesaid variables using the multivariable 
constrained optimization technique of Box, The details of 
this technique are given in references [32] and [35]. Relatively 
wors'Q initial guesses can be used for this optimization 
procedure* The values of different variables from the opti- 
mization procedure are then used as the initial guesses for 
solving the equations of first group. The solution of linear 
equations is done by standard matrix inversion procedure. 

The subprogrammes needed for solution of this group are ROHLIIT, 
BACK, MAXI, MIMl, FURZ, CORSX, CHECH, CEHTR, FUHC, CONST, 

MATH, SETl, LINEAR, GOLDEN, COM and OBJECT. 

The answer from first group of equations are then used 
for the solution of eqtiations in second group. While solving 
the radicca nmber distribution (i,e. second subgroup of group 
1) , the assumption made for solving a set of equations was 
immediately checked by solving an increased number of equations. 
Hence when the convergence was ultimately achieved, the 
assumptions made in the solution of last but one set of 



51 


equations (containing say^ *p* equations) was verified by the 
soluGion ox last set of equations. Hence these *p* equations 
should, in principle, be solved for live and dead polymer 
chain length distribution.* However, the solution of live 
polymer chain length distribution fox* f(o) type of pai’ticles 
is meaningless* Hence *p— 1’ equations are solved for live 
polymer chain length distribution. However, for dead polymer 
chain length distribution 'p* equations need to be solved* 

Ihe equations are linear and are solved by the same method 
employed earlier. The subprograiames required for this part 
are SET2, MATIH said CONI'R. 

The answer from 2nd group of equations is then used 
to solve the equations in 3rd group. The method of solution 
of these equations is similar to that employed for the 2nd 
group* The answers from this pajrt give the values of 
for different values of i and k. The subprogrammes required 
for this portion are SIi)T3, MATIN and COHPR. 

After the solution of modelling equations various ■ ' 
properties of output stream are calculated by using equations 
4*23'*4#3'2* Cumulative particle size distribution, live 
polymer chain length distribution, molecular weight distri- 
bution are then generated. The subprogrammes required for 
this portion include DISTR, DIPH2, COHB’R and QGIO. ^ 

The diagram showing the scheme to solve the modelling 
equations is given on the next page, 

CENTRAL LIBRAif 
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MENGLATURE 

(1) = nondimeiisionalizatioii constant for distribution 

functions, micelle concentration and droplet 
concentration, moles/lit. 

(2) [Iq]^= nondiniensionalization constant for various 

concentration terms, moles/lit. 


Definitiona of Dimensionl ess Variable s : 

1 , 
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CIIAPTER 5 


RESULTS Am DISCUSSION 

The modelling equations developed in the earlier 
chapters have been solved by using the approach outlined 
in the previous chapter. The model has been applied foi* 
the reaction conditions reported by DeG-raff and Poehlein [21], 
Their work is one of the most readily available and 
comprehensive work on emulsion pol 3 rmerizat ion in a CSTR, 

The experimental set-up used by these authors consisted of 
a jacketed stirred tank reactor fitted with a temperature 
controller. Nitrogen atmosphere was maintained inside the 
reactor. The reactants were fed in two different streams j 
first st real'll comprising of monomer (styrene) and the other 
comprising of aqueous solution of initiator (Ammonium per- 
sulfate) and sujrf actant (sodium laLiryl siilfate) , The effluent 
streaia was collected in a stirred vessel containing hydro- 
quinone inhibitor and was then used for experimental measure- 
ments. They have also developed a theoretical model, using 
residence time distribution (RTD) approach, to obtain particle 
size distribution, rate of polymerization, number of particlee 
and number average molecular weight of the polymer product. 
They have carried out experiments at various reaction 
conditions and have tried to compare the results with the 
values obtained frora their model. Computations using the 
model developed here will be compared with their results. 
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Vario^IS constants used for the solution of the modelling 
equations are listed on the next page* 

As has been stated earlier the chain transfer to transfer 
agent and the initiation in particles have been neglected 
in the present work, These mechanisms can, in general, be 
neglected if the monomer feed is free from transfer agent 


or trace impurities and if it is insoluble in water. Hence 
the values of k^^ and k^p are assuiiied to be zero, k^^ is 
also assumed to be zero, since for styrene the contribution 
due to termination by disproportionation is usually negligible 
compared to that with termination by combination [37]. It has 
found from earlier studies on batch systems that desorption 
from particles play a relatively small part in emulsion 
polymerization of styrene [15], Hence the value of k^ is 
assumed to be zei'o. It has been reported that the micelles 
swell and approxiBiately double their size when saturated 

r* 

With monomer [36] , This implies that approximately 75 per 
cent of the volume of a saturated micelle is filled with 
monomer which results in the value of [M]^^ as 0,0063, 

The reaction was carried out by DeG-raff and Poehlein [21] 
at several temperatures and surfactant concentrations. Only 
the runs at 70°C and with a surfactant concentration of 2.79 gms/ 
100 gms water were selected for simulation. The monomer 
to water ratio in the reactor feed was 0.54. From this 
ratio, is calculated to be 0.6494. The value of ( 9 ) was 

' W 
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Constant 

Unit s 

Value Ref . Uo , 

1. 


cm^/mole, sec 

2.650x10^ 

[21] 

2. 


CEis/ sec 

0.0 

assumed 

3. 

% 

1/ sec 

2.55x10"^ 

[21] 

4. 

^iP 

l/sec 

0.0 

assumed 

5. 


cm ^ /mole sec 

0.0 

assumed 

6 , 

^tc 

cm^/mole sec 

7.2x10^ 

[21] 

7. 

^ft 

cm^/mole sec 

0.0 

assumed 

8. 

Ir 

fm 

cm^/mole sec 

21.2 

[54] 

9. 

dp 

gms/cm^ 

1*0S^K 

[21] 

10. 

n 

gw5j cm^ 

0^8780 

[21] 

1 1 • 

V 

m 

cm^ 

65.4498x10“^ 

[15] 

12. 

"^mol 

cm^/gm.mole 

118.45 

[21] 

15. 

■^d 

cm^ 

4.1888x10“^^ 

assumed 

14. 


gm/mole 

10'^ 

- 

15. 

[Hln. 

noles/cm^ 

0,0063 

[3S] 

16. 

OE 

cm^/mole 

5.672 xlO"^^ 

[21] 

17. 

Y 

dynes/cms 

7.50 

[10] 

H 

00 

• 

y 

dimensionless 0,256 

[10] 

19. 

^■wc 

— fi 

moles/cm^ \mter 8,9x10"’ 

[35] 
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reported by DeG-raff and Poehlein [21] to be 0.615 and has 
been used for nost of the siuiulations. 

Calculations were made using the aforesaid constants 
at various reaction conditions and the results were compared 
with the theoretical predictions and experimental results 
obtained by DeGraff and Poehlein [21]. 

Particle Size Distribution ; 

The cumulative particle size distributions calculated 
at different conditions have been presented in Figures 2 and 
5. It can be seen from these plots that the results are 
practically identical to the theoretical results obtained by 
De Graff and Poehlein. 

As can .be seen from the figures, the breadth of the 
particle size distribution increases with increasing residence 
time. Wii^h the increase in the mean residence time, the 
distribution of times spent in the reactor broadens. Further 
the rate of voluiae growth (r^) will also increase due to the 
increase in the value of Both these effects cause a 

broader particle size distribution. 

In the RTD approach it is assumed that a fraction of 
the micelles fed to the reactor do not get activated and 
the rest of the micelles are converted into particles at the 
same instant. Hence all the particles should have the same 
birth time. The RTD approach followed by De Graff and Poehlein 
then requires the knowledge of only two mechanisms; (1) kinetics 
of volume groirbh of particles and (2) RTD of the reactor. 






mm: 
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In the present work this assumption of all particles being of the 
same age was not made to calculatei particle size distribution. 
However, the distributions obtained in this work coincide 
with those obtained by DeGraff and Poehlein. Such coincidence 
could be explained if the kinetics of volume growth of 
particles is either first order ar zeroth order. However, 

DeGraff and Poehlein used Stockmajrer ' s[2] expression for 
growth rate which is neither a first nor a zeroth order process. 
At the moment, an explanation fox* this coincidence cannot be 
offered. 

Rate of Polymerization and humbe;^ of Particles ; 

( a) Effect of Me a n Re s idence Time i 

The effect of mean residence tame on the rate of 
polymerization and number of particles has been presented 
in Figure 4 , With increase in the residence time, Ute average 
size of particles in the reactor increases and requires more 
amount of soap to stabilize thou, Phis soap is supplied by 
the micelles. Hence the number of micelles in the reactor 
decreases with increase in the mean residence time. The 
decrease in number of micelles leads to a decrease in the 
number of paarbicles as the residence time is increased as 
shown in Figure 4. 

DeGraff and Poehlein have assumed that soap is present 
only as micelles or as stabilizer for polymer particles. 

They have not taken into consideration the emulsifier present 



crization, moles/ lit. hr 



T=70°C 

(S)x = 2-79gms/100gms water 
(I)w=0-08gms/100gms water 


DeGraff 
Present work 


9. mean residence time (hrs) 


Effect of mean residence time on number 
of particles and rate of polymerization. 
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on the sui’face of monomor droplets and that present in 
aqueous phase upto critical micelle concentration, to calculate 
the free emulsifier available for generation of the new 
particles. Hence their values of the free emulsifier available 
for particle formation are high. Since they assimaed the 
number of particles to be proportional to the free emulsifier 
available, their values for number of particles should in 
general be on the hi^er side, Fiirther, DeGraff and Poehloin 
have assumed that all the particles are generated at the same 
time andhavQ the same birth time. However, in reality, the 
micelles may reside in the reactor for some time before 
getting converted into particles. Thus, actually, some of the 
micelles having low residence times, will not be converted 
into the particles. Hence the lamber of particles calculated 
by De Graff and Poehlein will be still hi^er than the actual 
number of particles. The assmptions made by DeGraff and 
Poehleini stated above, have not been made in the present 
work. Therefore, the present results should give better 
approximation for number of particles, particularly at low 
residence times. It can be seen from Figure 4 that at low 
residence times a considerably better fit to the experimental 
data has been obtained using the present theory, while the 
results obtained by BeGraff and Poehlein are far from the 
experimental points. The deviation from theoretical predictions 
at high residence times is however inexplicable. 
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The rate 01 polymerization will be directly proportional 

to bhe number of particles times i____ as can be seen from 

avg 

3.28 and 4.24, Hence an explanation similar to that given 
in the previous paragraph can be given about effect of mean 
residence tirae on the rate of poljnaei'ization also. The 
number of particles decreases while the average number of 
radicals per particle increases with increase in mean 
residence time. The variation of rate which involves the 
product of these two terms in given in Figure 4. The almost 
identical variation of the rate of pol 3 nnerization and number 
of particles implies that the average number of radicals in 
a particle does not vary much with increase in the mean 
residence time, 

("b) Effect of Initiator Concentration ; 

Figure 5 shows the effect of variation of initiator 
concentration on and M„. When the usual values of mean 

residence time and concentrations of initiator and soap are 
used the rate of initia,tion of radicals is considerably hi^ 
Compared to the free emulsifier available. Hence under these 
reaction conditions the number of particles in the reactor 
is considerably larger than the number of micelles. When hi^er 
concentrations of initiator are used hi^er rates of generation 
of radicals are obtained. However, most of these radicals 
will be captured by the particles since the total area of 
the particles is considerably larger than the total area of 
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the previous paragraph the value of number of particles 
is almost consta,ri.t "with respect to initiator concentration. 
Hence the rate of polymerization, ■which is proportional to 
the product of number of particles and the average number of 
radicals Toer particle, should not vary significantly with 
changes in initiator concentration. This behaviour can be 
seen from Figure 5. The present theory follows the experi- 
mental data reasonably well. 

Molecular Wei,shts ; 

(l) Effect of Mean Residence Time ; 

The effect of mean residence time on M^ and M.^ is shown 
in Figure 6. With increase in mean residence time there is an 
increase in the particle size and h.on:cG the total rate of 
entry of radicals into particles increases. This decreases 
the tiJiie of growth of radicals inside the particles. Moreover, 
the average nuraber of radicals present in a particle increases 
with size leading to the faster termination and lower times 
of chain gro-wth. Due to these reasons the live chains 
present inside the particles grow for shorter times resulting 
in lower molecular wei^ts. This effect can be seen from 
Figure 6 * 

D eG-raff and Poehlein [21^ did not consider the effect 
due to chain transfer to monomer for calculating M^, However, 
this phenomena in general will always be operative and will 
tend to bring down the molecular weights of the polymer. This 



9 , noGan residencG time, hours 


g. 6 -Effect of mean residencG time on 
molecular weights. 
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can be seen from Pigure 6 where the variation of molecular 
weights have been shown with and without considering chain 
transfer to monomer. The values obtained by the present 
theory are found to be better approximations to the experi- 
mental results than those obtained by DeG-raff and Poehlein’ s 
theory. However, the values still do not agree with the 
experimental results satisfactorily at low residence times. 
Better fit could be obtained by taking into consideration 
the chain transfer to impurities and the termination by 
disproportionation. These processes will tend to bring down 
the molecular wei^t further. 

DeGraff and Poehlein from the ETD approach could not 
obtain the wei^t average molecular weight of the product. 
However, the present theory allows the calculation of from 
moments of dead polymer- chain length distribution. These 
values are plotted in Figure 6 with and without considering 
chain transfer to monomer. It can be seen that these values 
compare fairly well with the experimental results. 

(2) Eff ect of Initiator Concentration i 

The effect of initiator concentration on molec-ular 
weights has been presented in Figure 7. As the initiator 
concentration increases the rate of radical absorption in 
the particles will increase. With the increase in the 
frequency of radical absorption in a particle, the live 
chains inside the particles are terminated faster, resulting 
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in a small molecular weight polsnaer as can be seen from 
Figure 7 . 

The results obtained by DeGraff and Poehlein and by 
the present theory with and without considering chain transfer 
to monomer are compared with the experimental data. It 
can be seen from the experimental results that the results 
obtained by the present theory considering chain transfer to 
monomer give better approximations to the data than DeGraff* s 
results and represent the experimental data well. 

The results discussed so far have been obtained by using 
the value of (9) to be equal to 0,615 as used by DeGraff 
and Poehlein. However, this value can be calculated by using 
equation 3,31, This method was employed to obtained the 
values of number of panbicles, rate of polymerization, average 
size and the molecular weight of the resultant polymer latex 
in simulating one run. The results are compared below with 
the results obtained by the earlier method at the same 
reaction conditions: 

1 , Temperature = 70 °C 

2, [S]^ = 2.79 gms/100 gas water 

3 , [I]^ = 0.35 gms/lOO gas water 

4, G =15 minutes 
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First Approach Second Approach 


1.9 

0.615 

0.801 (calculated 
from 4. 21 ) 

2. Average diameter 
(number basis i) 

729 

986 

3. Average diameter 
(volume basis i) 

817 

1104 

4. Bumber of i^articles 
per cm^ of emulsion 

1 . 09 x 10 ^^ 

6 . 12 x 10 ^'^ 

5, Rate of polymeri- 
zation moles/hr 

4.7 

3.47 

6 . 

5.46x10^ 

4.006x10^ 

• 7 . \ 

1.32x10^ 

9.392x10^ 

8 . Polydispersity 
index 

2.43 

2.344 

It can be seen from the table given 

above that the 

value of 9 obtained by 

the latex method ( 

i.e, 9 = 0.801) is 


different from the one used in earlier method (i.e, 9 = 0 . 615 ) « 
Because of the hi^aer value of 9 obtained in the latter 
method the rate of volume growth would be higher and hence 
the average particle size obtained is higher. Because of 
the higher size of particles, they require more soap for 
stabilization and henceless amount of soap is present as 
micelles, This results in lower value for number of particles 
by the latter method. Similarly, because of the larger size 
of the particles the frequency of radical entry in a particle 
is increased, leading to shorter growth times for the radicals 
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present in the part iclei..' This results in the lower values of 
molecular weights as reported in the table. 

The aforesaid comparison is made without considering 
chain transfer to monomer, 'The comparison was also made by 
taking chain transfer into consideration. The results are similar 
to the case where chain transfer is not considered. The 
molecular weights obtained by both the approaches are, however, 
smaller than the first case. 


Pirsi: 

: Approach 

' Second Approach 

1. 9 

0.615 

0.819 (Calculated 

2, Average diameter 
(number basis A) 

714 

from 4.21) 

986 

3 , Average diameterQ 
(volume basis A) 

799 

1104 

4. Humber of particles 
per cm^ emulsion 

1.14x10^^ 

6.12x10^'^ 

5, Rate of polymerization 
moles/hr 

4.78 

3.47 

6. 

3.925x10^ 

3.092x10^ 

7. 

9.215x10^ 

7.074x10^ 

8. Polydispersity index 

2.347 

2.287 


Apart from the results stated earlier, the present 
theory has the merit of ability to predict the entire molec\ilar 
weight distribution of the polymer product. The molecular 
weight distribution on weight basis for a particular set of 
the reaction conditions is plotted in Pigure 8, 



9gms/100gnns water 
35gms/100gms water 
inutes 




CHAPTER 6 


COHCLUSIONS 

1, The model presented in this work predicts the 
experimental data on nnmher of particles, rate of polymeri- 
zation and molecular weights better than the earlier model 
proposed by De Graff and Poehlein [21]. 

2, The results for particle size distribution are 
identical to those obtained by the RTE approach [21], 

3. The present model allows the complete moleculs-r 
weight distribution of the polymer product to evaluated. 

4. The model can be made more effective by including 
the phenomena like chain transfer to impurities and termi- 
nation by disproportionation. 
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THIS PROGRAM IS OSE0 TO SIMUIiATe A CONTINUOUS STIRRED TANK 

REACTOR FOR EMUDSION PODTMERIZATION.THE INPUT DATA FOR THE 
PROGRAM CONSISTS OF VARIOUS CONSTANTS, FEED CONDITIONS AND 
various reactor parameters, this PROGRAM IS ABDE TO CADCUltATB 
THE CONVERSION, RATE OF POLYMERIZATION, AVERAGE PARTICLE SIZE, 
NUMBER OF PARTICLES, PARTICLE SIZE DISTRIBUTION AND MOLECULAR 
WT.DISTIBUTION OF THE RESULTANT POLYMER LATEX, 

INITIALLY THE BOX'S OPTIMIZATION PROCEDURE IS USED TO 
CALCULATE INITIAL GUESSES FOR THE SOLUTION OF FIRST SUBGROUP,THE 
EQUATIONS or FIRST SUBGROUP Cl, E, MATERIAL BALANCE EQUATIONS, 
Particle size distribution and radical number distrlbutionlare 
THEN SOLVED, THE EQUATIONS OF SECOND GROUP ARE THEN SOLVED TO 
OBTAIN MOMENTS OF LIVE POLYMER CHAIN LENGTH DISTR, THESE RESULTS 
ARE THEN USED TO SOLVE THIRD GROUP TO OBTAIN MOMENTS OF DEAD 
POLYMER CHAIN LENGTH DISTR, VARIOUS OUTPUT VARIABLES ARE THEN 
CALCULATED AND THE REQUIRED DISTRIBUTIONS ARE GENERATED, 

THE DETAILS OF THE METHOD OF SOLUTION OF THE MODELLING 
EQUATIONS IS GIVEN IN THESIS tPAGE ], - 
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NOMENCLATURE, 

VARIABLES mmmmmmmmm 

XSTARTC 1 , 1 lii,Mf<ni«»,*f<r«,«i«iMi><««ii<p>«f«OTii«ai»iiKflp»«TOTAL SOAP CONCENTRATION 
XSTARTC 1 ,2) <pMni)i»ii«n„,<i*>»fn<<*iiii,p(K<iiti«r«K<ni<*«a»INITI ATOR CONCENTRATION 
XSTART C 1 f 3 )n«(Mwi,,iiir<va«ii,<»«,<p««t«,a,iir<»i*«*<pTOTAL MONOMER CONCENTRATION 
XSTARTC1,43 •i«<mi«i,»,*a,P!)mMtar<l,ara,M>«M>p<n>«*«i4*>i*MICELLE CONCENTRATION 
X5TARTC1,5)™-,,,„„.«*,,,««„««*».M0N0MER DROPLETS CONCENTRATION 
XSTARTC 1,63^,„„™«.,,,„.„«*.,.M0N0MER VOLUME FRACTION IN 
XSTARTC1,7),-,.*.™»*™,,.,»,**ZER0TH MOMENT W,R,T, VOLUME 
XSTARTCl,8),™„„,^.,*„«*,,*„„„riRST MOMENT W,R,T, VOLUME 
XST ART C 1 , 9 ) <„,,a„rwm,«<iKmmiwi««w<«»iww«i«ta*SECOND MOMENT W,R,T,' VCILUME 
X5TARTCI,10).™™,,„„.,^*,„,*,THIRD MOMENT W,R,T, VOLUME 


XSTARTClrll) 
XRADCI) 




FOURTH MOMENT W.R.f, VOLUME 


XLIVECK,!) 
XDEADCK,!) 
XFDCI) 


.RADICAL NUMBER DISTRIBUTION ^NUMBER OF 
PARTICLES HAVING I RADICALS, 

.MOMENTS OF LIVE CHAIN LENGTH DISTRIBUTION, 

IK TH MOMENT f 1 RADICALS IN PARTICLE,) 
.MOMENTS OF DEAD CHAIN LENGTH DISTRIBUTION, 

CK TH MOMENT f I RADICALS IN PARTICLE,) 

.FEED CONDITIONS FOR VARIABLES OF FIRST 5UB» 
GROUP, 

FFDCIJ 

rDLlV(K,I).,.,,,,*«.,.„rEED CONDITIONS OF XLIVECK,!) • 
FDDEDCK,I)«..,«,*,.«*,.FEED CONDITIONS Of XDEADiK,!) 

CONSTANT 

.PROPAflON RATE CONSTANT, 

.RATE CONSTANT FOR DISSOCIATION OF INITIATOR 
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IN AQUEOUS PHASE, 

CONSTANT FOR DISSOCIATION OF INITIATOR 

IN poltmer phase, 

AKZERO.,..,„»«»»,*,*DESORPTION RATE CONSTANT, 

AKTC.*w«».,*,,»,„^,,RATE CONSTANT FOR TERMINATION BT COMBINATIO 
AKTO«»»,^«...,„,„,,*.,.«RATE CONSTANT FOR TERMINATION 87 DISPROPOR- 
TIONATION* 

AKFT..«„**,«„«.«™.,RATE CONSTANT FOR CHAIN TRANSFER TO T 

TRANSFER AGENT, 

AKFM«**,„.««..,**»„^,,RATE CONSTANT FOR CHAIN TRANSFER TO MONOMER 

0^ POllYMERf 

DENSITT OF monomer, 

VMOItwmv^pwaittOTwwwwvtiivwaimpMOltAR VOUUMB OF MONOMER, 

SAI— •„««„„,,.«.*..„M0N0MER.P0U¥MER INTERACTION PARAMETER, 
XETA.«„.*»»,,„,.™«^.SURFACE TENSION , 
AMW,.„.,...,.«„«,.,,M01,ECUUAR WT. OF MONOMER, 
THETA,^«,„,„.,^,.«,^MEAN RESIDENCE TIME, 
TEMPw*i»!«,»»»»w>«w«*i(,,,ipr,,^,«»»i»««TEMPER ATU RE * 

VW,»»„*,«.»,„„»,.^».WATER!EMUI*SI0N ratio, 

AID C 10 7 N 1 SEC PAGE 53 
AMZERQj,„»i!p)„,iMi,<«<i,)*ipi,<in*v<mp C M 7 0 iSEE PAGE 53 
VDw«i,i»»i,*»i»*n,«ii««ii»ii,<w«,«ni*,,<»,ti*VOIjUME OF DROPUETS* 
VMi,nii,ip(i,*i»i(,»mm(iii,»mnp«pm,«ni,i«,VOI*UME OF MICEUIiES* 
CSii«„,,iHifi,«ri*>i«,HMi,Miiii,«imi,*,ii«f*,mpARCA OCCUPIED BY SOAP* 
AMM.».,.,,„».„..,^,N0N0MER CONCENTRATION IN MICELUES, 
FEE**.«,«...,.«,„«,,»ErFlCIENCy OP DECOMPOSITION OF THE 

INITIATOR, 

SNC™™«««™».,,«.CRITICAI4 MICEUUE CONCENTRATION, 


IPRIN««,«,«..«,.«,*.«PRINTING OPTION I IPRIN«¥l — INTERMIDIATE 

ITERATIONS ARE PRINTED | IPRIN:pO— • 0ND7 
DATA AND ANSWER FROM SIMULATION ARE PRINTE 
IPT.™.,,.™— ^.....PRINTING OPTION I lPT«sO— CONCISE OUTPUT 

IPfsi«*«DETAILED OUTPUT E ENTIRE DISTR, 

ARE CALCULATED/PRINTED? 

MX«„.™«*^™*««.*.«INPUT UNIT NUMBER TO BE TYPED BY THE USER 

NX. ...............RESULT UNIT NUMBER TO BE TYPED BY THE USER 

NX»l,20,2l,22,23,24—— RESULT ON DISK 
NX«5—— RESULT AT THE TERMINAL 


C 

DOUBLE PRECISION XSTART( 14, 167 ,R ( 14, 1 1 ) , FC 14) ,GC 167 ,HU67 
DOUBLE PRECISION XCCH7 ,X(2Q7 ,XTEMPC67 ,BC25,17 ,KRADC0l207 
DOUBLE PRECISION T(25,257 

DOUBLE PRECISION FFDC0S207 ,FDLIVCOl5,207 ,FDDED(08 S,O|207 ,XF0CI17 
DOUBLE PRECISION C50,C51 ,C52 ,C53 ,C54,C556 ,C56,C7 ,CTl ,C22,C33, C44 
DOUBLE PRECISION XLIVEtOlS, 207, XDEADCOiB, 0*207, XCORDC0I50) 

DOUBLE PRECISION YCOR0CO*5O7,GSMAL(O|5, 207, XLIVEfCOl57 
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14900 


DOUBI^E PRECISION C79 ,C78 ,C65, YETA , XZER ,OASHI ,DASHI1 ,D&SHI2 
DOUSltE PRECISION DirF,ONETH, fOOTH^TRP, PHI, AM, AIP,SU«l,S0M2 
D0U8l»C PRECISION BRACEl,iRACE2,8RACE3,SAI 

DOUBLE PRECISION THETR, DUMMY, THETVE,V,ANUMBR,AIW, XL1,PARI1 

DOUBLE PRECISION ZCORD(0f50 ) ,XVOLCOS5) , XDEADf (O: 5) ,XCUMCOl 503 

DOUBLE PRECISION XOl ,X02f C,DAK,X66,X667»XMAX 

INTEGER GAMMA 

ACCEPT »,MX 

ACCEPT »,NX 

C TO READ INPUT PARAMETERS FOR CONSTRAINED OPTIMIZATION 

READ(MX,») NOPT,M,K,ITHAX,IPRIN 
READCMX,»3 ALPHA, BETA, DELTA, GAMMA, 0ASHI 
READCMX,*) (XSTART(l,J),a*<l,NDPT3 
DO 100 II»2,K 

100 READ(MX,*3 CR C II ,03 ,0*1 ,NOPT) 

C TO READ INPUT PARAMETERS FOR LINEAR EQUATIONS AND 

C CONVERGENCE CRITERION, 

READtMX,*) MAXIT,NUMSIG,IPT,NUMEQ 
READCMX,*) EPS,DELNaN,PREC,ZETA 
C TO READ INPUT CONSTANTS FOR THE REACTOR, 

C 

C RATE CONSTANTS,...* 

READCMX, ♦) AKP,AKD,AKZERO,AKlP,AKTD,AMTC,AKFT,AKFM 
C MON0MER»POLfMER CHARECTARISTICS... 

READCMX, ♦) AMW,DP,DM,VM0L,SAI 
C REACTOR CHARECTARISTICS... 

READCMX, 4) THETA, TEMP, VW 
C OTHER CONSTANTS... 

READCMX, ») AI0W,AMZER0,VD,VM,CE,A«M,FEE,SWC,YETA 
C TO READ FEED CONDITIONS... 

READCMX,*) CFFDCKI3 ,KI»ia03 
DO 181 KI»0,5 

READCMX,*)CFDLIVCKI,IK),IK*1,10) 

181 CONTINUE 

DO 182 KI»0,5 

READCMX,*) CFDDEDCKI,IIC),IK*0,10) 

182 CONTINUE 

READCMX,*) CXFDCKI),K1»1,N0PT) 

ANAa6,02E+23 

RR»8,2764E+07 

C ■ ' 

0NETH»1, 0/3,0 
TW0TH«2,/3,0 

C CALCULATION OF DIMENSIONLESS CONSTANTS, 

BRACEI»C36,0*3,142)**0NETH 
8RACE2»VD**TW0TH 
BRACE3*VM**TW0TH 
C 

C2«AKIP*ANA*ANA*VM»VM*AI0W*DP/CAKP*DM) 

C3*A|CZER0*ANA*DP*BRACE3*BRACE1/(AKP*DM) 
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17900 
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18700 

18800 

18900 

19000 

19100 

19200 

19300 

19400 

19500 

19600 

19700 

19800 

19900 


C4ai(AKTC+AKT0)*DP/C2,O*AKP#0M) 

C50»AiaiilfCE/(BRACEl4BHACE3*ftMZER0*ANA) 

C51aCV«»SWC*CE)/CBRACEl»BRACEJ*AMZER0»AMA) 

C52aBRACE2/8RACB3 

C53eANA>l‘VM*OP#AMZERO/(AIOW*AMW5 

C54»AIOW»yMO0/(ANA*VD»AMZERD3 

C55SPVM /VD 

C556sAMM*VM*V«O0/VD 

C56aAKP»0M/(AKD*ANA»0Pf¥M3 

CT1»1,0 , 

C7»2,O*AKD*ANAf0P»AIQiif»VM»rEE4VW/CAKP»DM4AMZERO) 

THETR»THETAf AKP»DM/CAHA40P*VM3 

DUMMy»l,0E+O3 

C4O3*»AKTC*0P/£2,O*AKP»0«3 

C4O4=!AKTO*0P/C2,O*AKP»DM3 

fHETVEai.O/THETR 

C63aAKFT*AIOW4ANA40P»VM/CAKP4DH3 

ClOl*DP)|!VM4ANA/tDM*VMO03 

ClO2aAKFM*ANA»OP»VM/CAKP»0M«VMO03 

OU«MY2«1,OE»02 

C55»C£4,O*3,142/£3,O*VM3)**ONETH3*(2,»VMO04YETA/CRR»TEMP33 

OyMMySal.OE-OI 

0««My4*4,2985E+O3 

f « V f « « « « « y « « « i; « ^ f v # f f f i| 0 |< iK f « IK « jp; » f ^ IK « » f g|( g|^ « f in I|( « iK f IK K y j|( « JK K )K j|t f )|( J|I |f « 9)1 i|e ^ « 

C f *■ 

C . THE FO00O8IRS LINES SHOULD BE INSERTED FOR SECOND APPROACH, 

C X0l»0,0l 

C X02a0,99 

C CaO.OOOOl 

C iSRITECNX, 41893 

c 


C 


THE DATA 


NX, 3185, 

NX, 20013 
NX, 31853 
NX, 41813 
NX, 2002) 

NX, 20033 
NX, 2004) 

NX, 2005) 

NX, 2006) 

NX, 20073 
NX, 2008) 

NX, 2009) 

NX, 20103 

:nx,20ii) amw,dp,d«»vmol,sai 
:nx, 20123 


NOPT,M,K,ITMAX,IPRIN 
ALPHA, BETA, DELTA, GAMMA, DASHI 

(Cy,XSTA8T£l,J)3,Jal,NOpT) 



, MunKU 

ON,PREC,ZETA 

AKZERO , AKIP , AKTD , AKTC , AKFT , AKFM 

IM . VMnl. .SflT 
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22900 
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23100 
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23600 

23700 

23800 

23900 

24000 
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24200 

24300 

24400 

245U0 

24600 

24700 
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24900 


WRITE(NX#2014) 

WHlTECNXr20l5) AIOW,AMZERO,VD,VM,CErA«M»FEE,SWC,¥BfA 
WRITECNX,3015) 

WHITE(NX,3016) (CKI#rFB(KI)),Kl»l,10) 

WRITECNX,2113) ” 

OP 281 KlaO^S 

WRITECNX,2114) CKI,IK,rOliIV(KI,IK),IK*l,iO) 

281 CONTINUE 
WRITECNX,2115) 

DP 282 KI»0,5 

WKITECNX,2U6) (KI , IK, FD»ED(KI , IK) , IKal , 10) 

282 CONTINUE 
WRITECNX,2487) 

WRITECNX,2488) (KI,XF0(KI) ,KI»1 ,N0PT) 

c 

c 

2001 FORMATC//10X,'DATA for the problem'//) 

2002 rOR«AT(//5X,'NOPT « « ',I2,5X,‘'K * ',I2,5X, 

l^ITMAX a '',I5,5X,'IPRJN a %12) 

2003 FORMAT C/5X,'Ai:,PHA a * , riO,§, SX, •’BETA a ^ ,F10,6,5X, 'DElUTA a , 
iriO,6//SX, 'GAMMA a ' , 12 , 13 X , 'DASHI a »,FlO,6) 

2004 FORMATC//SX,'INITIAl. GUESSES-,,,.,,.*..' ) 

200$ F0RMAT(/2<:5X,'XSTARTC*',I2,') a %E1S,8)) 

2006 rORMATC//$X, 'INPUT PARAMETERS FOR LINEAR EOUATIONS AND 
ITHE CONVERGENCE CRITERION') 

2007 FORMATCXSX^'MAXIT a ' , 15 , 5X , 'NUMSIG a ',12,5X#'IPT a 
H2X5X,'NUMEQ a ',12) 

2008 FORMATC/SXr'EPS a % E15,8, 3X , 'DELNON a ',E15,8,3X, 

I'PHEC a ',E15,8/5X,'ZETA a ',E15,8) 

2009 F0RMATC/5X,'RATE C0NSTA»T5..,,,'/5X,'KP a',E15,8,5X, 

I'KD a ',E15,8,5X,'KZER0 a ' ,E15,8/5X, 'KIP a ',E15,8,5X, 

I'KTD a ',E15,8,5X,'KTC = ' >E15,8/5X, 'KFT a ',E15,8,5X, 

I'KFM a ',E15,8 ) 

2010 FPRMATC/5X,'MONOMER,POLfMER CHARECTARISTICS,,,' ) 

2011 ruRMATC/SX,'MW a ' ,E15,8,5X, 'DP a ' ,E15,8,5X, 'DM a *',£15,8/ 
15X,'VM0L a ' ,E15,8,5X, 'SAI a ',E15,8) 

2012 FaRMAT(/5X, 'REACTOR CHARECTARISTICS,.,') 

2013 F0RMATC/5X, 'THETA a ' ,E15,8,5X, 'TEMP a ',E15,8,5X, 

I'VW a '',K15,8) 

2014 F0RMATC/5X, 'OTHER CONSTANTS,,,') 

2015 rPRMATC/5X,'AlOW a ' ,E15, 8,4X , 'AMZERO a ',E15,8,4X, 'VD a 
IE15,8/5X,'VM a ' ,El5,8,4X, 'CE a ' ,El5,8,4X, 'AMM a 

1K15,8 /5X,'FEE a %Ei5,8,4X,'SWC » ',El5,8,4X,'yBTA a ',BlS,i) 

3015 FURMATC/SX,'FEED C0NDIT10NS,,.,'///5X, 'RADICAL NUMBER DISTRIBOTI 

ION,,,,.') 

3016 rnRMATC/2t5X,'FFDC',l2,') a ',riO.$)) 

2113 rOHMATC//5X,'LIVl POLTMER CHAINS....,,,,..',/) 

2114 FaRMATC/3t4X,'F0LlVC',I2,lH,,I2,') *',F8,4)) 

2115 F0RMATC//5X,'DEAD PQLXMBB CHAINS...,.,....',/) 

2116 FURMATC/314X,'FDDEDC',I2,IH,,I2,') a', £8,4)) 



25000 2487 r0RMATC/5X,''rEED COND, FOR VARIABliES OF INITIAL SET') 

25100 2488 rORMATC /3 C 5X , 12 , 4X , EU , 43 ) 

2S200 4181 F0RMATC/7X,'^DATA FOR CONSTRAINED OPTIKIZATlON,,«-»-«-_' 3 

25400 C 

25500 C THE FOLLOWING LINE SHOULD BE INSERTED FOR SECOND APPROCH, 

25600 C4189 FORMATCSX, 'OUTPUT OBTAINED B¥ EVALUATING PHI FROM EQUATION 

25700 C 4,20 

25800 C 

25900 C’*'’l‘»*******>l'**f’i'f»f f *4 

26000 C 

26100 C 

26200 CALL CONSXCNOPTiM^KOTMAX, ALPHA, BETA, GAMMA, DELTA, X 5TART,R,r, IT, 

26300 llEV2,NX,G,H,XC,IP8IN,DASHI,dNETH,TW0TH,BRACEl,BRACE2, 

26400 1BRACE3,C50,C51,CS2,C53,C54,C55,C556,C56,CTI,C7,THETR, DUMMY, XFD, 

26500 IC65,SAI3 

26600 IFCIT-ITMAX) 20,20,30 

26700 20 IFCIPRIN,ECJ,03 GO TO 4501 

26800 WRIT£(NX,2U73 PCIEV23 

26900 2117 rORMATCZ/SX, 'FINAL FUNCTION VALUE «',D18,10,//3 

27000 DO 300 a5»i,N0PT 

2tl00 NRltECNX,163a,XSTART(IEV2, J) 

27200 16 rORMATCaXr'XC'ia,') « '£20,83 

ataoo 300 continue 

27400 4501 GO TO 899 

27500 30 WRlTECNX,t73ITMAX 

27600 17 F0RMAT(2X, 'WARNING MAINlTHE NUMBER OF ITERATIONS HAS EXCEEDED' 

27700 12X,I43 

27800 C 

27900 C*********>lf****V!¥V****$*’$!f**3ft**$*ili«n‘*i$**$***$iH$M*ilti»**fli$****9*****i$!$$*tft* 

28000 C 

28100 899 DO 568 Kl=!l,NOPT 

28200 XCK13sXSTARTC1EV2,K13 

28300 568 CONTINUE 

28400 999 DASHII*DASHI 

28500 XZERaX(7) 

28600 DO 167 KIal,4 

28700 XTEMPCKI)sXCKI475 

28800 167 CONTINUE 

28900 NNON»NOPT+l 

29000 XCNNON3aCONrR(TWOTH,XTEMP,XZER3 

29100 CALL nonlincnnon,x,numsig,maxit,iprin,eps,dashii,nx,oelnon,prec, 

29200 1UNETH,TWOTH,BRACEI,BRACE2,BRACE3,C50,C51,C52,C53,C54,C55,C556, 

29300 1C56, C7,Cll,C65, DUMMY, THBTR,SAI,XFD,X663 

29400 DO mi IK«1,NN0N 

29500 ir(XCIK),LT,0,03 XCIK3s0,0 

29600 mi CONTINUE 

29700 V«Xt83*DUMMY4/XC7) 

29800 A1W«XC23 

29900 AM*X(43 



30000 


Alf»«0.0 

30100 


AMUMBR»XC7)/DIJMMY3 

30200 


r23aXCl2)»CDUMMlf4»4T«0TH)/DUMMy3 

30300 


C22»C2»V»AIP 

30400 


C33aC3/CV»»ONETH) 

30500 


C44«C4/V 

30600 


C78aC7»CV»»TWOTH)*AIW/(AM+F233 

30700 


C79aC74AM*AlW/CAM+F23) 

30800 

C 


30900 

C )|e f » )|t * )« * nt f JI f Ht )(( )(t * J|( ^ 1|; J|t ^ ^ ^ mj 5I ^ ^ Jjt |5 ^ ^ ^ JK 3|, ^ ^ SK ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 

31000 

c 


31100 

c 

INSERT THE FOIiltOWING LINES IN PLACE OF THE ENCLOSED SET OF LINES 

31200 

c 

GIVEN ABOVE FOR THE SECOND APPROACH, 

31300 

c 


31400 

C899 

DO 568 KIal,5 

31600 

c 

XCKI)sXSTARTCIEV2,KI) 

31600 

C568 

CONTINUE 

31700 

C 

X667aXSTART(IEV2r6) 

31800 

C 

DO 569 KI«7,N0PT 

31900 

c 

XCKl-lJaXSTARTCIEVarKI) 

32000 

c 

XrDCKI-l)«XFDCKI) 

32100 

C569 

CONTINUE 

32200 

C999 

DASHIlwDASHI 

32100 

C' 

XZERWXC6) 

32400 

c 

DO 167 KI»1|4 

32500 

c 

XTEMPCKI)aXCKI+63 

32600 

C167 

CONTINUE 

32700 

C 

DAKsC0NFR£0NETH,XTEMP,XZER3 

32800 

C 

CALL GOLD£XOl,XO2,C,X667,DAK,XZER,SAI,C65,0UMMy4,ONETH) 

32900 

C 

NNONsNOPT' 

33000 

C 

X(NNDN)aCONrR£TWOTH,XTEMP,XZER) 

33100 

C998 

X66*X667 

33200 

C 

GALL N0NLIN£NN0N,X,NUH5IG,MAXIT,IPRIN,EPS,DASHII,NX,DELN0N,PREC# 

33300 

C 

1ONETH,TWDTH,BRACEI,BRACE2,0RACE3,C5O,C51,C52,C53,C54,C55,C556, 

33400 

C 

lC56,C7,Cfl,C6S,DyMMT,THETR,SAl,XPD,X66) 

33500 

C 

DO im Klsi^NNON 

33600 

C 

irCX(KI),LT,0,0) XCK1)*0,0 

33700 

Cllll 

continue 

33800 

C 

XZERaX£6) 

33900 

c 

no 168 KIal|4 

34000 

c 

XTEMP£Kl)aX£KI+63 

34100 

C168 

CONTINUE 

34200 

C 

DAK«CaNFR£DNETH,XTEMP,XZER) 

34300 

c 

CALL GOLOtXOl,X02,C,XMAX,DAKrXZER,SAI,C65»OUM«T4|ONEfH) 

34400 

c 

X667»XMAX 

34500 

c 

DIFr«X667«X66 

34600 

c 

ir£IFRIN,EO,0) GO TO 4938 

34700 

c 

WR1TE£NX,4086) DIFF,X667 

34800 

C40B6 

F0RMATC/5X,'DIFFERANCE in X66 « %DIS,8/5X, "XeS *» »,Dl5.9/3 

34900 

C4938 

IFCABSCDIFF3.LT,0.0001) GO TO 987 



35000 

C 

GO fO 998 

35100 

C987 

VaX(7)*DUMM74/XC6J 

35200 

C 

AIWaXC2) 

35300 

c 

AMaXC4) . 

35400 

c 

AIPbO.O 

35500 

c 

Af|UM8RBXC6)/DUMMy3 

35600 

c 

r23aXCll)’''(0UMMy4»4TWOTH5/DUM«y3 

35700 

c 

C22*C24V»IP 

35H00 

c 

C33sC3/CV»4aNEXH) 

35900 

c 

C44«C4/V 

36000 

c 

C78«C7»CV»fXW0fH)*AIW/CAM+F23) 

36100 

G 

C7yaC7*AM»AIW/(AMtF23) 

36200 

c 


36300 

Q f icy f y f f f « f y i|e )|( |! y f f y ^ f j|e y y ^ ^ ^ It )|; f . 3|( jp ^ ^ JK ,K ^ ^ 3|, ^ ^ ^ ^ ^ ^ y ,|e ^ ^ y f ^ ^ ^ 

36400 


CAI^L IiINCNUMEO,FF0,ANUMBR,THETR,C22,C33,C78,C79,C44,NX,DASHI2, 

36500 


lNTeMP3,B,IPRIN) 

36600 


DIFr*0ASHl2«DASHIl 

36700 


lF(lPRl«,Ea,03 GO TO 4502 

36800 


WRlTECNX,4002) 0IFF 

36900 

4002 

rORMATC/2X,’'DIFF *',£15,8/3 

37000 

4502 

IFCABSC0IFF3,LT,ZKTA) GO TO 98 

37100 


0A$«I«0ASHI2 

37200 


GO TO 999 

37100 

98 

IFCiPRIM,eQ,03 60 TO 4503 

WOO 


i#RtTECNX,20l7) 

37500 

2017 

format (/*■ FUNCTION VAliUES AFTER FIRST SET*' 3 

37600 


RRITECNX,20183CCI,XCI33,Ib1,NNON3 

37700 

2018 

F0RMATC1C5X,'XC%I2,'3 « *',E15,8)3 

37800 

4503 

DU 187 IKb1,NTEMP3 

37900 


KIbIK 

38000 


XRADCKI-13bBCKI,13 

38100 

187 

CONTINUE 

38200 


IF(IPRIN,SQ,0) GO TO 4504 

38300 


WRITE(NX,2O193CCKI,XRA0£KI33»KIbO,NTEMP33 

38400 

2019 

FORMATC/rlC5X^XRADCSI2,'3 ** %E15,83) 

38500 

C 

LIVE POliYMER CHAIN LENGTH DISTRIBUTION STARTS, 

38600 

4504 

NTEMP2bnteMP3«1 

38700 


PHIsX(63 

38800 

C 


38900 

K ^ |( y K «« f |( f « K f i; f f i; It «« le f « It y ««««» f f « i«( «« |( f «* y « K |i « f » f y ]|( «««««« y «3|c f 

39000 

C 


39100 

c 

INSERT THE FOLLOWING LINE INSTEAD OP THE PREVIOUS LINE FOR 

39200 

c 

SECOND approach. 

39300 

c 

PHI»X66 

39400 

c 


39500 

(;;yyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyy 

39600 

^C 


39700 


DO 878 KI*l,NTEMP2 

39800 


XHVEC0,K13»KI»XRADCKI3 

39900 

878 

CONTINUE 



40000 

40100 

40200 

40JO0 

40400 

40&00 

40600 

40700 

40800 

40900 

41000 

41100 

41200 

41300 

41400 

41500 

41600 

417O0 

41800 

41900 

42000 

42100 

42200 

42300 

42400 

42500 

42600 

42700 

42800 

42900 

43000 

43100 

43200 

43300 

43400 

43500 

43600 

43700 

.43800 

43900 

44000 

44100 

44200 

44300 

44400 

44500 

44600 

44700 

44800 

44900 


857 


4250 

2499 

2500 
862 
C 

4505 


911 


912 

C 

C 


4251 


963 

4506 


871 

872 


no 857 MOHITsl,4 
M0Mitt;NT*M0Ji?T 


lMO«ENT,PHI,THETVE,rDLIV,XRAD,XLlVE,T,B) 

Mstl < w ■ ^ 


CALL MATIN CT,NTEMP2,B,M, DETER) 
DO 857 KI=l,flTEMP2 
X6TVECM0MENT,KI)=BCKI,1) 
continue 


DO 862 HONT*0,4 
IFCXPRIN,E0,O) GO TO 4505 
WftITECNX,4250) 

POLYMER CHAIN LENGTH DISTRIBUTION,. 
s4RTTE(NX,2499) MONT 
F0RMAT(/5X, 'MOMENT s%i2) 

WRITE(NX,2500) CKI,XLIVE(WONT,KI),KIsl rNTEMP2) 
F0RMAT(2C5X,I2,5X,E15.8)) 

CONTINUE 


DEAD POLYMER CHAIN LENGTH DISTRIBUTION STARTS, 

DO 911 MOMENT»0,4 
DO 911 IRADal ,NTEMP3 

GSMAL ( MOMENT, IRAD) =XLIVE (MOMENT, IRAD)/XRAD(IRAD) 

CONTINUE 

DO 912 MOHTa!0,4 

MOMENfaMONt 

CALL SEf3CNTEMP3,C22,C33,C78,C44,Cl02,C63,C403,C404 
ITRP, PHI, THETVE,XRA0,FDDEO,XHVE,GSMAL, MOMENT, T,B) 

Msl 


,V 




CALL MATIN(T,NTEMP3,B,M, DETER) 
no 912 Klal ,NTEMP3 


XDKADC moment, KI-1 )=BCKI,1) 

continue 


IF(lPRIN,EO,0) GO TO 4506 
WRIIE(NX,4251) 

FORMATC/SX, 'DEAD POLYMER CHAIN LENGTH DISTRIBUTION 
DO 963 MONTaO,4 
NRn'EtNX,2499) MONT 

WRITE (NX, 2500) CKI , XDEAD (MONT,KI ) ,KIaO , NTEMP2) 
CONTINUE 

no 872 MOMENT=0,4 
SUMlaO.O 

DO 871 Kl=l ,NTEMP2 
SUMiaSUMl+XLTVE(MOMENT,KI) 

CONTINUE ' 

XLIVET(MnMENT)=SUMl 

CONTINUE 

DO 874 MOMENTaO,4 
SUM2=0,0 


45000 

45100 

45200 

45300 

45400 

45500 

45600 

45700 

45800 

45900 

46000 

46100 

46200 

46300 

46400 

465O0 

46600 

46700 

46800 

46900 

47000 

47100 

47200 

47300 

41400 

47500 

47600 

47700 

47800 

47900 

48000 

48100 

48200 

48300 

48400 

48500 

48600 

48700 

48800 

48900 

49000 

49100 

49200 

1 49300 
49400 
49500 
49600 
49700 
49800 
49900 


00 873 KIsO,MTEMP2 
SUM2aSUM2+XDEAO(MOMENf ,KI) 

873 COMTIWUE ■ ■ • 

X0EADTCMOMENT)=SUM2 ' 

T!^PE »,SUM2 
87 4 COMIINIie 

A«l^wTsAMW»XDEADT( 1 3 /XDEADTCO) 

AMWWTs AMW*XOEADT(2)/XDEA0T(t) 

PDI=AMWWT/AMNWT 
?AVE»X C8 3 ♦l)UMMy4»VM/X ( 7 3 
C 

c 

DAVE2s((VAVE*6, 0/3, l423»»oifE7'H3#( 1 ,00+083 

0AV£l»CCO8FR(ONETH,XTEMP,XZEB3*(DUM«Y4»*O*^ETH3*C(6,0*VM/3,l42 3 

l+*08ETH)/XC73)»l,000+08 

PATNU(4«XC7 3 ♦AMZERQ*AMA/0UMMY3 

c 

c . ■ 

C INSERT THE FOLLOWING LINES IN PLACE OF THE SET OF T 4 INES 

C ■ ENCLOSED ABOVE FOR THE SECOND APPROACH, 

G VAVE*X(7)«D0MKY4*VM/X(63 

C 0AVE2«(CVAVE»6, 0/3, 1423»*0NETH)4(1. 00+083 

C OAVE1*CCONFRCONETH,XTENP,XZER)*COUMMY4*»ONETH)*C(6,0»VM/3,142) 

C l**0NETH3/X(6)3»l,00D+08 

C PATNU«bXC63»AHZERO=«‘ANA/DUN'MY3 

c 

C • 

C0NVER=(XFDC33-X(333/XFD(3) 

RATE= ( XFO ( 3 3 -X( 3 3 3 *AIOW43 ,6E+06/ (THETA*DUMMY2 3 
C 

c : • ' 

c 

c ■ 

C nuTPUT section BEGINS— 

WRIXE(NX,3186) 

318b FORMAT(lKl) 

WRITE(NX,3185) 

3185 FORMAT C IX » ' 

1 '3 
WRITECNX, 21473 

2147 FORMATC/lOX, 'OUTPUT FROM THE PROGRAM— «_—«—*/ 3 
WHITECNX, 31853 

WRITE(NXi 21483 PATNUM,VAVE,DAVE2rDAVEl,COMVER,RATE 

2148 FORMAT(/5X,'(U NUMBER OF PARTICLES C/CC EMULSION)*— »%EIS, 8 / 
15X|'(23 AVERAGE VOLUME OF PARTICLES CCC)— *— *%E15,8/ 
15X,'C33 AVERAGE DIA ? VOLUME BASIS ( "A)*_** ,E15 , 8/ 



60000 
50100 
50200 
50300 
50400 
50500 
SOfoOU 
50700 
50800 
5oyoo 
51000 
51 100 
51200 
51300 
51400 
51500 
51600 
51700 
51800 
51900 
52000 
52100 
52200 
52300 
52400 
52SOO 
52600 
52700 
528O0 
52900 
53000 
53100 
53200 
53300 
53400 
53500 
53600 
53700 
53800 
53900 
54000 
: 54100 
i 54200 
• 54300 
i 54400 
i 54500 
V 54600 
54700 
I 54800 
54900 


15X,'(4) AVERAGE DIA ;NU«BER BASIS f 

15X,'(5) CONVERSION obtained «-,.%E15.8/ 

15X,'(6) RATE (MOOES/ (LIT. HOORl)-..^ ,*%E15,8/) 

WRITE(NX,2149) AMNWT,AMWWT,PD1 

2149 FORMATCSX, '(7) NUMBER AVERAGE MOLECULAR WEIGHT., %E15,8/ 

15X,^C8) WEIGHT AVERAGE MOLECULAR WEIGHT fElS.R/ 

15X,'C9) POLIDISPERSITY INDEX™..— E15 . 8// J 
WRITF.CNX,318S) 

IFdPT.EQ.O) GO TO 4510 
WftITECNX,3186) 

WRITECNX,3187) 

3187 FURMAK /5X, 'DETAJ^LS OF THE DISTRIBUTIONS FOLLOW. V) 

WRITKCNX,3185) 

C 

c 

IWHICHSO 

XL1»1.0 

PART1»5000,0 

CALL DISTRCXLl , PARTI ,XLIVET.IWHICH.XCORD,ICORD,ZCORD,XCU«3 
WRITECNX, 21503 

2150 F0RMATC/5X, 'LIVE POLYMER CHAIN LENGTH DISTRIBUTION'/! 
WRITECNX,3185) 

, WRITBCNX, 21513 

2151 FORMATC/5X,'NO',5X,'MOL. WT. r,5X, 'NUMB. BASIS' ,5X, *WT, BASIS 
ICNORMAL)'/) 

SUMI=0.0 

PART1*XC0RD(2)-XCQRD(13 
DO 870 KIaO.50 
SUMlaSUMl+ZCORD(KI)»PABTl 
870 CONTINUE 

DO 863 XI=0,50 
ZCOR0(KI)aZCaBD(KI)/SUMl 
863 CDNIINUE 

DO 879 KI=0,50 

WRITECNX, 21523 KI ,XCORD(KI) .YCORDCKI) .ZCORDCKI) 

2152 FDRMAT(5X.I2.3X,EU,4,3X,Eli,4.3X,Ell,4) 

879 CONIINUE 

c 

c 

c 

CALL DISTR(XH .RARTl ,XDEADT,IWHICH.XCORO. YCORD.ZCORD.XCUM) 

SUMlxO.O 

PABTl»XCORD(2)-XCORD(l) 

DO 4851 KI»0.50 
SUMl=SUMl+ZCnRD(Kl) ♦PARTI 
4851 CUMTInUE 

DO 4853 KlaO.SO 
7X’nRDCKn»ZC0RD(KI)/SUMl 
4853 continue 

WKrTEtNX^3186) 



55000 

55100 

55200 

55300 

55400 

55500 

55600 

55700 

55800 

55900 

56000 

56100 

56200 

56300 

56400 

56500 

56600 

56700 

56800 

56900 

57000 

57100 

57200 

57300 

57400 

S7S00 

57600 

57700 

57800 

57900 

58000 

58100 

58200 

58300 

58400 

58500 

58600 

58700 

58800 

58900 

59000 

59100 

59200 

59300 

59400 

59500 

59600 

59700 

59800 

59900 


WRlTE(NXr-n88) 

;«!8TTE(NX,3185) 

3188 FORMAf C //) 

WRX^E(f^X,2l53) 

2153 FORM4T(/SX, *^DEAD POL7MER CHAIR LENGTH DISTRIBUTION V) 
WRIXECNX,3185) 

WRITECNX,2154) 

2154 FORMAT C/5X, »NO*,5X, *MOL,WT. » ,6X, 'NUMB.BASIS' ,5X, "NT. BASIS 
1 (NORMAL)*/) 

no 891 KIa0,S0 

WRXTECNX,2152) KI , XCORD (KI ) ,yCORD (KI ) ,ZCnRD(KI ) 

891 CONTINUE 

C . , ■■■■■ ■ 

c 

XLtel.O 

parti«ioo,o 

INHICH=l 
00 892 KI«0,4 

XVOLCKr)»X(KI+7)l'(DUMM74»»KI)/I)UMMyi 

892 CONTINUE 

CALL OISTRCXLl, PARTI, XVOL»lWHICH,XCORU,yCORD,ZCORD,XCUM) 

DO 2176 KI«0,50 

• XCORO(KI)»(CXCORD(KI)*VM»6,0/3.142)»*ONETH)*C1.0D+08) 

2176 CONTINUE 

6fRlTE(NX,3186) 

WRITE(NX,3188)' 

WftITE(NX,3185) 

WHITE(NX,2155) 

2155 FORHATC/5X, 'PARTICLE SIZE DISTRIBUTION (CUMULATIVE)*/) 
^RnE(NX,3185) 

!»IRTTE(NX,2156) 

2156 FaRMAT(/5X,*NO*,7X, *DIA("A)',9X, 'FRACTION* /) 

DO 890 KI=Q,50 

WRTTE(NX,2157) KI , XCORD (KI ) , XCUMCKI ) 

890 continue 

2157 F0RMAT(5X,I2,5X,E1U4,6X,E11,4) 

c 

4510 WRirE(NX,3186) 

STOP 

END 

C 

c 

SUBROUTINE NONLIN (N ,X ,NUMSIG,MAXIT ,IPR1N , EPS,DASHI , NX,DELNON , 

iPREC,0NETH,TWnTH,BKACEl,BRACE2,BRACE3,C50,C51,C52,C53,C54,C55, 

lC556,CS6,C7,C’n,C6S,0UMMI,THETR,SAI,XFa,X66) 

c 

C THIS SUBROUTINE SOLVES A SET OF NONLINEAR ALGEBRAIC 

C EOUATIONS.IT IS USED TO SOLVE THE EQUATIONS OF FIRST SUBGROUP 

C OF GROUP 1, OTHER SUBROUTINES DIRECTLY CONNECTED WITH THIS 



60000 

60100 

60200 

60300 

60400 

60500 

60600 

60700 

60800 

60900 

61000 

61100 

61200 

61300 

61400 

61500 

61600 

61700 

61800 

61900 

62000 

62100 

62200 

62i00 

6240'0 

42500 

62400 

62700 

62800 

62900 

63000 

63100 

63200 

63300 

63400 

63500 

63600 

63700 

63800 

63900 

64000 

64100 

64200 

64300 

64400 

64500 

64600 

64700 

64800 

64900 


subroutine are BACK,AMAXl,AMINf,FUNX. 

SEE REF 31 FOR DETAILS, 

DISCRIPTION OF PRINCIPAL ARGUMENTS-, 

XCD— — — — — — — HItHE vector of VARIABLES. 

NUMSIG— — — — ^^..-CONVERGENCE CRITERION . 

EPS— ——————..CONVERGENCE CRITERION, 


—MAXIMUM NUMBER OF ITERATIONS ALLOWED BY 
THIS PROCEDURE. 

—MINIMUM ALLOWABLE STEP SIZE FOR CALCULATING 
PARTIAL DERIVATIVES, 

—MINIMUM ALLOWABLE ABSOLUTE VALUE OF PARTIAL 
DERIVATIVES, 


C 
c 
c 
c 
c 
c 
c 
c 

C MAXIT- 

c 

C PREC— . 

c 

c delnon. 

c 

c 

DOUBLE PRECISION X C 123 ,PARTC1 2) r TEMPI 1 2 ) ,COEC 12 , 1 3 ) , XFDC 1 1) 
DOUBLE PRECISION FMAX,FPLUS,F, ABSFF, HOLD, FACTOR, ETA, H 
DOUBLE PRECISION ABB , DERMAX , TEST . AA , CC , X66 , DASHI , QNETH , XWOTH 
DOUBLE PRECISION BRACEl ,BRACE2,BRACE3 ,TaETR, DUMMY, SAI 
DOUBLE PRECISION C50 ,C5l ,C52 ,C53 ,C54 ,C55 ,C556 , C56 ,C7 , CTl ,C65 
DIMENSION ISUBC123,LOKP(12,12) 

RELCOMslO,OE+00*»(-NUM5lG3 

UTESfst 

00 700 «*1,MAXIT 
DO 1111 KI»1,N 
IK=KI 

IF(XCIX),LE,0.03 XClR)sO.O 
till CONTINUE 

L0UIT=0 
FMAX=U.O 
M1=M-1 

IFdPRIN.NE.l) GO TO 9 
WRlXEtNX,49)Ml,(XCI),I=i,N) 

49 FORMATC5X, 'Ml= * , 15/5 C5X,E1 1 ,4) ) 

9 DO 10 J=1 ,N 

10 LOKPfl,J)sJ 
DU 5555 K=1,N 
rF(K-l3 134,134,131 

131 KMINsK-I 

KTEMPsK 

CALL BACK CKMIN,N,X, ISUB,COE,LOKP,KTEMP) 

134 JJJSK 

CALL FUNXCX,r,JUJ, DASHI, 

lONETH,IWaTH,BRACEl,BRACE2,BRACE3,C50,C51 ,C52 ,C53 ,C54 ,C55 ,C556 , 
1C56,C7,CT1 ,C65, DUMMY, THETR, SAI, XFD,X66) 

ABSFFsABSCP) 

FMAXaAMAXlC,FMAX,ABSFF3 



65000 

65100 

65200 

65i00 

65400 

65500 

65600 

65700 

65800 

65900 

66000 

66100 

66200 

66300 

66400 

66500 

66600 

66700 

66800 

66900 

67000 

67100 

67200 

67100 

67400 

67500 

67600 

67700 

67800 

67900 

68000 

68100 

68200 

68300 

68400 

68500 

68600 

68700 

68800 

68900 

69000 

69100 

69200 

69300 

69400 

69500 

69600 

69700 

69800 

69900 


C 

c 

rrCABSCt'),GE.EPS) GO TO 1345 
OOUITkEOUXT+I 
If (LOUIT.WE.N) GO TO 1345 
GO 10 725 

1345 FACTOKaO.OOlE+00 

135 TTALLlfsO 

DO 200 laK,N 
nEMP*LOKPCK,I) 

H0I.D*XC1TEMP) 

ET&«FACTOR»ABSCHOLD) 

HaAMIMlCF«AX,ETA) 

CCC H»ETA •. 

TF(H,LT,PREC) HsPREC 
XCITEMP)sHOOD+H 
IF(K-U 161,161,151 
KTEMPsK 

151 GALE BACKCKMIN,H,X, ISUB, COE, EOKP,KTEMP) 

161 JJjssK 

CALO FU«X(X,FPLyS,Jaa,DASHI, 

10NETH,TW0TH,B8ACE1 ,BRACE2,BRACE3,C50,C51 ,C52,C53 ,C54,CS5,r556, 
1C5§,C7,CT1,C65, DUMMY, THETR,SAI,XrD,X66) 

PARrClTEMP)«CFPEUS-F)/H 

XCITEMP)*flOl.D 

IFCABSCPARTCITEMPD.LT, DELTA) GO TO 190 
C 

ABB=A0S(F/PART(ITEMP)3 
IF(AB6.LE,1 ,E+15) GO TO 200 

190 TTALL¥=ITALLY+1 

200 COMTIfeiUE 

C 

MK~N-X 

TF(1TALLY,LB.NK) GO TO 202 
C 

FACTnR=FACTOR*10,OE+00 
TF(FACTOR,Giai.) GO TO 775 
GO XO 135 

202 TF(K,LT,N) GO TO 203 
IFCABSCPARTCITEMP)) ,LT, DELTA) GO TO 775 
NpLlaN+l 

CQE(K,MPLns0,0 
F(ViAX*ITEMP 
GO TO 500 

203 KMAX*LOKP(K,K) 

DERMAXsABSCPARTCKMAX) ) 

KpLUSaK+l 

Dp aiO XaKPLUS,N 
JSUBafiOKP(K,I) 



70000 


Tys'r*4ftscyART(asiJB)) 

70100 


TrCIKST,LT*DERMAX3 GO fO 209 

70200 


Ofe:R«AX=TEST 

70i00 


f,aKP(KPLU5,I)sKMAX 

70400 


KMAXSJSUB 

70500 


GU ft) 210 

70600 

209 

GOKPCKPLUSa)=JSWB 

70700 

210 

CONl'INUE 

70800 


irCABS(PART(KMAX)),EQ,0,0) GO TO 775 

70900 


1S(JBCK)»KMAX 

71000 


C£JE CK,M+l)sO,0 

71100 


DO 220 JB((PLUS,N 

71200 


asoB=LOK:PCKPLus,a) 

71300 


COECK,JS!J83**PARTCJSU0)/PARTCKRAX) 

71400 

C 

71500 


NPLlsM+1 

71600 


cobck,npli)=!Cok(k;,nplu+partcjsor)*x(jsub) 

71700 

220 

CONflMUE 

71800 


NPLl=Mtl 

71900 

500 

CDE CK, 8PI<n a CCOBCK, NPLl)-r) /PART (KMAX)+X(KMAX) 

72000 

C 

|«RITE(NXr6666) K, CISUBCIK) , IKsJ , 1 2 ) , ( (l*OKP( I J, JI ) , I Jsl , 12 3 

72100 

c 

1JI»1,123 

72200 

C6666 

P0RMATCX2/5(2XrB15,7)/5(2X,El5.7) J 

72100 


CONXIflOE 

72400 


?IPBl»it||+l ' 

tasoo 


X(KMAX)*COECMrMPi:iU 

72600 


IFCM.EQ.U GO TO 610 

72700 


Wl=M-i 

72800 


CALL BACK (N1,M,X,ISU8, COE, LOKP,K) 

72900 

610 

IFCM-U 650,650,625 

73000 

C 


73100 

625 

DO 630 1 = 1, fl 

73200 


AA=ABS(XEMPCI)-X(I)) 

73300 


CC=ABS(XCl))*RELCON 

7 3400 
73500 


TF(AA.GT,CC) GO TO 649 

73600 

630 

COOTTNUE 

73700 

C 


73800 

c 


73900 


-ttest=jtestfi 

74000 


lFCJTEST-33 650,725,725 

74100 

649 

OTESTsl 

74200 

650 

DU 660 1=1 ,M 

74300 

660 

TEMPCDaXCI) 

74400 

700 

CONTTWtIE 

74500 


•IP.TTEIMX, 1753) 

74600 

1753 

FORMATC/2X, 'WARNING NONLINlNO CONVERGENCE MAXIMOM NO OF 

74700 


UTERATION USED') 

74800 


TFCiPRIN,NE,l) GO TO 800 

74900 


WHITRCMX, 1763) 



7500Q 1763 FOR«fiT(/'riINCTION VALUE AT LAST APROXIMATIOM FOLLOW’) 

75100 TFLAGaJ 

75200 GO TO 7777 

75300 725 IF C iPRIN . NE. I ) GO TO 800 

75400 TFLAGsO 

75500 7777 DO 750 lKsl,N 

75600 KITalK 

75700 

75800 CALL f UNX(X,PARTCKri3 ,KII ,DASHI , 

75900 lUMETH,TWOTH,BRACl£l,BRACE2,8PACE3,C50,C5l,C52,C53,C54,C55rC556, 

76000 1CS6,C7,CT1,C65, DUMMY, THETR,SAI,XFDrX66) 

76100 750 CONIINUE 

76200 IFCiFLAG.ME.U GO TO 8777 

76300 MHTXECWX, 7788 ) (PARTCK) ,Ksl ,1) 

76400 7788 FORMATC3F15.8) 

76500 GO ID 800 

76600 8777 WRITECNX, 751) 

76700 751 F0RMAT(//5X, ’CONVERGENCE HAS BEEN ACHIVED .THE FUNCTI’) 

76800 WRITECNX, 7515) (PARTCK) ,K=1 .N) 

76900 7515 F0RMAT(5X, ’VALUES AT FINAL APPROXIMATION FOLLOW*// 

77000 1C2X,F16,8,3(2X,F16,8))) 

77100 MRITECNX, 8888) (XCKI) ,Klsl ,N3 

77200 @088 F0RMATC5X, ’ANSWER ’/3 C2X,F16,8) ) 

77300 ‘ GO TO 800 

77400 775 WftITE<NX,752) 

ffSOO 752 F0RMATC//5 X, 'warning NONLINjMODIFIED JACOBIAN IS SINGULAR TRY *) 
77600 WRirE(NX,7525) 

77700 7525 FORMAKSX, 'DIFFERENT INITIAL APPRUXIMATIOS ’) 

77800 DO 4005 Isl ,N 

77900 nPE 4500, (COE(I,J) ,J=1 ,N) 

78000 4005 CUNIInUE 

78100 TYPE » , (ISUB(I) ,Ial ,N) 

78200 4500 FURmAfC / , 3X, 6 (2X,E1 I . 4) ) 

78300 C max1T=M+1 

78400 800 return 

78500 END 

78600 C 

787 00 C 

7880U SURHOUTlNfi BACKCKMIN , N , X , ISUB,C0E ,LUKP,K) 

78900 c*^**************'********^***^***************************************** 

7 9 000 C ' ' 

79100 C THIS SUBROUTINE IS CALLED FROM NONLIN.IT IS USED FOR 

79200 C BACKSOLVIWG THE SETS OF EQUATIONS OBTAINED IN NONLIN AT VARIOUS 

79300 C STAGES OF ITERATION, 

7940U C , 

79500 *41 ******* **$*******$*******$***^**$**********************$ 

79600 Duimr.E PRECISION XC12) ,C0E(12,13) 

79700 DIMENSION ISUBC 12 ) ,LOKP( 1 2 , 12 ) 

79800 DO 200 KKsl ,KMIN 

79900 KMaKMrN«KKf2 . . 



80000 

80100 

80200 

80300 

80400 

80500 

80600 

80700 

80800 

80900 

81000 

8U00 

81200 

81300 

81400 

81500 

aifooo 

81700 

8J800 

81900 

82000 

82100 

82200 

82100 

82400 

82500 

82600 

82700 

82800 

82900 

83000 

8310U 

83200 

83300 

83400 

83500 

83600 

83700 

83800 

83900 

84000 

84100 

84200 

84300 

84400 

84500 

84600 

84700 

84800 

84900 


C ' 

KM1»KI4-1 
Ki^AXalSUBCKMl J 

K(»fi4Ax)s0,0 
DO 100 JJSKM,N 
JSU6=L0KPCKM,iJ) 

KMl sKi4""l 

X(KMAX)»X(KMAX)+C0E(KM1,JSUB)*X(JSIJB3 
100 coMtis^ue 
KMIsKM-1 
NPLis?< + l 

X(KMAX)=X(KMAX)+COECKMi,NPLl> 

200 COiaTINUE 

RETURN 
END 
C 

c 

FUNCTION AMAXKA^C) 

c 

C THIS SUBROUTINE IS CALI-ED FROM NONl<IN,lT REPORTS THE 

C GREATER OF THE TWO NUMBERS REPRESENTING THE ARGUMENTS C I ,E, A AND 

C C3. 

C 

DOUBLE PRECISION A,C 
1E(A,GT,C) C=A 
AMAXisC 
RETURN 
END 
C 

c 

: function AMINl CP,«) 

c 

C THTvS subroutine IS CALLED FROM NONLIN.IT REPORTS THE 

c smaller of the two numbers representing the arguments(I,e,p and 

C 0). 

c 

nouHr.F,: pretision p,g 

TF fP.GT.O) F=0 • 

AMTN1 aP 

PETURN ; 'i ^ 

END 

c ■ ■ : " 

c * 

SUBROUTINE FUNXCX, RES, K,DASHI, 

1UNETH,TW0TH,HRACE1 ,BRACE2, BRACES, C50,C51,C52,C53, C54, CSS, C556, 



85000 

85100 

85200 

85300 

854O0 

855O0 

85500 

85700 

85800 

bsyoo 

86000 

B6100 

86200 

86300 

86400 

86500 

86b00 

86700 

86800 

86900 

87000 

87100 

87200 

87300 

87400 

87500 

87600 

87/00 

87800 

87900 

88000 

88100 

88200 

88300 

88400 

88500 

88600 

88700 

88800 

88900 

89000 

89100 

892O0 

89300 

89400 

89500 

896O0 

89700 

89800 

89900 


tC56,C7,CTl ,C65, DUMMY, THETR,SAI,XFD3 

c 

C THIS SUBPOUTINE IS CALLED FROM NOhiLIU.IT REPORTS THE 

C VALUE OF K TH EyUATlON OF THE SET OF ALGEBRAIC KOUAITOMS TO 

C HE SOLVED BY NQfSiLIN, 

C 

(',i^f.!lt^i^^t**********$3^^*if$t*>^*$lli*^ilf:*i^1^il/.ilf:ii*$ilfil^t***^********************it^*** 
DOimLE PRECISION XC 1 2 ) ,XFD( 12) ,F( 6) ,FFD (6) ,FnLD(6 ) 
double precision RES,DASHI,ONETri,TWOTH 

double precision bracei,brace2,brace3, dummy, thftr, sat 

DOtlSLE PRECISION C50 ,C5 1 ,C52 ,C53 ,C54 ,C55 ,C556 ,C56 ,C7 ,CT1 ,C.65 
C DOUBLE PRECISION rZERO,FFDEEP 

C COMMON UN£TH,fMOTH,BRACEl,8RACE2,BRACE3,C50,C51,C52,C51,C54,C55 

e , lC556,CS6,CTi ,C7,THETR»aUMMY 
NMAX«4 

c 

c 

c 

S!/1«8,08b-04 
M2*47,69 
i3»3,807E-02 
«4*2%937E+06 
' . ■' H5af7.0998E»0l 

«6»0,63 
W7»0,l99r 
F0T4D(1)58.5111E+02 
F0LD(2)37,3E+06 
F0LD(3)»9,432E+10 
FOLD(4)=t.622E+15 
C 

DUMMYlsl.OE+06 
DUMMY2 =s 1 ,0E-02 
DUMMY3=1*OE«03 
DUMMY4=4,29B5Ef03 
C 

c X ARE THE OUTPUT VARIABLES FROM THE REACTOR 

C XFD ARE THE FEED CHARECTARISTIC OF AFORESAID VARIABLES 

FZEHOsXCV) 

FFnZF,HsXFl)(7 3 
no 210 ENtEMPsl ,MMAX 

ntfmsintemp 

FCUTF.M) = XCNTEM+73 
FFr)CNTEM)sXF0(NTJBM+7) 

210 COUTTNUE ' ■ , 

Tf (K.GT.7) GO TO 205 
r,u 10 (1,2,3,4,5,6,7)*K 

1 CC."‘iH=C0?4rRCMuTH,F,FZERQ) , 

RESaXC12)-CUNh ' 



90000 


RETURi^ 

90100 

2 

R£S« ( XED ( U -XC 1 ) ) /THETR 

90200 


RESsRES/Wl 

90300 


RETURi'i 

90400 

3 

RES=C56*CXFDC2)*X(2))/THETR-XC2) 

90500 


RES=RE5/W2 

90500 


RETURM 

90700 

4 : 

ReS=(XFP(3)-‘X{3))/{THETR*0UMMy2)-C5i't'XC5)*(rZERO/DlJMMX3)*DASHI 

90800 


RESsRES/Wi 

90900 


RETURN 

91000 

5 

RKS«C504X(l)-XC4)-CSl-X(12)*(OUMMf4*#TNOTH)/ 

91100 


lDlI«My3-C52*X(5) 

91200 


RES=RES/W4 

91300 


RETURN 

91400 

6 

RESeC544XC3)/DUMMy2»£XC5))-C55»(F(l)*DU«MY4/DUMMy33* 

91500 


1XC6)-C556*XC4) 

9lfo00 


RES»RES/W5 

91700 


RETURN 

91800 

7 

RES«X£6)-0,615 

91900 


ReS»RES/W6 

92000 


RETURN 

92100 

205 

IF £K,GT*8) go to 206 

92200 


RES»07»XC23*XC4)/(X(43 + CX(12)*(OlJMMy4**Ti«OTH)/ 

92300 


lUimWYS J J + CFFOZER-FZERO? / iTHEf R*DaMMY3 3 

92400 

ccc 

RES*RES/W7 

92500 


RETURN 

926O0 

205 

TFCK.GT.9)GO TO 207 

97700 


RE5=C7*XC2)*XC4)/CX(4)+CXC123*fDUMMY4**T^OTH)/ 

978O0 


1UUMMY3j)+CT 1 »DASHI=»N*X(6)»FZERQ/tDUMMY3*(l-X(5))) + 

97900 


1 tFFD(N)-F(N))41DUMMY4»*N3/tTHETR»DUMMY3) 

93U00 

CCC 

RES=RES/FOLO£N) 

93100 


RETURN 

9 3200 

207 

RES=C7*XC2)!«‘Xt4)/CX£43t(X(l2)»{DUMMY4*»TWOTH)/ 

93300 


10UHMY31 )+CTl*OASHl*N*X(6)*F(N-13*{UUHMY44=t'CN-l))/(DUMMY3» 

93400 


1 ( 1 -X ( 6 ) ) ) + C FI* 0 C N ) -F C N ) ) ♦ ( DUMMY4««' =*'fn / ( TriETR*DUMMY3 ) 

93500 

CCT 

RES=PES/F0LD(K) 

93500 


RETURN 

93700 


FNn 

93800 

C 


93900 

c 


94000 


function CONFR£FF,yY,yYO) 

94100 


4; « 4. « 4; f ^ 4: t )|i :|C 1); « f; f f f f f f f f f f )|i: fi: « ^ « f 4( t f f f JK » 1; ^ f 4; « « «: i|( « $ « 

94200 

C 


94300 

C 

THE INPUT FOR THIS SUBROUTINE ARE YY AND YYO.YYCD 

94400 

c 

REPHESFNT THE MOMENTS OF THE DISTRIBUTION ( I«l , 2 J .YYO 

94500 

c 

REPRESENTS ZEROTH MOMENT OF THE DISTRIBUTION .THE FUNCTION 

94500 

c 

RETURNS THE FF TH MOMENT OF THE DISTRIBUTION .FOR DETAILS OF 

94700 

c 

THE procedure. F.MPliOYRU SEE CHAPTER 4, 

i 94800 

c 


94900 




95000 


nuUBOe PRECJSIOM Y(6),YYC6),CCN(6J,SSN(6) rCAOPHACS) ,SALPHAC6) 

95100 


DDOBOE PRECISION YYO , AA, A, ALPHA , ALPHAl ,ALKKA,ALPH.41 , ALPHPl 

95200 


DOUBLE PRECISION ALPHFJ^ALPHJl rf'F, RATIO! 

95300 


NmAX»4 

95400 


AAs!(Y(2)=i'YYO/CYY(U*YY(l)3 ' 

95500 


YCU»l,0/(AA-1.0) 

95600 


yc2)*AA*yci3nci) 

95700 


ALPHAsYCU-l, 

95800 


ALPHAlsfCl) 

95900 


TF(ALPHAl.LT.O) GO TO 450 

96000 

c 


961 00 

c 


96200 


DO 118 1K®3,NMAX 

96300 


K I ® I K 

96400 


y(K13 sCCALPHAl*yyO/YY(t))*»KI)»YY(KI)/YYO 

96500 

118 

CONTINUE ■ 

96600 


ALPHFlsALPHA+rF+1,0 

96700 

C 


96800 


DO 100 NMMM«3,NHAX 

96900 


NaNMMM 

97000 


CALPHACNjal 

97100 


DO 100 KKJJ»1,N 

97200 


KK»KKkIJ^;* 

9IJ0O 


nkic«»n*kx 

ft400 


ALKKAaAliPHA+KK+1. 

ttsoo 


CCNCKK) = CC*1)»»KK)*FACT0(N)»GAMM<AT,PHAI )/tFACTOCKK)>l' 

97600 


lFACTa(NKK3*GANlACAl.KKA) ) 

97/00 


CALPHACNlsCALpHACNl+CCNlKKl^ Y(KK) 

97BOO 

100 

CONTINUE 

97900 


no 201 NMMM=3 ,NMAX 

98000 



98100 


SAI.PHACN) = 0,0 

98200 


AIUPriNl=N + ALPHA + l, 

98300 


SAf.PHACN)sGAM^(ALPHNl)*GAMMCALPHFl)/(FACT0CN)»GAMMCALPHA13) 

98400 


DU 201 JJKsl ■ 

98500 


iJ J s J J K 

98600 


ALPHFJ=ALPHA+FF+JJ+i ,0 

98700 


AbPttJl=Ar,PHA+JJ + l ,0 

988O0 


Al.PriNlKAl,PHA+N+l,0 

98900 


NUJaN-.JJ 

OOlOO 


SSM C J J ) = ( C - 1 . 0 ) 4* J J ) ♦GANM ( ALPHN 1 3 ♦GAMM C ALPHF J) / (FACTO ( N J J 3 ♦ 

00200 


1FACTU(JJ)»GAHM(ALPHJ133 

OOJOo 


SAT.PHA(N) = SA0PHACN)+SSN(JJ3 ■ 

00400 

00500 

00600 

201 

CUNTINUE 

00700 


AsGAMmC ALPHF 13 x ■ 

00800 


DO 400 NMWMs3,NMAX 

00900 


NaNMMM 



01000 
01 100 
ouou 

01300 
0140U 
01500 
OloOO 
01 700 
01800 
01900 
02000 
02100 
02200 
023O0 
02400 
02500 
02600 
02700 
02800 
02900 
03000 
03100 
03200 
03300 
03400 
OlSOO 
03600 
03700 
03800 
03900 
04000 
04100 
04200 
04300 
04400 
045 00 
04600 
04/00 
04800 
04900 
05000 
051(H) 

: 05200 

•; 05300 

i; 05400 
05500 
0560U 
05700 
05800 
059O0 


4k&+CALI>HAC!<)^SALPHA(N) 

400 CUMlTMipj 

A«4/GAMM(ALPHAU 

COMFRsC C¥Z (1 )/ ClfyO’f'ALPHAl ) )»4FF)4A*irY0 

RETURN 

450 RATiaiayyCU/YYO 

COHFRaYyO«CRATI01*#FF) 

RETURN 

END 

C 

c 

FUNCTION FACTO(N) 

C$ t $*^$!^***** t**^^***^***>^f tilfistf *t*>t^*!¥***** ********** t*** 

c 

C THIS FUNCTION REPORTS THK VALUE OF FACTORIAL OF AN 

C INTEGER VARIABLE N. 

C 

c********** ******** *************************************************** 

FACTOal 

IFCN.EO.O) GO TO 2 
00 i OKscl,N 
rACTOairACTO*JK 
I, CONTINUE . 

2 RETURN 

END ' 

C 

C • 

function GAMMCXX) 

c********************************************************************* 

C ■ ' . : 

C THIS FUNCTION REPORTS THE VALUE Of GAMMA FUNCTION OF 

C THE VARIABLE XX, 

C 

c.* *************************************** ****************************** 
OOUBLE; PRECISION XX,X,¥ 

TFCXX-5/,0) 6,6,4 

4 7b;R=2 

o:amm=i ,e+3Q 
RETURN 

6 XsXX 

ERRal ,0E-06 

Te;R«0 

GAMrtsil.O 

rF(A«2,0)50,50,15 
10 lK(X-2,0) 110,110,15 

15 ' XaX-1 , 

GaMMsGAmM 4!X 
Gi) 1:0 10 

SO TF(X-l,01 60,120,110 



06U00 

c 

RKG IF X IS near negative INTEGER OR ^ERO 

OfilOu 

60 

TF(X-ERR) 62,62f80 

06200 

62 ' 

Y=FLnATCXNT(X))-X 

06JOO 


IF tABS(y)-ERH) 130,130,64 ' 

0640u 

6 4 

1FC1,0-X-ERR) 130,130,70 

06500 

C. 

X NUT NEAR A NEGATIVE INTEGER OB ZERO 

06600 

70 

IF(X-l,0) 80f8O,110 

06700 

80 

GAMhaGAMM/X’ 

06800 


x*x+l ,0 

06900 


GO TO 70 

07000 

110 

YaX-1,0 

07100 


GY»1.0+Y*C-0. 577 1017+Y*C+0,9858540tY* (-0,876421 8+ 

07200 


1 C+y.83282t2+Y*C-0,5684729+Y*(0,2548205+Y+(-, 05149930)5)1)1) 

07300 


GAMrtaGAHM*GY 

07400 

129 

RETURN 

07500 

130 

lERal '■ ' : 

07600 


return 

07700 


END 

07800 

C 


07900 

c 


08000 


SORROUTINB SETl (N,FFD,ANUMBR,THETR,C22,C33,C78,C79,C44,T,B) 

08100 


08200 

c 


08300 

c 

THIS SUBROUTINE CALCULATES THE ELEMENTS OF COEFFICIENT 

08400 

c 

MATRIX OF THE SET LINEAR EOUATIONS IN SECOND SUBGROUP OF THE 

08500 

c 

GROUP 1 .COEFFICIENT MATRIX IS REPRESENTED BY T,B(I,1)IS THE 

08600 

c 

VECTOR OF CUMSTANTS , 

0870U 

c 


08BOU 

c******* 


08900 


noUHLE PRECISION T ( 25 , 25 ) ,8(25 , 1),FFDCOX20) 

0900U 


nuUBt-E PRECISION ANUHBR,THETF ,C22 ,r33 ,C78 ,C79 ,C44 ,THETVE, BRACES 

09100 


MlaN-l 

09200 


TMF,i'VE=l ./THETR 

09300 


DU 104 IROWsl ,N 

09400 


no 104 IC0L=1 ,N 

095OO 


TCIRrjwaCOL)=0,0 

09f>no 

104 

Cn‘'l fwOpj 

09700 

c 


09WOU 


IIU J05 lRUwi=l,N 

09900 


II- (iROt^.EQ.N) 00 X0 433 

lOoOO 


TCnLaiRQW 

10100 


lCnLM2a|CUl,.-2 

10200 


TCnbMl=lCOL-t 

10300 


TCOLPiaiCOLtl 

10400 


IC06P2=iC0L+2 

loboo 

c 


10600 


NESTalCQLM2 

10700 


TrrNFST.LE.-l) GO TO 102 

10800 


TF(NFST.!,E.O) GO TO 103 

10900 


T(IRnN,MFST)aC22 



IJOOO 

103 

NES'iaiCDOMl 

moo 


T(I8nt^,NEST)=C78 

11200 

102 

WE51'«1C0L ^ . 

11300 


T(IKOw,NEST)s-C78-^33*£IROW-U-C44*(IROW-1)*(IROW-2)-THETV|' 

11400 


1-.C22 . ■ 

11500 


MESlsiCOLPl 

11600 


TFCMKST.GT.K) GO TO 105 

11700 


T£IftOW,NE5T)sC33»lCOL 

11600 


MES’r«IC0Ijp2 

liyoo 


XFCftiEST.GT.N) GO TO 105 

12000 


T{IB0i(^,NEST) = ICQL*IC00Pl*C44 

121O0 


GO TO 105 

12200 

43 3 

DO 434 Klal^N 

12300 


IKsKl 

12400 


T(IRnw,IK)*1.0 

12500 

43 4 

CONTIMUE 

12600 

1 05 

COMTI^^UE 

12700 

C 


12800 


BRACE5aTHETR*ANUMBP 

12900 


DO 439 IR0Ws1,N 

13000 


irCIROW,EQ,N) GO TO 438 

13100 


irCIR0W,EQ,2) GO TO 437 

13200 


BCIROW,l3a-FFDClROW-13/BRACE5 

13100 


GO TO 439 

13400 

437 

BCIR0«#i)s-FFDCIROW-l)/BPACE5-C79/AKFJMBH 

tssoo 


GO TO 439 

13600 

438 

RClHOW,l)s+l,0 

13700 

439 

CCJNTIt'IUE 

13B00 


RETURN 

13900 


END 

14000 

C 


14100 

c 


14200 


SUBKHUTINE LIN(NUME0,FFD,ANUMBR,THETR,C22,C33,C78,C79,C44,NXf 

14300 


tDA6H12,NrEi^P3,ri,IPRIN3 

14400 


1450m 

■ c 


14600 

n 

THIS SUBROUTINE SOLVES THE EQUATIONS OF SECOND SUBGROUP 

14700 

c 

nr THE FIRST GROUP OF EQUATIONS , INITIALLY 'NUMEQ' EQUATIONS ARE 

14800 

c 

SOLVED BY CALLING SETl lO OBTAIN COEFFICIENT MATRIX AND MATIN 

14900 

c 

TO THIS EQllAlUNS.THE ANSWER IS OBTAINED AS BCl.l) 1 = 1 , — -Nimfi:Q, 

15000 

c 

THE LAST rWU VALUES OF THE SET ABE TESTED TO BE LESS THAN 

15100 

c 

1,0t>07 AND THE PROCESS IS REPEATED BY INCREASING THE NUMBER OF 

15200 

c 

EQUATIONS BY 2, FOR DETAILS SEE CHAPTER 4. 

15300 

c 


15400 


15500 


double precision T£25,25),B(25, l),FFDCO:20) 

15600 


DOUBLE PRECISION ANUMBR, THEIR, C22, C33 ,C78 ,C79,C44 ,DASHI2 

15700 


nOfIBLE PRECISION SUM3rSUM4 

15800 


NTEMP'siMUMEO 

15900 

.334 

GAT.L SETl (NTEMP,FFD,ANUM8R,THETR,C22,C33,C7 8,C79,C44,T,B) 



160<'0 


Hal 

16100 


CAIO i^AftilCT,NTEMP,B,M,DETEH) 

162UO 


MTEinPlaNTEMp-1 

16300 ■ 


M'X'EMP2a/^Tfi:MP-2 

1640U 


irCbCe^TEMPl ,1) ,GE,i,0E-07) GO TO 338 

16600 


TFC8CNTEMP2,U,GE.l,0E"07) GO TO 338 

16bOU 


MTE«P3=f^TEMP-3 

16700 


TFCIPBIN.EQ.O) GO TO 339 , , 

16&00 


WRXTEtnX,lU9) CNTEMP3,(RCIK,l>,IKal,NTEHP3j) 

16900 

1119 

FORMATC/SX, »ANSi4E8 FROM LiN V/5X , *NXEMPi= ' , 12// , 2 (5X ,F15 , 8 1) 

I700u 


GO TO 339 

17100 

338 

^XKMPaMTE«Pt2 

17200 


GO TO 334 

17300 

339 

SUM3»0,0 

17400 


,SUM4»0,0 

17500 


DO 33S IK=1 ,NTEMP3 

176O0 


KI»IK 

17700 


SUM3=Sim3+B(KI,n^CKI-l) 

17800 


SUM4»SUM4+B(K1,1) 

17900 

335 

CQNTIllUE 

18000 


0ASHI2«SUM3/SUM4 

18100 


RfifORH 

182O0 


EUn ' ' '■ ■ 

18100 

C 


11400 

C ■ 


■■ilsoo 


SUBROUTIME MATINCA,R,B,W#0EfER) 

18600 


18700 

C 


188O0 

r 

THIS SUBROUTINE INVERTS THE COEFFICIENT MATRIX OF A SET 

18900 

c. 

OF OTNEAR EQUATIONS. A REPRESENTS THE COEFFICIENT MATRIX ,N IS 

19000 

C! 

THE MUMBEP OF EQUATIONS , PUT Mai IF ONL.I SOLUTION OF EQUATIONS IS 

19100 

c 

MS.rrKii.}HF SPLUIIUN IS REPRESENTED BY BCI , 1) ,1 = 1 — ,N , 

19200 

c 


10 30 0 


19400 


n.dlhMl PRECISION BC25,U 

19500 


nnupLE PRECISION AC25,25) 

1960u 


DOUBLE PRECISION AMAX, SWAP, PIVOT, DEfRM,T 

19/00 


Dli'^ENSTON lPIVT(4a),INDEXC40,2) 

19800 


FQnH/ALENCE(IROW,JRaW) ,(ICOLM,JCOLM) , (AMAX, T, SWAP) 

19900 

10 

DETLPal.O 

2000U 

15 

no 20 -lal ,N 

00 100 

20 

rprvTCJiao 

00200 

30 

nu bsu (=i,N 

003OO 

40 

AHAXaU ,0 

00400 

45 

nu 105 3 sl,N 

00500 

50 

ff- ( tPlVXC n-l j 60,105,60 

OOtjOO 

60 

!)U lOU 6=1 ,N 

007OO 

10 

H (IPIVT(X)-U 80,100,740 

00800 

HO 

TFf AMAX-A0S(A(.1,K))) 85,100,100 



ooyuu 

85 

iKOw»a 

oaoou 

90 

lCOtiM»K 

oiiuy 

loo 

CONTINUE , ■ . . 

01200 

105 

cuwnMUE 

OIJOO 

110 

iPivTtjcoii«jiBiPivrcicoiiH7+i 

0l400 

130 

irClROW-lCOOM) 140,260,140 

oisoo 

i4o 

OifgRp-PEXER 

oifeoo 

150 

OU 200 

oi*/oo 

IbO 

S#AP»A€IR0«,lil 

oiwoo 

ifo 

A t i OB ii ) » A C 1 COfcM , 5 J 

oiyoo 

200 

A(JtC00M,liJ«8WAP 

02000 

205 

260;250,210 

02100 

210 

O0’25O 5»1,M 

02200 

220 

5wAF»B(|R0if,5) 

02^00 

230 

8UK0w;5)»B(IC05M,5J 

02400 

250 

iCJCOOM,l,j»5WAP 

02500 

250 

lliOlXtI,lS»IROB 

02000 

270 

iNDI;xCl,2i«lcoiiR 

02/00 

310 

PlVOTiA < leOiiH , ICOOM ) 

02000 

320 

OiXRRfOiTgRfPIVOI 

02900 

310 

ACi4‘o5«;ieo5M)«i,o 

01000 

340 

00 'IfO ■ ■ 

OSlOO 

350 

AaOO0lt^'fcJ«A(lCO0R,O)/PJ¥OT 

0l200 

355 

I#«»J 38O;380,35O 

.''.'iOilOO 

35o 

0O“3to'5»liM 

ol4oo 

170 

8(ieO0M,07«BClCO0M,0)/PIVOT 

03500 

380 

00 550 lil«l,N 

04b00 

3yo 

irCtl-lCOOMJ 400,550,400 

oS/oo 

400 

T«A(01,1CO0M) 

OJBOO 

420 

Al0i,iCO0M)*O,O 

03900 

430 

DO 450 0P1,» 

04000 

450 

ACi:*l,L3iACH,0)-A(ICO0M,0)4f 

04100 

455 

IFCM) 550,550,460 

04200 

450 

DO 600 Oal,M 

04300 

500 

BCl4l,l.)aB(:i*l,L)»BCICO0M,0)*f 

04400 

550 

CONTIBOK 

04500 

600 

DO 710 1*1, N 

04600 

610 


04700 

620 

tF(INDEXCl.,l3*INDEXCt,2)) 630,710,630 

04800 

630 

JR0W*INDEXCD,1) 

04900 

640 

aC00Mair«DEX(l/,2) 

05000 

650 

DO 705 Kal,N 

05100 

660 

SWAP«ACK,aROW) 

05200 

670 

ACK,vTROB)*ACK,aCOl.«) 

05300 

700 

ACK,JCOLM)*SBAP 

05400 

705 

CONTINUE 

05SOO 

710 

CONTINUE 

05500 


DO U K*l,N 

05700 


irCIPIVTCK)-!) 12,11,12 

05800 

1 1 

CONTINUE 



10900 


SUMlnO.O 

11000 


DO 467 Kll«0, MOMENT 

11100 


' KIWKII 

11200 


IMONT«MQMENT«KI 

11300 


SUMIb5UM1+FACT0(M0MENT)»SSMAI^CKI,IRAD)»GSMAI.CIM0NT,IRAD)/ 

11400 


lCFACTOCIMONT)*rACTOCKI)? 

11500 

467 

continue 

1160O 


TERMCI84D3*»C4034SUM1*XRAD(I8AD3/V 

11700 

468 

continue ' . ' ; 

11800 


DO 9 IR0W*l,NTEHP3 

11900 


IFCIR0W,GT,1) GO TO 7 

12000 


BCIROW,U»'»TERMC23»C4044XI,IVECMOMENT,2)/V 

00100 


IwFDDEDtMOMENT, 03*THETVE 

00200 


GO TO 9 

00300 

7 

BCIRON,13»-TERMCIROWtl)-C404*lROW4 

00400 


1XLIVE(M0MENT,IR0W+13/V'*C63»IRP4XUIVE(M0MENT,IR0W-13 

00500 


1»C102*PHI»XI4IVE(MOMENT,IROW»*1 3-rDDED (MOMENT, IROW-1 3 ♦THETVE 

00600 

9 

continue 

00700 


RETURN 

00800 


END , 

00900 

C 


OlOOO 

c 


otioo 


SUBROUTINE SET2CNTEMP3,C22,C33,C44,C78,C79,C63,C102,C101 ,fRP, 

01200 


i«0MENT,PHI,THEtVE,FDliIV,XRAD,XUlVE,T,i3 

0t300 


01400 

C 


01500 

c 

THIS SUBROUTINE CAI<CULATES THE COEFFICIENT MATRIX OF THE 

01600 

c 

SET OF SET LINEAR EQUATIONS OF GROUP 2.(I,E, LIVE POLYMER CHAIN 

017 00 

c 

LENGTH DISTRIBUTI0N3 COEFFICIENT MATRIX REPRESENTED T, VECTOR OF 

01800 

c 

CONSTANTS IS BCI,13 > 

01900 

c 


02000 


02100 


DOUBLE PRECISION T(25, 253,6(25,1), XRA0COI2O3 

02200 


DOUBLE PRECISION FDLIVC5,203 ,XLIVE(0s5,203 ,fERMC20) 

02300 


DOUBLE PRECISION C22,C33,C44,C78,C63,C101 ,C102,PHI,TRP,THETVB 

02400 


DOUBLE PRECISION SUMl,C79 

02500 


NTEMP2aNTEMP3»l 

02600 


DO 4 IROWal ,NTEMP2 

02700 


DO 4 iCOLsl ,NTEMP2 

02800 


TCIROW,ICf3L3a0,0 

02900 

4 

CONTINUE 

03000 

c 


03100 


DO 5 IROWal, NTEMP2 

03200 


ICOLalROW 

03300 


ICOLM2*ICOL-2 

03400 


ICOLMlaiCOL-l ■ . , 

035OO 


ICOLPlalCQL+1 

03600 


IC0LP2»IC0L+2 

03700 

c 




03«yf) 


NESTwlCOIiMa ; 

03900 


ir(NE5T,LE,-13 GO TO 2 

04000 


ir(»iEST,LE,03 00 TO 3 

04100 


TCI«0W,NEST)*C22 

04200 

3 

NESTbICOLMI 

04300 


■ TC1R0W,NE5T3»C78 

04400 

2 

nest»icoi:» 

04500 


T(IROW,NEST)*»C78»C22-C33»ICOIi*C63*T8P»Ci02»PHI-C44*ICOL*ICOLMl 

04600 


1*THETVE 

04700 


ffEST*lCOliPl ; ! 

04800 


IFCNEST,GT,NTEMP2) GO TO 5 

04900 


T(IR0WrNEST)aC334IC0t. 

05000 


NEST»lC0IiP2 

05100 


IFCNteST,GT.NTEMP2) GO TO 5 

05200 


TC1R0W,NEST)» C44*IC0LPl*IC0I. 

05300 

S 

CONTINUE 

0S400 


DO 468 IRAD1»1,NTEMP2 

05500 


IRA0*1RAD1 

05600 


SUM1«0»0 

05700 


DO 467 KIl»l, MOMENT 

05800 


KlaKIl 

05900 


IMONT»MO«ENT»KI 

06000 


Simi«SUMl4rACTOCMOMENT)»XDIVEClMONf,IRA0)/CrACTO(KI34 

06100 


IFACfOdMONT 3) 

06200 

467 

CONTINUE 

06300 


TERMCIRAD3*Ct014PHI»6UMl 

06400 

468 

CONTINUE 

06500 

C 


06600 


DO 9 IR0W»1,NTEMP2 

06700 


IFCIRQW,GT,l) GO TO 7 

06800 


BCIROW^l) a-C78»XRADCQ3"Cl024PHI41,04XRAD(l)*TMETVE» 

06900 


IFDLIV (MOMENT, 13 "TERM Cl3"C634TRP4i,04XRADCi)»C79 

07000 


GO TO 9 

07100 

7 

8ClROW,l)a«C78*XRADCIROW"l3-2,O»C224XRA0(IROW-23"ClO2»PHl» 

07200 


1 IROW4XRA0 (IROW 3 -THETVE»FDLI VC moment, IRON) "TERMC IRON! 

07300 


1-C63*TRP»IROW4XRADCIROW3 

07400 

9 

CONTINUE 

07500 


RETURN 

07600 


END 

07700 

C 


07800 

C 


07900 


FUNCTION DIFF2(AMV,V) 

08000 


08100 

C 


08200 

c 

THE INPUT FDR THIS FUNCTION CONSISTS OF AMV AND V,AMVCI) 

08300 

c 

ARE THE MOMENTS OF A DISTRIBUTION, CI«0,1 ,2— "•) THE FUNCTION 

08400 

c 

RETURNS THE VAUUE OF DISTRIBUTION FUNCTION EVALUATED AT V, 

08500 

c 

FOR DETAILS SEE CHAPTER 4, 

08600 

c 


08700 




08800 

08900 

09000 

09100 

09200 

09300 

09400 

09500 

09600 

09700 

09800 

09900 

toooo 

10100 

10200 

10300 

10400 

10500 

10600 

10700 

10800 

10900 

11000 

moo 

11200 

11300 

11400 

11500 

11600 

11700 

11800 

11900 

12000 

12100 

12200 

12300 

12400 

12500 

12600 

12700 

12800 

12900 

13000 

13100 

13200 

13300 

13400 

13500 

13600 

13700 


DOUBLE PRECISION AMVCOjS) ,GXC6) ,CALPHA(6) ,Al,PHAL(65 
C AMV ARE PIRST FOUR MOMENTS OF DISTRIBUTION 

C AMVZER IS THE ZEROTH ABSOLUTE MOMENT OF DISTRIBUTION 

C X IS THE X CO-ORDINATE OF WHICH 

C THE FUNCTION IS TO BE EVALUATED. 

C RESULT IS THE FUNCTION VALUE CALCULATED AT 

C AFORESAID X CO-ORDINATE MVi, 

C ■ 

C INPUT PORTION— — 

C CONVERT MV INTO GXCI.E, MOMENTS OF NORMALISED 

C DISTRIBUTION.) 

C 

DOUBLE PRECISION AMVZER, PDI .ALPHA, ALPHAl .ALPITI . ALPIKl .XMINU . SUM 

DOUBLE PRECISION V 

AMVZER^AMVCO) 

PDI«AMV(2)*AMVZER/CAMVCU*AMVC1)) 

GXC1)*1,0/(PDI-1.0) 

GX(2)s»PDI»GX£l)*GXCl) 

ALPHA*GXC1)-1,0 

ALPHA1«GXC1? 

NMAX*4 

DO 1S8 KI«3,NMAX 
IKi»Kl" 

GX£lK)aCCALPHAlMAMVZER/AMV£ 13 )4*IK)*AMV(IK)/ AMVZER 
118 CONTINUE 

C CONVERT ABSOLUTE CO-ORDINATE V TO THE C0RRES5P0NDING 

C NORMALISED CO-ORDINATE X, 

X*£ALPHAl»AMVZER/AMV(l))*V 
C CALCULATE CALPHA, 

PO 100 K1»3,NMAX 
IKbKI 

CALPHA(IK)3l,0 
DO 100 ITAM*1,IK 
. ITEMPbITAM 

IKITAMbIK-ITEMP 
ALPITIbRIiPH ai + itemp 

CALPHAClK)sCALPHACIK) + £ C-1)»*ITEMP)*FACT0CIK)*GAMM£ALPHA1)» 

IGX C ITEMP) / CFACTOC IKITAM)*FACTO(ITEMP)*GAMM(ALPITl ) ) 

100 CONTINUE 

C CALCULATION OF MODIFIED LAGURE8 POLVNOMIAL 

DO 200 KIa3,NMAX 
IKbkI 

ALPIKIbALPHAI+IK 

ALPHALC1K3bgammcaLPIK1)/CGAMMCALPHA1)»PACTOC1K33 

DU 200 ITAM*l,IK 

ITEMPbITAM 

IKITAM«IK-ITEMP 

AltPTTl bALPHAI + ITEMP 

ALPHALCIK3BALPHAL£IK)+CC-13f»ITEMP)4GAMM(ALPIK134 

lCX*»ITF.Mp)/CFACTO(IKITAM)fGAMMCALPITl)»rACTO(ITEMP)) 



13800 

200 


CONTISfOE 

13900 

C 


CAOCUtiATlON or NORMAliISED DISTRIBUTION PUNCTON, 

14000 



SUM«1 ,0 

14100 



DO 300 KI*3,NMAX 

14200 



IKsKI 

14300 



SUMj#SUM+Cai4pHA(IK)»Al,PHAI,ClK3 

14400 

300 


CONTINUE 

14500 



XWINU»»X 

14500 



GN0RMa9SUM»EXPCXMlWU)»CX»4A6PHA)/GAM«(Ai:»PHAl) 

14700 

c 


CAUCUUATION Of VALUE OP TJE ABSOLUTE FUNX, 

14H00 



DIFF2«*GN0RMfALPHAl4AMVZBR4AMVZER/AMVCl) 

14900 



RETURN 

15000 



END 

15100 

c 



15200 

c 



15300 

c 



15400 

c 



15500 



SUBROUTINE CENTR (N,M,K,IEVi , I , XC ,X .Kl ) 

15600 


15700 

c 



15800 

c 


THIS SUBROUTINE IS REaUlRED FOR OPTIMIZATION, IT 

ISfOO 

C 


CALCULATES THE COEFFICIENTS OF THE CENTROID OF THE COMPLEX, 

16000 

c 



16100 

^ l|e f f « 9K « f f f f f; f y f f f f it{ f f f f f f f f f f f II y ;|C f y « f f f f f i|t f f f )|C f f f « « V « f f ;|c « « 9); « i|t « f )|C IK ]|C 

16200 



DOUBLE PRECISION XCH,M),XCCN3 

16300 



DO 20 J»l,N 

16400 



XCtJJ«0,0 

16500 



DO 10 IL«l,Kl 

16600 



XCCa)aXCCJ)+X(IL,J) 

16700 


10 

CONTINUE 

16800 



RKaKl 

16900 



xccu)®(xcca)-x(iEvi,a)>/(RK*i,0) 

17000 


20 

CONTINUE / 

17100 



RETURN 

17200 



END 

17300 

c 



17400 

c 



17500 



SUBROUTINE CHECK CN, M,K,X,G,H, I, KODE,XC, DELTA, Kl) 

17600 

(;;iK;K¥«K’K)Kl!K4»t« iK««3KiK«!«)K444*f *f f 4«4«)Kf iK44>K44)K««f 

17700 

C 



17800 

c 


THIS SUBROUTINE IS REQUIRED FOR OPTIMIZATION, IT CHECKS 

17900 

c 


IF A GIVE9 POINT SATISFIES THE GIVEN SET OF CONSTRAINTS AND 

18000 

c 


AND CALCULATES A FEASIBLE VERTEX FOR THE COMPLEX, 

18100 

c 



18200 

(;;iKK!KV)K»f )KK)K«!4)K44)l‘f 4 f 44f«4t4!4f f 4**4* 

18300 



DOUBLE PRECISION XCK,M) ,GCM3 ,H(M) ,XCCN) 

18400 


10 

KTbO 

18500 



CALL C0NSTCN,M,K,X,G,H,I) 

l«600 



DO 50 Jal,N 

18700 



IFCXCI,J)-GtU)3 20,20,30 



18800 


20 

XCI,a)»GCJHDEI<fA 

18900 



60 TO 50 , , ^ : 

19000 


30 

irCH(J)-XCI,J)5 40,40,50 

19100 


40 

XCI,J)wHCjJ-DEI*TA 

19200 


50 

CONTINUE ^ ; 

19300 



N0»5 . ' ■ ■ ■ 

19400 

C 


WRITE CNo,43 tcj,xci,a),xcca),x( 2 ,a)) ,j=i,?#j 

19500 



irCKODE) 110,110,60 i 

19600 


60 

NMwM+i ' , ; 

19700 



DO 100 i 

19800 



CAIiI^ G09STCN,M,K,X,G,H,I5 

19900 



iF(X(i,a)-G(a) j 80,70,70 ; 

20000 


70 

ircHCJ)-x(i,a)> 80 , 100,100 

20100 


80 

IC?t«I ■ 

20200 



KT*1 

20300 



CALt CENTR(N,M,K,IEyi,I,XC,X,Ki3 i 

20400 



DO 90 \ 

20500 



xci,Ja)»(xci,ja)txccaa))/ 2 ,o 

20600 


90 

CONTINUE ' ■ ' ■ 

20700 


100 

CONTINUE . , 

20800 



ircKf) iio,uo/,io 

20900 


110 

RETURN 

21000 



END 

2H0O 

c 



21200 

c 



21300 



SUiROUTINE C0NSXCN,M,K,1TMAX,AUPHA, BETA, GAMMA, DEUTA,X,R,r, 

21400 



1IT,IEV2,N0,G,H,XC,IPRIN,DASHI,0NETH,TW0TH,BRACEI,BRACE2, 

21500 



1 BRACES, C50,C51,C52,G53,C54,C55,CS56,C56,CT1,C7,THETR,DUMMT,XFD, 

21600 



iC65,SAT) 

21700 

c; ¥ f « 3^ y 3|C ]«c # f )); :|c f Ip ](( f f f f f f f f f f i|t f f f f « 9 :|i !|E f f; ifc « f f )|( J|S « « y f f f f !|e I|c « « f f ft « « f f 

21800 

C 



21900 

C 


THIS SUBROUTINE IS USED TO FIND MAXIMA/MINIMA OF A 

22000 

c 


CONSTRAINED MULTIVARIABUE FUNCTION BY BOX^S METHOD, IT IS USED 

22100 

c 


TO FIND INITIAL GUESSES TO SOLVE THE EQUATIONS OP FIRST SUBGROUP 

22200 

c 


or GROUP UFOR DETAILS SEE CHAPTER 4 AND REF 32 

22300 

c 



22400 

c 


DISCRIPTION OF PRINCIPAL VARIABLES,.*, 

22500 

c 


ALPHA„„,„.,,.,,,,„,„.,»,,,,CQMPLEX EXPANSION PARAMETER, 

22600 

c 


CUSUALLY 1,3) 

22700 

r 


Bill, CRltl^SHXOl^ 

22800 

c 



22900 

c 


DLI,^T TO IFXND ^ f*l!A.SXBIiE FOXHT 

23000 

c 


mm THE CftWUtATED VERTEX Of TH 

23100 

c 


COHPtEX IS FOUND TO BE OUTSIDE 

23200 

c 


mj CONSTRAINfi ^ 

23300 

c 


M PWMIMWmilWMmfiMlIplillW QF VARIABLES || 

23400 

c 


H OF COMSTRAWTSf 

23500 

c 


iC OF VERTICES# 

23600 

c 


I TM HUMBER OF ITERATIONS 

23700 

c 


hhhmm fm 



2J900 

24000 

24100 

24200 

24300 

24400 

24500 

24600 

24700 

24800 

24900 

2S000 

25100 

25200 

25300 

25400 

25500 

25600 

25700 

25800 

25900 

26000 

26100 

26200 

26300 

26400 

26500 

26600 

26700 

26800 

26900 

27000 

27100 

27200 

27300 

27400 

27500 

27600 

27700 

27800 

27900 

28000 

28100 

28200 

28300 

28400 

28500 

28600 

28700 



10 

20 


30 

40 


SO 


51 

52 
18 


159 

55 

56 
65 


70 


K0DB»0 

irCM-N)20,20,10 

K0DE»1 
CONTINUE 
DO 40 II«2,K 
00 30 a»l,M 
xciifjmo^o 
CONTINUE 
CONTINUE 
00 65 1H*2,K 
00 50 

Ckhh CONSTCN,H,K,X,S,H,I) 
xcii, jjbgcj) 48(11,0) 

CONTINUE 

Klsll 

DELTA, Kl) 

IFClI-2351,51,55 # / 

irCIPRIN)52,65,52 

WRITE(N0,018) 

F0RMATCIH1//2X, ^COORDINATES OF INITIAL COMPLEX') 
WRITE (NO, 159) CIO, J,X(lO,a) ,0»1 ,N) 

F0RMATC/3(2X,2HXC,I2,1H,,I2,4H) a ,E13,6)) 

T O' f T t>OT HI > iSt gg r rr , » , A ^ , w / ^ 


p UKWATt/4(2X,2H 
IFCIPRIN) 56,65,30 
WRITE (NO, 159) cii,j,xcii,a),a 
CONTINUE 
KlaK 

on 70 ILai,K 
I ® 1 L 

CALL rUNCCN, 

1BRACE3,''»'« 

1C65,SAI 
CONTINUE 
K0IINT«1 
I AaO 

IFCIPRIN) 


■«lfN) 


NCCN,M,K,X,F,I,DASHI ,BRACE1 ,BRACE2, 

,C50,C5l ,C52,C53,C54, CSS, C556,C56,CT1 ,C7,THETR, DUMMY, XFD, 

E 





20800 

72 

WRITBCN0,2n 

28900 

21 

F0RMATC/2X, 'VALUES OF THE FUNCTION') i 

29000 


hritecnd, 22) ca,Fca),jai,K) : 

29100 

22 

F0RMATC/,3C2X,2HFC,12,4H} s ,E13.6)5 

29200 

80 

IEV1*1 

29300 


DO 100 ICMa2,K 

29400 


IPCF(IEVl)-r(ICH)) 100,100,90 

29500 

90 

IEV1»ICM ' : 

29600 

100 

CONTINUE 

29700 


IEV2BI i 

29800 


DO 120 ICMa2,K i 

29900 


ircrciEvai-^rcicw)) 110 , 110,120 

30000 

110 

IEV2WICM 

30100 

120 

CONTINUE 

30200 


ircrCIEV23«(F(lEVl)+BETA)) 140,130,130 i 

30300 

130 

KOUNTbI I 

30400 


GO TO 150 ; 

30500 

140 

KOUNT«FOUNT+1 i 

30600 


IF(K0UNT-GAMMA3 150,240,240 j 

30700 

ISO 

CALL CEMfR(N,f<l,K,XEVl,I,XC,X,Ki) ( 

30800 

' " i ' ■ ■ 

DO 160 < 1 

30900 


XCIEVl ,*3i3)wC1.0+ALPHA)4XC(Ja)«ALPHA#XClEVl,iJJ) i 

31000 

160 

CONTINUE i 

31 too 


IbIEVI ■ i 

31200 

' 

CALL CHECKCN,M,K,X,G,H,I»KODE,XC, DELTA, Kl) 

31300 


CALL FUNC C N , M , K , X , F , I , DASHI , BRACE 1 , BRACE2 , 

31400 


1BRACE3,C50,C51,C52,C53,C54,C55,C5S6,C56,CT1,C7,THETR, DUMMY, XFD, | 

31500 


1C65,SAI) i 

31600 

170 

IEV2*1 ■ i 

31700 


DO 190 ICMa2,K 

31B00 


IF(FCIEV2)-F(1CM)) 190,190,180 

31900 

180 

IEV2=ICM 

32000 

too 

continue 

32100 


TP(IEV2-1EV1) 220,200,220 

32200 

200 

DO 210 JJal,N 

32300 


x(ievi,aj)a(xciEvi,ao3+xccuj))/2,o 

32400 

210 

CONTINUE 

32500 


lalEVl 

32600 


CALL CHECK(N,M,K,X,G,H, I rKOOE,XC, DELTA, Kl) 

32700 


CALL FUNC (N,M,K,X, F, I, DASHI, BRACEI,BRACE2, 

32800 


lrtRACE3,C50,C5l,C52,C53,C54,C55,C556,C56,CTl,C7,THETR,DUMMY,XFD, 

32900 


1C65,v5AI3 

33000 


GO TO 170 

33100 

220 

CONTINUE 

33200 


IFCIT/50*50,NE.IT) IPRIN«0 

33300 


mtype«i 

33400 


irCIT/50»50,NK,IT) MTYPE«0 

33500 


IFCMTYPE.EQ.O) go to 4511 

33600 

C 

IFCTT.GT, 300) IPRIN*1 

33700 


TYPE *,IT 



3Jeoo 
33900 
34000 
34100 
34200 
34300 
34400 
34500 
34500 
34700 
34800 
34900 
35000 
35100 
35200 
35300 
35400 
35500 
35600 
35700 
35800 
35900 
36000 
IftOO 
.'■W 3100 
36300 
36400 
36500 
36600 
36700 
36800 
36900 
37000 
37100 
37200 
37300 
37400 
37500 
37600 
37700 
37800 
37900 
38000 
38100 
38200 
38 300 
38400 
38500 
38600 
38700 


4511 

230 

23 

24 


25 

26 
228 


IFCIP8IN) 230,228,230 
W8IfC(NO,233lT 

NOMSERMS) 

WRITECN0,24) 

corrected point- 3 

WRlTECNo'^fr 

WRITE(N0;22)CI,F(I),I*1,K} 

WR 1 TECN 0,253 » m / 

CEMROID') 

JRITECNO,26)CJC,XCCaC),JCal,N3 

IPRIN«IPRIN1 
irCIT-lTMAX) 80,80,240 
RETURN 
END 


240 

C 
C 

SUBROUTINE CONST(N,M,K,X, a, H, I) 

PRECISION X ( 1 4 f 16) * GC i 6 ) f H C 1 63 

c COSIRAiJs” “"I"''® VARIOUS EXPLICIT AND IMPLICIT 

r constraints, 

r Tir"*"::’ -—-.——-.-opper constraints, 

t I.N+1, — ....implicit constraints' 

HC1)=!4,9 

GC2)=0,0 

H(2)3!1,o 

GC3)b0,0 

HC3)»l,g2 


GC4)=1 ,OE-06 
HC4)*1 ,0E+Q6 
GC5)»1 ,0E-06 
H(53sl,OE'f06 
GC6)s,0l 
HC6)*l ,0 
G(7)»1,0E-.06 
HC73al,0E+03 
GC8)»1,0E-06 
HCai«l,0E+04 
GC9)«1 ,0E«0h 
HC9)»1 ,OE+05 



38HU0 

38900 

39000 

39100 

39200 

39300 

39400 

39500 

39500 

39700 

39800 

39900 

40000 

40100 

40200 

40300 

404U0 

40500 

40600 

40700 

40800 

40900 

41000 

41100 

41200 

41300 

41400 

41500 

41600 

41 700 

41 800 

41900 

42000 

42100 

42200 

42300 

42400 

42500 

42600 

42700 

42800 


GC10)»1,OE*06 
H(10)al ,0E+06 
G(ll )*1 ,0E««07 
HCU)*1,0E+07 
IFCN-M) 1,2,1 

XCI,12)aXCI,7)*X(I,9)/XCI,83»»2 

GC123^1^0 

Htl2)*sl0,0 

XCI,13)«Xfl,8)/XCI,7} 

GC13)«0,5 

HC13)»10,0 

XCI,14)*XC1,9)/XCI,8) 

G(14)«1,0 

HC14)»10 

Xa,15)»XCI,10>/XCI,9) 

0(153*1,0 ■ 

HClS)al0,0 

XCI,16)*XCI,113/XCI,10) 

GC163*1.0 

H(163*10,0 

RSfURN ' ■■ T 

ENB . ^ ' 


,C52,C53,C54,C55,C556.C56.CT1 C7 THPTR nriiiMV ¥crk% 

THIS SUSR'OUTINE RETURfIS THE VALUE OE THE OBilKM*T¥!*’ 

? 0Ba.ciivS'‘!!L?T,r,°'’M'Dl'rR0« 

^ R,H,S, OF EQUATIONS OF FIRST SUBOROUp OF GROUP 1, 

nnURr'fc' fF(K3 ,XFD(1 1 3 ,Z(6) ,PFDC6) ,FQI#DC6) 

DOlIRLh PRECISION XC253 ,RU«y (6) 

“SIImf PRFrJ^Jn» ™f5'!''°r'®S»':EJ,BRACE3,BR«CM,TBETR,DUKMI 

njUBljR PRECISION G50,C51 ,C52,C53,C54»C55.C5SS .C 56 07 OT* 

OOUbS PRFCT5TnN plf-oio^r® ' 

OOUBLE PRECISION FOURTH, FIFTH, SIX, SEVENT, SUM 
w ■'# A X s 4 


42900 C 


43000 

ONETHal,0/3, 

431 00 

TwnTH«2,0/3, 

43200 C 


43300 

*^l*8,08E-04 

43400 

W2*47,69 

43500 

WJa3,807E..02 

43600 

W4«2,937E+06 

43700 

W5a 7,O998E-0l 



43B00 


W6»0,63 ■ ■ ■ : 

43900 


W7*0,1981 ' ■ 

44000 


F 0 LDCn«e 8 *SUlE +02 

441.00 


FDIjD(2)*7,3B+06 

44200 


r0LD(3)tt9,432E+10 

44300 


FDIjD(4)al,622E+lS 

44400 

C 


44600 


l5UMMlfl»l,0E4-06 ■ 

44600 


DUMMy2=»l,0E-02 

44700 


DUWMTr3»1.0E-03 

44800 


D0MMy4**4,2985E+03 

44900 

c 

■ ■■ ■■■ '.V 

4S000 

c 

X ARE THE OUTPUT VARIAB1,ES FROM THE REACTOR 

45100 

c 

XFD ARE THE FEED CHARECTARISTIC OF AFORESAID VARIABLES 

45200 


DO 450 IK*1,INN 

45300 


K1«IK 

45400 


X(KI3»XDUMTCa,KI) ■ 

45500 

450 

COHTINUE 

45600 

C 


45700 


H«K»7 

45800 


rZERO»XC73 

45900 


rfDEER»XFD(73 

46000 


DO 210 I»TEMP»l,HMAX 

46S00 


«TE«»INTEMP 

46200 


ZCNTEM3aiXCMTEM+7J 

46300 


FrDCNTEM)*XFDCNTEM+7) 

46400 

210 

COHTIHUE 

46500 

1 

RES«CXrDCl)-XCl3)/THEfR 

46600 


FIRSTfiRES/Wl 

46700 

2 

RESaC56*CXFDC23-XC2))/THETR-XC23 

46800 


SECnND=RES/VI2 

46«00 

3 

BESaCXFD(33-XC3))/CTHETR»DUM«y2}-C53*Xe6)*CFZERO/DUMMY33*DA5HI 

47000 


THIRDaRES/W3 

47100 

4 

RESaC504X(l)-X(4J-C51»C0NPR(TW0fH,Z,FZER0J*CDUMMY4*»TW0TH)/ 

47200 


lDUMMY3»C525XC53/DUMMyi 

47 300 


FnURTHsRES/W4 

47400 

5 

RESaC54»X(3)/DUMMy2«CXC53/DUMM¥U«CS5*(ZCl)*DUMMY4/DUMMY3)* 

47500 


tXf6)-C5561'X(4) 

47600 


FIFTHaRES/HS 

47700 

6 

RF:SaXC6)*0,615 

47800 


RIX3RES/W6 

47900 


RKSaC7»XC2)»XC43/tXt43tCCOMFR(TWQTH,Z,rZERO)4(DUMMY4^4TWOTH3/ 

40000 


IDUMMYIDf CFFDZER-rZER03/CTHETR#DUMMY3) 

48100 

C 


48200 

^ f ;f! f f f ^ y s(! f f « « f f 4; f f f f 4 f f f )|! f f f f f f :(e 4 f 4: 4 f 4 f 4 4 4 4 t « f f V « f * f f * 

40300 

c 


48400 

c 

INSERT THIS LINE INSTEAD OF PREVIOUS LINE FOR SECOND APPROACH 

48500 

C6 

RESal,0«.XC6)+DLnGCXC6))t5AI*(l,0«X(63)+C654XC7)/ 

48600 

c 

UCilNFRtONETH,Z,FZERO)4(DUMHY4MONETH)3 

48700 

c 




48B00 

4B90Q 

490O0 

491U0 

49200 

49300 

49400 

49500 

49600 

49700 

49H00 

49900 

50000 

50100 

50200 

50300 

50400 

50500 

50600 

50700 

50000 

50900 

51000 

51100 

51200 

51300 

51400 

51500 

51600 

51700 

51800 

SI 900 

52000 

52100 

S2200 

52300 

52400 

52500 

52600 

52700 

52000 

52900 

53000 

531 00 

53200 

53300 

53400 

53500 

53600 

S3700 


C«* 


SEVBNI»RES/W7 
DO 650 9*1,NMAX 




IFCN,GT,1) GO TO 207 

RUMyCN3«RES/rOI.DCN5 C f«ETR»0UMMf 33 

GO TO 650 

CONTINUE 
SUMwO.O 

DO 550 IK*1,NMAX 
kh»ik . . 

SUM»SUMfABS(RUMT(KI3) 

CONTINUE ■ ' 

peu3»*Fca3 

RETURN ' 

END"’' 

SOSROUTINE 0G10CXI:4,XU,Y,ARV3 

C! Iff t.« ■»*' « , :#W f « *-!» Wlfc B . 


207 


650 


550 


FUNGTIfra^PvaT RETURNS THE VALUE OF INTEGRAL OF A 

Fu XriOM EVALUATED IM THE RANGE XL TO XU. IT IS nSFO f^n» 

distribution!”® CUMULATIVE DISTRIBUTION FROM A DIFFERENTIAL 


C 

c 
c 
c 

DOUBLE PRECISION AMV(0r53 
DOUBLE PRECISION A,B,C 

aso^b^^cxu-xd+xl 

BaCXU-Xfa) 

C * 0 ^ 4 8 6 9 5 3 3 0. 



53800 
53000 
54000 
54100 
54200 
54300 
54400 
54500 
54600 
54700 
54800 
54900 
55000 
55100 
55200 
55300 
55400 
55500 
55500 
55700 
55800 
55900 
56000 
56100 
56200 
56300 
56400 
56500 
56600 
56700 
568 00 
56900 
57000 
57100 
57200 
57 300 
57400 
575U0 
57600 
57700 
57800 
57900 
58000 
58100 
' 58200 
58300 
58400 
, 58500 

I 58600 
i 58700 


Y«B*CY+O,H776211'CDIFF2CAMV,Ai.C)+0IFr2(A«V,A-C))) 

RETURN ■ ; 

END . ^ 

c ' ■ ; 

c : 

SUBROUTINE DISTRCXIil ,PARTt , AMV, IWHICH,XCORD,yCORD,ZCORD ,XCU«) 

C ^ : 

C THIS SUBROUTINE SENERATES THE DISTRIBUTION FUNCTION FROM 

C ITS MOMENT. INITIAUW THE FUNCTION IS EVALUATED AT STEPS OF 

C PARTI STARTING WITH THE LOWER LIMIT OF XLl.THE USEFUL RANGE OF 

C DISTRIBUTION IS THEN CALCULATED .THE USEFUL RANGE IS THEN 

C DIVIDED INTO 50 INTERAVALS AND THE DISTRIBUTION FUNCTION IS ' 

C CALCULATED AT THESE 51 POINTS. XCORDCDIS THE SET OF X 

C CO-ORDINATES AT WHICH THE FUNCTION CALCULATED, YCORDCI) IS THE 

C SET OF FUNCTION VALUES AT THESE POINTS ON NUMBER BASIS, i 

C ZCORDtIJ IS THE SET OF FUNCTION VALUES ON WT, BASIS, XCMM I 

C REPRESENTS THE SET OF NORMALISED CUMULATIVE DISTRIBUTION 

C FUNCTION CALCULATED AT THE AFORESAID POINTS 

C IWHICH«0—— NORMALISED CUMULATIVE DISTRIBUTION IS NOT CALCULATED j 

C { 

c I 

DOUBLE PRECISION XCORDC0I50) ,TC0RD(0»50) #ZCORDCOl50) .XCUMfOsSO) 
DOUBLE PRECISION STARC50) .AMVCOlS) i 

DOUBLE PRECISION XL.XLt ,XU1 .XU.PARTl ,PART2,¥,YMAX 
DOUBLE PRECISION RATIO, AREA, TOTAL 1 

XL«XL1 ■ ■ ! 

C TO FIND USEFUL RANGE OF DISTRIBUTION, 

YMAXsO.O 

1 X=DIFF2(AMV,XL) 
irCY.GT.YMAXl GO TO 15 
RATIOaY/YMAX 

IFCRATID.GT.O.OOOUGO TO 16 
GO TO 2 

15 YMAXaY 

16 XLsXL+PARTl 
GO TO I 

C UPPER USEFUL LIMIT FOR THE DISTRIBUTION HAS BEEN OBTAINED 

2 XUlsaXL 
PART2a(XUl-XLl)/50,0 

XL«XL1 X 

C TO Firm DIFFERENTIAL AND CUMULATIVE DISTRIBUTION 

XCDROCO)*XLl 
YCOROCO)»DtFF2CAMV,XLl) 

ZCORDCO)aXCORD(0)*YCORDCO) 

XCUMCO)«0,0 
AREAan.O 
DO 3 KI*1,50 
XU*XL+PART2 



58800 

58900 

59000 

59100 

59200 

59300 

59400 

59500 

59600 

59700 

59800 

59900 

60000 

60100 

60200 

60300 

60400 

60500 

60600 

60700 

60800 

60900 

61000 

61100 

ilSOO 

61300 

61400 

61500 

61600 

61700 

61800 

61900 

62000 

62100 

62200 

62100 


6 

3 


4 

5 

C 

C 


XC0RDCKI)«XU 

YCO80fKI)»DIPr2(AMV,X[J) 

1 F ( I WHICH * EO 0 ) 00 yQ 0 
CALI, QG10CXL,XU,y,AMV3 
STARCKIlay 
AREAatAREA+STARCKI) 

XL»XU 

CONTINUE 

irCIWHlCM,EO,0) GO TO 5 

TPTAW?REr^ normalised CUMULATIVE. DISIRIBUIIOM 

area«o,o 

DO 4 KI«1 ,50 
AREA*AREA+STARCKI) 

xcumckdsarea/total 

CONTINUE 

RETURN 

END 


SpROUTINE FUNXXCX,RES,K,0ASHI, 

‘cm 0S5 • CS5« . 




r AMu>.,Mn UIS SUBROUTINE IS SIMILAR TO FUNX.FOR SECOND APPRQJ1<'H 

C CHANGE the NAME OF SUBROUTINE TO FUNX. SECOND APPROACH 

nZE"^ ’"^’.NFocuj.rcsi.FrMsnmDnr ’ 

DOUBLE PRECISION RE5,DASHl,0NETH,TW0TH 

DOUHLF PRrr?f?nw !Rj^®i!®’^*^®2'®^*^®^3'0UMM¥,THETR,SAI 

HMAXs4 


62400 C 

62500 C 

62600 C 

62700 

62800 

62900 

63000 

63100 

63200 

63300 

63400 

63S00 

63600 

63700 


Wl*s8 ,08E"04 

W2a47,69 

W3«3,807E-02 

W4*2,937E+06 

W5«7,0998E-01 

W6»0,63 

W7a0,1981 

FOltD(l)«8,SlllE+02 

FnLD(2)«7,3E+06 

FnLD(3)»9,4321+10 

FnLD(4)«l,622E+15 



63800 

c 


63900 


DUMMYl*l,0B+06 

64000 


DU«My2«»l,0E-02 

64100 


DUMMY3«U0E-03 

64200 


DyMMY4«4,2985E+03 

64300 

c 


64400 

c 

X ARE THE OUTPUT VARIABLES FROM THE REACTOR 

64500 

c 

XPD ARE THE PEED CHARECTARISTIC OP AFORESAID VARIABLES 

64600 


H«K*7 ■ ■ 

647Q0 


FZER0«XC6) : 

64800 


FPDZER«XFDC6) 

64900 


DO 210 INTEMPB!1,NMAX 

65000 


NTEM»IHTEMP 

65100 


FCNTEM)ftXCHTEMf6) 

65200 


rFDCSTEM)*XF0CNTE«+6) 

65300 

210 

CONTINUE ' 

65400 


irCK,OT*63 GO TO 205 

65500 


GO TO (1,2,3,4,5,6),K 

65600 

1 

CaNH»C0NFR(TW0TH,r,FZER03 

65700 

1 

RES»XCU)-CONH 

65B00 


RETURN 

65900 

2 

RESPCXrOCl)*XCU>/THEfR 

66000 


RES*RES/Wl 

66100 


RETURN 

66200 

3 

RES»CS6*(XFDC2)*XC2))/THETR«X(2) 

66300 


RES*RES/«2 

66400 


RETURN 

66500 

4 

RES»CXFDC3)-XC3))/(THETR«DUMMy2)*C534X66*CFZER0/DUM«Y3)*0ASHI 

66600 


RES»RES/W3 

66700 


RETURN 

66800 

5 

RESaC50fX(l)»X(4J"*C51-XCll)*C0UMMY444TNQTH3/ 

66900 


10UMMY3»C524X(53 

670 00 


RES*RES/W4 

67100 


RETURN 

67200 

6 

RES=C54#X(3 3/DUMMY2-(X(S3 3-C554(rCl)4DUHMY4/DUWY3)* 

67300 


1X66-C556’»'X(45 

67400 


RESaRES/WS 

67500 


return 

67600 

205 

IF (K,GT,73 GO TO 206 

67700 


REv5BC7*XC2 3 4XC43/CXC4)tCXCll5*(DUMMY4*»‘4TWOTH)/ 

67800 


lOUMMY3 5 3 + CFF0ZER-PZERO5/’CTHEfR»DUMMy33 

67900 

ccc 

RESaRES/W7 

68000 


RETURN 

68100 

206 

IFCK,GT,8)G0 TO 207 

68200 


RES»C7»XC2)4XC43/CXC43TCXCll3*CDUMMY4»fTWOTH)/ 

68300 


inUMMY3)3+CT14DA5HI!<‘N4X664PZER0/CDUMMY34(l*X66)3 + 

68400 


lCFrDCN)-rCN5)*(DUMMY4»4N)/CTHETR»DUMMy33 

68500 

CCC 

RES»RES/F0LD(N5 

68600 


RETURN 

68700 

207 

HES«C74XC23*X(43/CX(43 + (XC113>t«(DUMMY4»»TWOTH)/ 



68800 
68900 
69000 
69100 
69200 
69300 
69400 
69500 
69600 
69700 
69800 
69900 
70000 
70100 
70200 
70300 
70400 
70500 
70600 
70700 
70800 
70900 
71000 
71100 
ft 200 
71300 
71400 
71500 
71600 
71700 
71800 
71900 
72000 
72100 
72200 
72300 
72^00 
72500 
72600 
72700 
72800 
72900 
73000 
73100 
73200 
73300 
73400 
73SOO 
73600 
73700 


lDl|MMY3n+CTl*DASHIfN*X66*FCN*l)*(0[JMMy4*»(PI-l))/CDUMMy3* 
lCl->X66))4CFPD(N5«FCN))*CDUMMY4f»N)/CTHi:TR*DUMKy3 3 
CCC RKv5««RES/PO0DCN) 

RETURN ■ 

END 

C 

c 

SUBROUTINE GOXi0CXOl,XO2,C,XMAX,O,X7,SAl,C65, DUMMY#, ONETHJ 

C ' ■ : ■ 

C THIS SUBROUTINE SODVES EQN,4,2l BY GODDEN SECTION METHOD 

C FOR OPTIMIZATION, 

DOUBLE PRECISION X01,X02,C,XMAX, D,X7,SAI,C65, Cl, ONETH 
DOUBLE PRECISION X3,X4,TEMPX2, 13, y4, XI, X2 
Ci»C/10,0 

irCABSCX0l»X02),LT.C) GO TO 100 
/ Xl»XOi 
X2« X02 

X3»Xlt0,3820*CX2»X13 

X4»X2-0,3820*(X2*X13 

CALL 0BaECTCT3,X3,D,X7,SAI,C65,DUMMy4,0NETH) 

CALL 0BaECTCY4,X4,D,X7,SAl,C65,DUMMY4,0NETH) 

50 IFCY3,GE,y43 GO TO 30 

Xl*X3 

IFCABSCX1«X23,LT,C3 GO TO 100 
TEMPXTeXUO.asaO^CXT-Xl) 

IFCABS(TEMPX2-X43,GT, C13 GO TO 40 
• y3ay4 
GO TO 45 

40 X3aTBMPX2 

CALL 0BJECTCY3,X3,D,X7,SAI,C65, DUMMY#, 0NEfH3 
45 X4aX2-0,38204tX2-Xl) 

CALL 08aECTCy4,X4,D,X7,5AI,C65, DUMMY#, 0NETH1 
GO TO 50 

30 X2=X4 

IK(ABSCXl-X2) ,LT,C3 GO TO 100 
TEMPX2aX2«0,3820»(X2-X13 
IF(ARSCTF.MPX2-X33,GT,C13 GO TO 60 
Y43Y3 
GO TO 65 

60 X4a TEMPX2 

CALL 0BJECT(Y4,X4,D,X7,SAI,C6S, DUMMY#, ONETH) 

65 X3aXl+0,3820 *CX2-X13 

CALL nBJECTCY3,X3,0,X7,SAI,C65, DUMMY#, 0NETH3 
GO TO 50 

100 XMAXaexil-X23/2,0 

CALL DHUECTCYMAX, XMAX,D,X7,SAI, C65, DUMMY#, QNETH3 
YMAXXaYMAX 



RETURN 

END 


73800 
73900 
74000 C 
74100 C 
74200 C 
74300 r 

74400 SUBROUTINE OBJBCTCTOPf ,XOPT,D,X7,SAI,CS5,0UMMy4,ONETH5 

74600 C 

74700 C THIS SUBROUTINE REPRESENTS THE OBJECTIVE FUNCTION TO BE 

74B00 C OPTIMIZED FOR THE SOLUTION OF EQN, 4,21. 

74900 C 

75000 C’*'*44****itt44»4***»***»4***»*«i|5f*#f»*»»****4»4»»<)|ef »♦♦»♦♦♦♦»♦»♦♦*♦♦*»♦»#♦♦ 

75100 ,e 

75200 DDUBI.E PRECISION yOPT,XOPT,0,X7,SAI,C65,ONETH 

75300 yOPT«l,0-XOPT4DIiOG(XOPT)+SA14Cl,0«XOPf>4Cl,O-.XOPT)+C65*X7/ 

75400 lCD*CDUMMy4T*0NETH)) 

75500 yOPT»wABS(YOPT) 

75600 RETURN 

75700 END 
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UATA F0»» THE PRHHLEM 

data FHR cilNSYRATNED OPYtMlZfl'yin M^-l / * 

K s 14 ITMAX 5= 2000 IPRIN s o 

BETA * 0.000100 DEliTA a 0,001000 

BASHI a 0,900000 


NOPT all M a 16 

ALPHA a 1.300000 
GA^K A a S 


TNITTAL GUESSES, 


XSTARK n 
XSTARTC 3) 


0,40000000E+01 

0,16300000B+01 


XSTARTf R) a 0*99000000E+00 
XSTARTC 7) a 0. 99000000E+00 
XSTAPTC 9) a 0,30000000E+01 
XSTAPTCIU a 0.24000000E+02 


XSTARTC 2) 
XSTARTC 4) 


0,99000000E+00 

0,99000000£400 


XSTARTC 63 a 0 ,64000000£+00 


XSTARTC S3 
XSTARTC 103 
XSTARTC 


0.14000000E+Ot 

0,60000000£+01 


INPUT parameters EOR LINEAR EQUATIONS AND THE CONVERGENCE CRITERION 


MAXIT a 200 
NUMEO a 6 


NUMSIG 


EPS a 0,IOOOOOOOE-07 
ZETA a O.lOOOOOOOE-08 


IPT a I 


DELNON a 0,10000000E-06 PREC 


0.50000000E-05 


rate constants^—-. 

Kp a 0,26500nOOE+06 
KIP a 0,OOOOOOOOE400 
KFT a O,000O0O00E400 


KD a 0.23300000E-04 
KTD a O.OOOOOOOOE+00 
KFM a 0.21200000E+02 


WUNDMfe;R^PO|,jMER CHARECTARISTtCS-.. 


vmol 


0. t040n000E + 03 
0.1 1R45000E403 


REACTOR CHARFCTARISTTCS. 


0,10526000E+01 

0.25600o6oE+00 


KZERO a O.OOOOOOOOE+00 
KTC a 0,72000000E+10 


DM a 0,87800000E+00 


TriFTA 


0,9000UOOOE+03 


TEMP a 0.34300000E+03 


VW s 0.64940000E+00 


other nuMSTANTS»„„ 

Ainw a 0,1S400000E-04 
VM » O,6544OSO0E-l 9 
FEE » 0,10000000E+01 


AMEERO a 0.16fillOOOE-12 VD a 0 ,41888000E-U 
CE a 0.367200001+10 AMM a 0 .63000000E-02 
SWC a 0,89000000E-05 YETA a 0 ,75000000E+01 


feed conditions. 


RADICA?' number DTSTHTBUTION 


FFDC n * 0,000000 

rroc 31 » 0.000000 

Ef'DC 53 a O.OOOOOO 

FFDC 71 a 0,000000 


FFDC 2) a 0,000000 
FFOC 43 a 0,000000 
FFDC 6) a 0,000000 
FFDC ft) a 0,000000 



T.IVE priLYMJrR chains, 


FDLIVf 

D, 1) * 

0,0000 

FDLIVf 0< 

FPLTVC 

f), 4) « 

0.0000 

FDLIVf 0, 

FPLTVC 

0, 71 « 

0,0000 

FDLIVf 0, 

fOLlVC 

0,10) * 

0,0000 

FDLIVf 

FHLTVf 

t » 1 ) * 

0.0000 

FDLIVf t, 

FDLTVC 

1, 4) s 

0.0000 

FDLIVf I, 

FPLIVf 

1, 7) * 

0,0000 

FDLIVf 1. 

FHLTVf 

1.10) ■ 

0,0000 

FDLIVf 

FPLTVC 

2, 1) « 

0,0000 

FDLIVf 2. 

FOLTVf 

2 , 4 ) • 

0,0000 

FDLIVf 2, 

FHLTVr 

2, 7) * 

o.oooo 

FDLIVf 2, 

FDLTV( 

2,10) * 

0,0000 

FDLIVf 

FDLIVf 

3, 1) a 

0,0000 

FDLIVf 3, 

FDLIVf 

3 , 4 ) ■ 

0,0000 

FDLIVf 3. 

FDLIVf 

3, 7) « 

O.UOOO 

FDLIVf 3. 

FDLIVf 

3,10) a 

0,0000 

FDLIVf 

FDLIVf 

4, 1) ■ 

0,0000 

FDLIVf 4. 

FDLIVf 

4, 4) * 

0.0000 

FDLIVf 4. 

FDLIVf 

4, 7) « 

0,0000 

FDLIVf 4, 

FDLIVf 

4.10) * 

0,0000 

FDLIVf 

FDLIVf 

5# 1) » 

0,0000 

FDLIVf 5, 

FDLIVf 

5. 4) * 

0,0000 

FDLIVf 5, 

FDLIVf 

5, 7) * 

0,0000 

FDLIVf 5. 

FDLIVf 

5.10) » 

0,0000 

FDLIVf 


r^ERD polymer chains 


FD0FO f 

0, 1) 

s 

0,0000 

FDDEDf 0, 

FDDFLf 

0, 4) 

s 

0,0000 

FDDEDf 0, 

FDDFDf 

0, 7) 


0,0000 

FDDEDf 0, 

FDDFDf 

0.10) 

s 

0,0000 

FDDFDf 

FDDFDf 

1. U 

a 

0,0000 

FDDEDf 1, 

FDDFDf 

1. 41 

35 

0,0000 

FDDEDf 1, 

FDDFDf 

1,7) 

» 

0,0000 

FDDEDf 1, 

FDDFDf 

1,101 

9 

0,0000 

FDDEDf 

FDDFDf 

2. 11 

S5 

0,0000 

FDDEDf 2, 

FDDFDf 

2, 41 

S 

0,0000 

FDDEDf 2, 

FDDFDf 

2, 71 

« 

0,0000 

FDDEDf 2, 

FDDFDf 

2,101 

9 

0,0000 

FDDEDf 

FDDFDf 

3, 11 

m 

0,0000 

FDDFDf 3, 

FDDFDf 

3, 41 

a 

0,0000 

FDDFDf 3, 

FDDFDf 

3, 71 

a 

U ,0000 

FDDEDf 3, 

FDDFDf 

3,101 

a 

o.oooo 

FDDEDf 

FDDFDf 

4, n 

a 

U.UOOO 

FDDEDf 4, 

FDDFDf 

4, 41 

M 

0 .oooo 

FDDEDf 4, 

FDDFDf 

4, 7 ) 

9 

0,0000 

FDDFDf 4, 

FDDFuf 

4,101 

9 

0,0000 

FDDEDf 

,^,,FDDFDf 

5, n 

• 

0,0000 

FODEDC 5^ 


2) » 

0,0000 

FDLIVf 

0. 

31 

s 

0,0000 

5) *■ 

0^0000 

FDLIVf 

0. 

6) 

a 

0,0000 

83 ■ 

0,0000 

FDLIVf 

0. 

9) 

3 

0,0000 

2) • 

o '. 0000 

FDLIVf 

1 . 

3) 


0,0000 

5) a 

0,0000 

FDLIVf 

1 . 

6) 

3 

0,0000 

8) a 

o.oooo 

FDLIVf 

1. 

91 

3 

0,0000 

2) a 

0,0000 

FDLIVf 

2. 

3) 

3 

0,0000 

S ) SB 

0.0000 

FDLIVf 

2 . 

6) 

S 

0.0000 

8) a 

0,0000 

FDLIVf 

2 . 

91 

3 

0,0000 

2) a 

0.0000 

FDLIVf 

3, 

3) 

S 

0,0000 

5 ) a 

0,0000 

FDLIVf 

3. 

6) 

3 

0,0000 

8) a 

0.0000 

FDLIVf 

3. 

9) 

3 

0,0000 

2) a 

O.OOOD 

FDLIVf 

4. 

3) 

3 

0.0000 

5) a 

0.0000 

FDLIVf 

4. 

6) 

3 

0.0000 

8) a 

0,0000 

FDLIVf 

4. 

9) 

S 

0,0000 

2) a 

0,0000 

FDLIVf 

5 . 

3) 

S 

0,0000 

5) a 

0,0000 

FDLIVf 

5. 

6) 

3 

0,0000 

8) a 

0.0000 

FDLIVf 

5 . 

9) 

3 

0,0000 


2) 

3 

0,0000 

FDDEDf 

0. 3) 

3 

0,0000 

5) 

3 

0.0000 

FDDEDf 

0. 6) 

3 

0.0000 

81 

mm 

0,0000 

FDDEDf 

0. 9) 

3 

0,0000 

2) 

mil 

0,0000 

FDDEDf 

1. 3) 


0.0000 

5) 

3 

0.0000 

FDDEDf 

1 , 61 

3 

0.0000 

8) 

3 

0,0000 

FDDEDf 

1. 9) 

3 

0,0000 

2) 

3 

0.0000 

FDDEDf 

2, 3) 

8 

0,0000 

5) 

3 

0,0000 

FDDEDf 

2. 6) 

S 

0,0000 

8) 

3 

0.0000 

FDDEDf 

2, 9) 

a. 

0.0000 

21 

3 

o.oooo 

FDDEDf 

3, 3) 

3 

0,0000 

5) 

3 

0,0000 

FDDEDf 

3. 6) 

8 

0.0000 

8) 

m 

o.oooo 

FDDEDf 

3, 91 

3 ’ 

0,0000 

21 

3 

0,0000 

FDDEDf 

4. 3) 


0,0000 

5) 

3 

0.0000 

FDDEDf 

4, 6) 

3 

0,0000 

8) 

3 

0.0000 

FDDEDf 

4, 9) 

3 

0.0000 

2) 

m 

O.unoo 

FDDEDf 

5. 3) 

3. 

0.0000 




5. 4) a 0,0000 


FDDEDf S, 5) a 

0^0000 

FD0BD( 5, 

FHOFDf 

77 * 0,0000 


FODEDC 5, 8) a 

0.0000 

FODEDf 5, 

FHOFOr 

‘'.10) a 0,0000 


FDOEOt 



FEFO 

cnwD. priR VftRTABOES 

OF INITIAL SST 



1 

0.4080E+01 

2 

O.iOOOE+Ot 

J 

0.1920E+01 

4 

O.OOOOE'fOO 

b 

O.OOOOEfOO 

6 

0,0000E+00 

7 

0,OOOOE+0O 

a 

0.OO0OE4OO 

9 

0,0000B+0O 

10 

Q.0000&>00 

li . 

0,0000E+00 





0,0000 

0.0000 


UtITPUT PROM the PROCRaM l_^:.M.:.,..:,.^^ • * ' ^ “ 

$$.*********^*****$'***$*****$****1^$*3t^*^*$*!$'$*!$$!^***t**’¥**>$**)¥*$*^*$***tt**** 

(1) NHMRER hr PARTICLES (/CC EMULSIONJ^^ 0. 1 1466412E+16 

(2) average VOLtmE nr particles (CC),-«-_ 0.26767$43E-15 

(3) average DIA iVOLiJME BASIS r*a 7, ] 0,79956I24 Bt03 

f4) average DIA iNUMBER BASIS f "A) — I— 0.71403866Et03 

(5) CnwVERSION ORTATiapn. _ * 0.40438988E+00 

C6) RATE (MOLES/ (LIT. HOUR) O.ATSaSOOOE-^Ol 

(/) MUMBER AVERAGE MOLECULAR 0.39254813E+06 

(0) WEIGHT AVERAGE MOLECULAR WEIGHT^.,^^ 0.921S7680E+06 
fy) POLYDISPERSTIY I R.23476785E+01 











DETATL.*! HF the DTSTRTBITTIONS F0LLQW_1_'_1 ‘ 



T,1VE 

pni.YMFR CHAIN 

IiRngth distribution 




Mhf,. «T, 

NUMB. BASIS 

WT.BASISCHORMAI.) 


0 

o.ionoE+oi 

0.U34E-03 

0. 60371-07 


1 

0,7010^403 

0.tll6E-03 

0.41641-04 

2 

O.i40ir+04 

0.941 2E-04 

0.70201-04 

( 

3 

0.210,IE404 ■ 

0.7868E-D4 

0,88011-04 


4 ' 

0.280tF*^04 

O.6S48E“04 

0,97651-04 


S 

o,isoir^n4 

O.S434R-04 

0.10131-03 


6 

0,42'"ir+O4 

O.4S01E-O4 

0,10071-03 



0.A90|r^ft4 

0,3722E-fl4 

0,97121*-O4 


8 

0,i::ftO),E^04 

0.3074^*04 

0,91681-04 

f, 

9 

O,630ir«-O4 

0.2537E-04 

0,85111-04 

ll. 

10 

0,7of»lE*04 

0.2092E-04 

0,77991-04 

i§. 

11 

0.770!Ff04 ■ 

0.t725E*O4 

0.70711-04 : 

II 

12 

‘0.8401^*04 

0,1421E"O4 

0,63561—74 

It 

13 


0.1 171E-04 

0.56731-04 

p;; 

14 


0,0646E~05 

0,503 3'?'-04 


15 

0,i050Et05 

0,7947R*05 

0,44431-04 

te 

1 fi 

O^lUOFtOS 

0,6549E-O5 

0.39061-04 

p 

17 

O.J120E4>05 

0.S399E-05 

0,34211-04 


S 18 

0*l26ORtO5 

0.4452E-05 

0.29871-04 



0,U30E+05 

0.3b73E-05 

0.26011-04 


1:20 

0.1400E+05 

0.3032E-05 

0,2260E-04 


' 21 

0.147OE+O5 

0.2504E-05 

0.19601-04 


22 

0.l540Et05 

0,2069E-05 

0,1697E-04 


23 

0,1610E+05 

0.171 1E»05 

0,14671-04 


24 

0.16R0E+05 

0.1416E-05 

0.12661-04 


25 

0.1750F+05 

0.ll72E-0'5 

0,10921-04 


26 

0,t820E+05 

0.9712E-06 

0.94121-05 


27 

0.1890E+0S 

0,9052E-06 

0.81031-05 


28 

0,1960E-»-05 

0 ,6680E">06 

0.69711-05 


29 

0.2030E+05 

O,5545E-06 

0.59931-05 


30 

0,2100F+05 

0.4605E-n6 

0.51481-05 


31 

O,2l70F+05 

0,3826E-06 

0.44201-05 


32 

0.2240F+05 

0,3l80E-06 

0.37931-05 


33 

0.7310Et05 

0.2644E-O6 

0.32521-05 


34 

O.?380E+0b 

0,2l99E-06 

0,27871-05 


35 

0.2450Ft05 

0.18291-06 

0.23861-05 


3t, 

0.2520F+05 

0,15221-06 

0,20421-05 


• / 

O,2500F+05 

0,12661-06 

0.17461-05 


*4 8 

0,'?6«‘'OF4 05 

0,10531-06 

0.14911-05 


'‘4 V 

f'. >7 

0.87581-07 

0.12731-05 



•I . 

0,72811-07 

0.1 0851-05 


1 

ft ’;87aFi-0S 

o,6050«'-0? 

0.9245F-06 




0.50251-0? 

0.78651-06 


- „l 

^C'1 

0,41701-01 

0.66831-*^6 



^C'^'rr,v''3 

0. 34581-0 s 

0.567lW“,.0fc 


A 

t\. !.J, 5 iVj 

SU 28651-0'’ 

C,48O51«-06 


i 

ft *'•■'- 

0,/237 11-01 

0,. 40641 -06 


^ i 

ft _ ft'ft 

ft , ^yi^vF-O'’ 

0.34321-06 


< , 

'V Mf.ri*<-0 5 

ft . 1 61 71“0 ’ 

0. 28931-06 


%•! . ' 

ft , >41 ft - '■' i 

0.1332V>ft^ 

0.24331-0® 


r. s 


ft . \ 0 961-0 ? 

C.2042S'-0®. 


r*/;flD Pf'i?,rKKk C’'';ATN LKNCTH siISTRlBUTlOf« 


Wt 1 



WT. BASIS (NORMAL 

, > 


0„<l5798+f)0 

0,42178-06 

1 

0.389a«:403 

0,71958-01 

0,40608-04 

2 

0. tl77R + 04 

0,5342E«01 

0.5790E-04 

J 

0.1165P;t04 

0,4136E-01 

0.67238-04 

4 

0.2353K+04 

0.3326E-01 

0,72068*04 

5 

0,2941F+04 

O,2739E«0l 

0,74178-04 

6 

0,3529Rt04 

0,2295E«01 

0.74608*04 

7 

0,4ll7R+04 

0.19S1E-01 

0.73968*04 

8 

0,470SR+04 

O,1677E-0l 

0,72678*04 

9 

0,5293Rt04 

0,1456E-01 

0,70958*04 

10 

0,58R1R+04 

0,1273Et0l 

0,68968*04 

11 

0.6469R404 

0.1121E»01 

0.66808*04 

12 

0.70S7Et04 

0*9928E*02 

0,64528-04 

13 

0,7645Rt04 

0,8828E-02 

0.62158-04 

1 4 

0.82331:404 

0,7877E-02 

0,59728*04 

15 

0,882lBt04 

0,7046E-02 

0,57238-04 

U 

O,94OtB404 

0,631SE-02 

0.54718*04 

17 

0.9997!!;404 

0,56668*02 

0,52168-04 

18 

0a089B+O5 

0.5087E*02 

0.49598-04 

19 

0,11178405 

0. 45691*02 

0,470lE*Q4 

20 

0.1t76e405 

0,4102E*02 

0,44438*04 

21 

0.1235E405 

0,3682E-02 

0.4187E-04 

22 

0.1294R405 

0 , 3301E*02 

0.39338-04 

23 

0.1352R405 

0,2957E*02 

0.36828-04 

24 

0.14tlR405 

0.2645E*02 

0,3437E-04 

25 

0.1470R+0b 

0.2362E-02 

0,31988-04 

26 

0.1529R405 

0.2106E*02 

0,29658-04 

27 

0.15888405 

0,l874E-02 

0,27408-04 

28 

0,1647R405 

0.1664E-02 

0.2523E-04 

29 

0,1 705R+05 

0,1475E-02 

0.23168-04 

30 

0.1764R405 

0.13048*02 

0.21188-04 

11 

0.18238405 

0.1 150E-02 

0,19308-04 

12 

0,18828+05 

0,1011E*02 

0,17528-04 

33 

0,1 9418405 

0,8870E-03 

0. 15858-04 

34 

0,1989R+05 

0,77578*03 

0,14288-04 

3b 

0.20588+05 

0,67 628*03 

0,12828-04 

36 

0,21178+05 

0.5B75E*03 

0,11458-04 

17 

0,21768+05 

0. 50868-03 

0,10198-04 

38 

0,2234R+05 

0,43858*03 

0,90248-05 

39 

0 , 229 38 +05 

0,37668*03 

0,79528-05 

40 

0,23528+05 

0,32188-03 

0,69718-05 

4 1 

0,241 18+05 

0,2737E-03 

0,60768-05 

42 

0,24708+05 

0,23148*03 

0.52638-05 

4,3 

0,25298+05 

0,1944E-03 

0.45278-05 

44 

0,25878+05 

0,tb22E-n3 

0.38658-05 

45 

0,26468+05 

0.13428*03 

0.32708-05 

46 

0,27058+05 

0,11008-03 

0.27398-05 

47 

0,27648+05 

0,89068-04 

0,22678-05 

48 

0.28238+05 

0,71168-04 

0,18508-05 

49 

0.28818+05 

0.55868.-04 

0,14828-05 

50 

0,29408+05 

0,42878*04 

0.11618-05 



PAPTICr.fc. ST 2 K DIST»1WUT10|« tCUMULAf IVEl 


NO 

DTA( "A) 

FRACTION 

0 

0.8000r+02 

0,0000E+00 

1 

0,4S83Kf03 

0.1683E+00 

2 

0,5735F:+03 

0.3083E+00 

3 

O.A5fi4!>:+03 

0.4248E+00 

4 

n.7225r+03 

0.S216E+00 

5 

0,77B2K+03 

0,6022E+00 

6 

0,R270Ktn3 

0,6§92E+00 

/ 

0,fi706ft03 

0,7249E+00 

a 

0,9102f+03 

0,77l3E+00 

9 

0*9486f+03 

O,8a98E+0O 

10 

0,9a0SRt03 

0*84i9E+00 

11 

0.10!2Et04 

0,8685E+00 

12 

O,lO42Rt04 

0*89a7E+00 

13 

0.1070Ct04 

0.909tE+00 

14 

, 0a097E-*‘04 

0,9244E+00 

IS 

o.iiaartoi 

0.9I73E+00 

16 

0,1147E404 

0.f478E+00 

i1 

0.1l7O|!tO4 

O,9866E+O0 

la 

O.U93e+04 

0,96l9E+00 

19 

0,12t4IEt04 

O,970OE+00 

20 

0,J235E4-04 

O.975lE+a0 

21 

0.1256Et04 

0,9793E+00 

22 

0,I2751*:+04 

0,9828E+00 

23 

0.1294E+04 

O.98S7E+0O 

24 

0.13l3Et04 

0,9881E+00 

25 

0,133lK+04 

O.99O1E+O0 

26 

0.134aEt04 

0,9918E+00 

77 

0.1365P+04 

0,9932E+00 

?a 

0.13R2E+04 

0,9944E+00 

?9 

0.13P8EtA4 

0,9953E+00 

30 

0,1414^+04 

0,996iE+00 

31 

0.1430F+04 

0,9968E+00 

32 

0,1445F+04 

0,9974E+00 

33 

0.1460E+04 

0,9978E+00 

34 

0,1474F+04 

0,9982E+00 

35 

0,14fi9E+04 

0,9985E+00 

36 

0.1503R+n4 

0,99a8E+00 

37 

0.1516E+04 

0,9990E+00 

3a 

0. 1530E+04 

0,9992E+00 

39 

0.1543F+04 

0.9993E+00 

40 

0,1586P+04 

0,y998E+00 

41 

0.156yE+04 

0,9996E+00 

42 

0,15a2R+04 

0.9997E+00 

43 

0.1504F, + 04 

0.9997E+00 

44 

0. 1607E+04 

0.999flE+00 

4S 

0.1619F+04 

0,9998E+00 

46 

0,1631K+04 

0.9999E+00 

47 

0.1642F+04 

0,9999E+00 

4a 

0.1684F+04 

0,IOOOE+Ql 

49 

0,l6fi5F+04 

0,lO0OE+0l 

80 

0.1676F+04 

O,10OOE+O1 
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